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ABSTRACT OF THE DISSERTATION 



Random Delaunay Triangulations, the Thurston- Andreev Theorem, 

and Metric Uniformization 

by 

Gregory Leibon 
Doctor of Philosophy in Mathematics 

University of Cahfornia San Diego, 1999 
Peter Doyle, Chair 

In this thesis a connection between the world of discrete and continuous conformal ge- 
ometry is explored. The world of discrete conformal geometry is related to disk pattern 
construction and triangulation production. In particular we discover a generalization of 
the Thurston- Andreev theorem to any angles in (0,7r] (when x(M) < 0).. The proof of 
this theorem relies on an energy measuring how "uniform" the angle data of a triangu- 
lation is, where by uniform angle data I mean the data of a geodesic triangulation of a 
hyperbolic surface. The connection to continuous geometry is via averaging this energy 
over geodesic triangulations, a process which forms an energy measuring how "uniform" 
(constant curvature) a metric is. In fact, the entire discrete energy proof carries over to 
produce a proof of the metric uniformization theorem of conformal geometry. The ran- 
dom triangulations used to produce this averaging are random Delaunay triangulations, 
and they are explored in some detail. In particular the averaging techniques developed to 
construct the energy on metrics can be used to explore other geometric issues, including 
the production of a probabilistic interpretation of the determinant of the Laplacian and 
a new probabilistic proof of the Gauss-Bonnet theorem. 
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Chapter 1 



Introduction 



It has been known for a while now that the world of disk patterns is intimately 
connected to the world of geometric uniformization. To orient ourselves to what this 
statement means it is useful to mention the most famous case. The disk pattern is the 
Koebe theorem which guarantees that a specified (reasonably nice) graph embedded in 
the closed unit disk can be realized as the nerve of a disk pattern, which is in fact unique 
up to Mobius equivalence. In [^] Thurston presented the idea of (nicely) packing a 
bounded simply connected planar domain, recording the nerve of this packing, and then 
using a suitably normalized solution to the Koebe disk pattern with respect to this nerve 
in order to approximate the conformal mapping of this region to the disk. That such a 
mapping exits is our first example of a geometric uniformization problem, the Riemann 
mapping theorem. Thurston's procedure has been verified to work in several cases. In 
fact this conjecture has been solved in several cases. The case of using a finer and 
finer hexagonal disk pattern for the approximation is very well understood; it was first 



shown to uniformly approximate the Riemann mapping in |16| and is now seen in |10| 
to approximate uniformly all of the Riemann mapping's derivatives. 

In this thesis another connection between these two stories is developed. For 
now I'll give a brief outline which can be viewed as an introduction to the introduction. 
Precise statement of the below discussion can be found in the elaborated introduction, 
sections 1.1 - 1.3. 

In this thesis a generalization of a disk pattern theorem often called the Thurston- 
Andreev theorem is presented, which is the solution to a disk pattern problem on a 
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compact surface. To be precise, in chapter two of this thesis we prove... 

Result 1 Given a geodesic triangular decomposition of a hyperbolic surface one can 
associate combinatorial data: the topological class of the triangular decomposition and 
to each edge the angle between the circumscribing circles of the two triangles sharing the 
edge. If the angle data is in (0, vr] then this combinatorial data uniquely determines both 
the hyperbolic structure and the geodesic triangulation. Furthermore, if one starts with 
a topological triangular decomposition of a compact surface with xi^) < 0; ^'^'^ angles 
with values in (0, vr] assigned to each edge satisfying certain necessary linear conditions 
(conditions (ni) - (77,4) of section 2. 3.J\ ), then there is a hyperbolic surface and geodesic 



triangular decomposition realizing this data. 



The angles in (0, ^] case is well known, see [17|. This theorem has several inter- 
pretations in terms of not only a disk pattern problem, but also as a convex hyperbolic 
polyhedron production theorem, and perhaps most importantly (from the uniformiza- 
tion point of view) as a geodesic triangulation production theorem. This triangulation 
production can be viewed as the second step in the proof of the above result. The first 
step is a completely linear max-flow-min cut-type problem guaranteeing that given such 
edge angle data one can produce triangle angle data which "formally" has the same edge 
angle data associated to it. Given such triangle angle data one can measure how uniform 
this data is by using an energy on the set of possible conformally equivalent angle data. 
This energy is larger when one is closer to the angle data of a geodesic triangulation of 
a hyperbolic surface (so for aesthetic reasons I've chosen an energy which one tries to 
maximize, my apologies to the physicists). The proof becomes to guarantee the existence 
of a unique set of angle data maximizing this energy and to show it corresponds to the 
unique solution of the above disk pattern result. 

The geometric uniformization theorem which this result will be shown inti- 
mately related to, is the uniformization for geometric surfaces, namely: 

Result 2 Every Riemannian surface is conformally equivalent to a Riemannian surface 
with constant Gaussian curvature. 



This theorem is very well known. The heart of what takes place here is that in 
chapter 4 there is a proof of the above uniformization theorem which is directly related 
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to the above disk pattern theorem and its proof. This relationship is to first (as above) 
solve the completely linear problem of finding a varying negative curvature metric, then 
we can measure how distant this metric is from a constant curvature metric by averaging 
the energy associated to data provided by randomly selected triangulations. In the end 
one finds the same arguments used to uniformize the discrete angle data can be applied to 
this averaged energy on the space of metrics, providing a nice proof of the uniformization 
theorem of surfaces. This is not the first energy method proof of this theorem, in fact 
there are many. One particularly interesting energy is the det(Ag) as used in Q. In fact 
these random energies are related to the det{Ag), and in section 1.3 this relationship is 
made explicit. 

The main tool in articulating the averaging procedure is a detailed understand- 
ing of random Delaunay triangulations. The details of this understanding form chapter 
3. I think these triangulations are quite interesting independent of their connection to 
metric uniformization, and I hope to convince others of this by beginning the precise 
introduction with a probabilistic proof on the Gauss-Bonnet theorem. 

It is worth noting that these ideas may be explored in three-dimensions. In fact 
much of what takes place here can be viewed as realization in two dimensions of an idea 
by Peter Doyle to prove that a three manifold accepting a varying negative curvature 
metric also accepts a constant sectional curvature metric. If proved this result would 
shed considerable light on Thurston's uniformization conjecture. In section 1.3 I sketch 
the three dimensional strategy, as related to the two dimensional strategy presented in 
1.1-1.2. 



1.1 A New Proof of the Gauss-Bonnet Theorem 

The Gauss-Bonnet formula in its most primitive form is stated for a compact 
boundaryless surface endowed with a geometry as: 



^ [ kdA = x{M). 

Jm 



The quantity, k , being integrated over the surface, M, is the Gaussian curvature, and on 
the right-hand side is the Euler Characteristic. The Euler characteristic is a topological 
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invariant of the surface which can be computed from triangulation as 

X(M) = V-E + F, 

where V is the number of vertices, E is the number of edges, and F is the number 
of faces in the triangulation. In a somewhat unintelhgible nut shell, the proof here is 
accomplished by randomly triangulating the surface and then noting while the Euler 
characteristic is constant the other (now random variables) E, F, and V have expected 
values which can be computed and compared to this constant. As the density of the 
randomly distributed vertices goes to infinity one finds these expected values produce 
the Gauss-Bonnet formula, along with a probabilistic interpretation of curvature. 

This proof's initial need is an articulation of what a random triangulation of 
a surface is. The first step is to ignore the fact that anything random is going on here 
and to simply attempt to construct a a geodesic triangulation in a fixed metric g from 
a given set of points, p = {pi, To do this one examines all the triples and pairs 

in {pi, ...,Pn} and decides whether or not to put in a face for a given triple or an edge 
for a given pair; this decision procedure will be relative to a certain positive number 6 - 
the decision radius. The procedure (with its origin found in Delaunay's "empty sphere" 
method, see |5|) is to formally put in a face for a triple or an edge for a pair if the triple 
or pair lies on a disk of radius less 6 which has its interior empty of points in {pi, ...,pn}. 
Clearly every pair on a face forms an edge (though not necessarily the converse), hence 
at this moment we have an abstract two complex. Call this abstract complex K^p^^ p^^y 
and denote a polyhedral realization as |-f^{pi,...,p„}| (for a realization of it place the points 
at (0, . . . , 0) and (0, . . . , 0, 1, 0, . . . , 0) in R""-*- and take the convex hull of the subsets 
found in -fC{p,,...,p„}). 

Our goal now is to realize this polyhedron inside the surface, which will involve 
restricting 5. To get going we must recall a well known geometric constant associated to 
the surface, the injectivity radius, which is the largest number i such that if d{p, q) < i 
then there is a unique geodesic between p and q of length less than i. If one assumes 
that 5 < g one can (see the proof of ||) produce a continuous map 

R-\K{p,_p„}\^M, 

with the edges parameterizing geodesies in the g metric. That our procedure forms a 



Figure 1.1: Part of a Typical Delaunay Triangulation 



triangulation is now equivalent to R being a homeomorphism; and in such a case we will 
call the resulting triangulation a Delaunay triangulation. Such triangulations typically 
will look locally like figure 1.1 . 

At this point it is useful to find a geometric criterion on a set of points {pi, p„} 
guaranteeing that it forms a triangulation. To accomplish this it is useful to recall a sec- 
ond geometric constant associated to the surface: the strong convexity radius r. On a 
compact surface it is the largest number with the property that in a disk of radius less 
than T the interior of the unique minimal length geodesic connecting any pair of points 
in the disk's closure is in the disk's interior. We may now introduce our criterion for 
triangulation detection: we will call a set of points {pi, ■.■,Pn} (5-dense if each open ball 
of radius S contains at least one pi and {pi, ■■■,Pn} contains no four of its points is on a 
circle of radius less than mm{|,r}. With this concept we will prove in section 2.2 the 
following theorem... 

Theorem 1 If {pi, ...,pn} is 5-dense with 6 < mm{|,T}, then {pi,...,pn} forms a 
Delaunay triangulation. 

From here on out always assume the decision radius 6 satisfies 6 < 

min{|, r}. 

Now comes the second step in this proof: to put down random sets of points. 
The points will be distributed with a density A via a Poisson distribution relative to g 
(developed in the section 2.3), and it is convenient to also randomly select an ordering 
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of the points after placing them down. In technical terms, any such point configuration 
will be denoted p and will live in a measure space P with a probability measure 
and specified measurable sets B - all of which are constructed in section 2.3. For now 
all we need to understand are certain basic properties and terminology. Let U C M 
be measurable with area A(U). The first property is the characteristic property of the 
Poisson distribution, namely the set of point configurations in P with exactly n points 
in U is an element of B with measure ("^(^)''^) ^ 

n! 

Recall a random variable on P is simply a B measurable function. We will 
encounter several in this proof, namely E{p), V{p), -F(p) , and the characteristic (in- 
dicator) functions of various measurable sets (denoted lu with U measurable). The 
expected value of a random variable L{p) is the integral of the function over P with re- 
spect to the measure Pa, and will be denoted Ea(L). As an example of such an expected 
value, using A to denote the area on M, we have from the above characteristic property 

n=0 ■ n=l ^ ^' 

which reveals the use of the term density for A. The most important set for us in P is the 
set C P consisting of configurations which have Delaunay triangulations associated to 
them. By theorem |^ when our points are put in with a high density we expect they will 
typically form Delaunay triangulations. Allowing 0(A^°°) to mean a quantity decaying 
faster than any polynomial in A, in section 2.3 we prove this with: 

Theorem 2 If L is any one o/It^-, E, V, or F, then L is measurable and has expected 
value Ea(L) = Ea(L1tJ + 0(A-°°). 

Now one observes that the Euler characteristic of any configuration in is the 
Euler characteristic of an actual triangulation, hence the constant xi^)] so along with 
the above theorem |2| one has 

X(M) = Ea(x(M)1tJ + 0(A-°°) = Ea((V - E + F)1tJ + 0(A-°°). 

Applying the mathematical triviality, yet philosophical miracle, that expected values add 
gives us 

X{M) = Ea(V1tJ - Ea(E1tJ + Ea(F1tJ + 0(A-°°). 
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In an actual triangulation we have = F so 



X{M) = Ea(V1tJ - |ea(F1tJ + Ea(F1t,) + ©(A" 



As explored above, we have the expected number of vertices is XA; this along 
with theorem gives us 



X(M) = - ^Ea(F1tJ + 0(A- 
So finally using theorem |^ one last time we find: 
Formula 1 (Euler-Delaunay-Poisson Formula) 

X(M) = hm (aX-Iex{F: 



X—>oo \ 2 

The goal now becomes to compute Ea(F). Now a triple of points occurring in 
a Poisson configuration, y G M x M x M, forms a face exactly when one of its associated 
disks is empty of points. The first thing we'd like is that we are in fact talking about a 
unique such disk. This is the most sophisticated fact needed, and occurs throughout the 



proof. I'll record it here for future reference; for a couple of proofs see |12]: 



Lemma 1 (The Small Circle Intersection Lemma) // a triple of points lies on a 
the boundary of a disk with radius less than 6, then this disk is unique among disks of 
radius less than 6. (Recall 6 < min{^,T}.) 

Back to the computation, by the characteristic property of the Poisson dis- 
tribution the probability any area a disk is empty is e"^". So letting Vs C x^M be 
the set of ordered triples living on circles of radius less than 6, we now see the proba- 
bility that y has an associated face should be fy^ e~'^"'^y^dA^ , with a{y) the area of 
y's uniquely associated disk. This together with the fact that the expected number of 

triples is Yl'^=o (3) ^'^^''nl ~ ^^e^ leads one to hope, as is confirmed in the beginning 

of section 3, that 



\3 I- 

Ea(F) = — / e-^^(y)dA=^ 
6 Jv, 



(1.1) 
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To compute this explicitly it is necessary to put coordinates on Vs- First one 
chooses a way to discuss directions at all but a finite number of tangent planes of M 
(via a orthonormal frame). Then one can parameterize a full measure subset of Vg with 
a subset of (6*1, ^2, ^3, r,p) £ x x x (0, 6) x M hy hitting the triple described by 
sitting at the point p £ M and moving a distance r in each of the three directions Oi, 
02, and ^3 . That this is a legal parameterization is easily deduced from lemma ^, and 
discussed in section |3.2.2 . Note that when fixing p and varying 9i in these coordinates 



one produces the Jacobi field Ji{r); whose norm I will call ji. Using this notation, letting 
d9 = dOi A (102 A d9^, and letting 1/(9) be the area of the triangle in the Euclidean unit 
circle with vertices at the points corresponding to the {9i}, we will see in section 3.1 that 
equation (1) expressed in these new coordinates is 
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e"^<'''-^PheJeJeA^)dOdrdA. (1.2) 



Both ji and the area of a ball function, a{9,r,p), are directly related to the 



curvature; with this relationship equation 1.2 can be expressed as 

— [ /V e-^^'-\{9) (r^- -r^ + + "©"(r^)^ d9drdA, (1.3) 

6 JmJo 7x35-1 V 2 12 J 

with the "0"(r^) a well controlled function of 9, r, p and A. In fact the "0"(r^) term is 
so controlled that there is a pM > associated to the surface such that if 6 < pM after 



integrating one can reduce equation 1.3 to 



Ea(F) = 2AA - - /" kdA + 0(A-^' 



(1.4) 



with the |0(A"2)| < CA~2 (see section 3.2.3 for the details concerning 1.3 and 1.4 and 
the related estimates). 

In particular we may now plug this into the Euler-Delaunay-Poisson Formula 
(formula |l]) to give simultaneously a probabilistic interpretation of curvature and new 
proof of the Gauss-Bonnet theorem. 
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Formula 2 (Euler-Gauss-Bonnet-Delaunay Formula) Using a decision radius 6 < 
Pm we have 

X(M) = lim (^A - ^E;,(F)) = ^ [ kdA. 

Note this allows us to interpret the Gaussian curvature as the defect in the 
expected number of faces in a random Delaunay triangulation in the surface's geometry 
from what would be expected in Euclidean space. 



1.2 Conformal Geometry and Uniformization 

The techniques used in the previous section to prove the Gauss-Bonnet theorem 
can be used to examine certain aspects of conformal geometry. The ideas will involve the 
comparison of conformally equivalent metrics. In terms of Riemannian metrics we say 
two metrics on M, g and h, are conformally equivalent if h = e^'^g for a smooth function 
We will always be thinking in terms of a fixed background metric called g and will 
label its associated geometric objects like its gradient, Laplacian, curvature, norm, area 
element, or area as V, A, A;,| • | , dA or A. For the h = e^'^g metric we shall denote these 
object with an h subscript. 

The key to using the ideas of the previous section is the development of a dis- 
crete analog of a conformal structure and a conformal transformation relative to a Delau- 
nay triangulation. The most important observations in guessing this discrete structure 
concerns two triangles ti and t2 in a Delaunay triangulation sharing an edge e, and 
the intersection angle 9^ between the circles in which ti and t2 are inscribed (see figure 
|1.2D. The facts of use to us are that 0*^ G (0,vr) and 6'^ is preserved under a conformal 



deformation of the metric up to order 0(r) (see formula ^ in section 4.1.1). 

To use this information it is necessary to note 9^ can be written down accurately 
in terms of the angles within the triangles. 



Formula 3 Letting A^', and C" be the angles in ti as in figure l.i we have 



r = + + O(r^). 

It is worth noting this formula is exactly true on a surface with constant cur- 
vature. One proof of this formula in the negative curvature case utilizes an object, an 
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Figure 1.2: The Notation of Neighboring Triangles 



Figure 1.3: The Ideal Prism 

ideal hyperbolic prism, that will turn out to be fundamental to all that takes place here; 
and it is useful to construct and discuss it now. That it shows up at such a fundamental 
point in understanding the triangulations will hopefully inspire its otherwise unexpected 
appearance in the next section. 

The prism is constructed from the angle data of a hyperbolic triangle, namely 
a set of positive angles {A, B,C} such that A + B + C — it < 0. To construct it first form 
a hyperbolic triangle with the {A,B,C} data, then place the triangle on a hyperbolic 
plane and then place the plane in hyperbolic three space. Now union this triangle with 
the geodesies perpendicular to this 2-plane going through the vertices of the triangle. 



The prism of interest is the convex hull of this arrangement, See figure 1.2. 

With these conformal geometry facts in mind we may set up our discrete con- 
formal geometry. 
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1.2.1 Discrete Uniformization 

The discrete object replacing a metric is the information naturally associated 
to a Delaunay geodesic triangulation living on the surface in this metric; namely a 
topological triangulation, T, and the angles in all the triangles. Given T there are 3F 
slots {oi} in which one can insert possible triangle angles, which we will place an order 
on and identify with a basis of a 3F dimensional real vector space. With this basis choice 
we will denote this vector space as R^^, and denote vectors in it as a; = ^ A*aj. Let a* 
be a dual vector such that Q*(aj) = 6j. Sometimes we will refer to the angles in a fixed 
triangle, and will abuse notation in order to refer to figure |1.2| by letting a^{x) = A^, 
(3'{x) = B\ and = C\ 

We will require that these angles satisfy certain requirements that the angles 
in an actual Delaunay triangulation would satisfy. In particular we require that all the 
angles are in (0, vr) and that the sum of all the angle at a vertex is 27r. 

Motivated from the previous section we define the intersection angle at an edge 
to be the function 

^ {x) = + , 



with the notation inspired from figure 1.2. The Delaunay quality is refiected in the fact 
that the "intersection angles" must satisfy 9'^{x) G (0,7r) for each edge e. It is worth 
noting that as an immediate consequence of the fact that the angles sum to 2tt at a 
vertex we have that the angle discrepancies at a vertex tt — 9^ (x) sums to 2tt at a vertex 



as well, see section 2.3.1 



We will call angles satisfying the above requirements a Delaunay angle system. 
With the previous section as motivation, we define a pair of Delaunay angle systems x 
and y to be conformally equivalent if they arise form the same triangulation and if for 
each edge e the 6'^{x) = 9^{y). 

Motivated by the Gauss-Bonnet theorem we define the curvature of a triangle 
t relative to the Delaunay angle system x to be k^{x) = A + B + C — tt, where {^4, C} 
is the angle data corresponding to angle slots in t. On a surface with < there is 

a very simple set of linear equations which will gaurentee that a Delaunay angle system 
is conformally equivalent to a Delaunay angle system where each triangle has negative 
curvature. 
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Let 5 be a set of triangles in P, and denote the cardinality of S as \S\. Let 
condition (*) on a set of angles {9^{x)} associated to all the edges of P be the condition 
that for any set of triangles S we have 

W Eeesn^)>^\S\- 
As an immediate consequence of theorem ^ we have the following lemma. 

Lemma 2 (The Discrete Teleportation Lemma) A Delaunay angle system x has 
{9^{x)} satisfying (★) if and only if x is conformally equivalent to a negative curvature 
Delaunay angle system. 

We will only be concerned with Delaunay angle systems satisfying the the con- 
ditions of this lemma, which note includes the angle data associated to a Delaunay 
triangulation of surface with varying negative curvature. Denote the convex bounded 
set of negative curvature angle systems conformal to a given one, x, as Nx- 

The pleasure derived from the Delaunay angle systems comes from a beautiful 
Energy which lives on Nx- Let V^{x) denote the volume of the ideal hyperbolic prism 
constructed from t's angle data relative to x, as described in the previous section. Now 
simply let the energy be 

E{x) =Y,Vt{x). 

tGP 

The points where E attains a maximum will be of great interest to us. It is 
worth noting that E is a. continuous function on the compact set N^;, so attains its 
maximum. In fact it is an immediate consequence of lemma |5| in section 2.1.4 that the 
angle system with maximal energy is in fact in N. An understanding of this as well 
as the energies behavior at critical points comes from understanding its differential. To 
interpret it let a, h and c denote the edge lengths opposite to the angle A, B, and C in 
the hyperbolic triangle determined by the angle data {A, B, C} (which once again exists 



since k (x) < 0). In formula g of section 2.1.4 the following formula is produced: 

Formula 4 dE^ = eT ^i(^)'^^ ^^^^ 

„ , . 1 / / cosh(a) — 1\ , / cosh(6) — 1\ , / cosh(c) 
Ei{x) = - [In I V -l'^ V 



2 

where a*(x) = A, ai £ t, t's angle data is {A,B,C}, and the a,b, and c are determined 
as above. 
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In order to best exploit this formula it is necessary to note (see the proof of 
corollary |l] in section 2.1.4) that the tangent space at any point of Nx can be described 



explicitly as the translation to that point of the span over all edges of the vectors 

W^e = 71 + /?! - (72 + /32) 

as in figure |1.2| . The loveliness of this energy can now be expressed in terms the following 
observation about its critical points (i.e. where dE = 0). 

Observation 1 From the above formula and the above description o/Tp(Nx) we see at 
a critical point x of the energy in a conformal class satisfies 

= dE{w^) = In ( ^^^H^^)-^] _ In " ^ 



2 J \ 2 

for all edges e with the notation coming form as usual from figure [7]J. Hence at a critical 
point of E we have that oi = 02 and the set of —1 curvature triangles formed from the 
given angle data fit together to form an actual constant curvature surface. 

I will call the Delaunay angle system of such a critical point a uniform angle 
system. The question becomes: how many (if any) uniform structures can be associated 
to a given angle system? From above we know there is at least one internal maximum 



and in fact E is strictly concave down (see lemma |^ in 2.1.4 ). So any critical point is a 
maximum and unique, as needed. 

Theorem 3 (Discrete Uniformization Theorem) Ifx{M) < andx is a Delaunay 
angle system where {9'^{x)} satisfies (*) then x is conformally equivalent to a unique 
uniform angle system. 

It is worth noting that there is a rephrasing of the above uniformization con- 
struction in terms of a solution to a disk pattern problem, namely as a corollary of 



theorem U in section 2.3.1 



Theorem 4 // one is given a topological triangulation of a surface with x(^) < '^^^ 
a set {9^} associated to the edges satisfying both condition (★) and that at each vertex 
^{e&v} — 9^ = 2vr, then there is a uniquely associated constant curvature surface on 
which the given triangulation is realized as a geodesic triangulation and the circles in 
which the triangles live meet with the specified 9^ angles. 
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Solutions to problems similar to this are well known. The first solution to such 
a pattern problem goes back to Koebe and pattern problems closely related to the above 
corollary were implicit in the work of Andreev Q and then rediscovered and articulated 
in this language by Thurston (see p!7||).The above theorem (actually a generalization of 
it found in p.3.1 ) is a generalization to arbitrary angles of the convex ideal case of the 



Thurston-Andreev theorem. The first solutions to these pattern problems using energy 
methods (as far as I'm aware) can be found in Q, and energy methods using hyperbolic 
volumes (in the toroidal case) have there earliest versions in and later in p!^ . 

1.2.2 Continuous Uniformization 

Now we'd like to mimic the uniformization procedure in the previous section 
for metrics. For metrics by uniform structure I will mean a metric e'^'^g with constant 
curvature. It is useful to note that = e^'^(— + k). 

For starters let us note in the metric world we still have 

Lemma 3 (The Metric Teleportation Lemma) Every metric on a surface is con- 
formally equivalent to a metric with either negative, positive , or zero curvature. 

Proof: This is a consequence of the Fredholm alternative, which says if / € C°°(M) 
satisfies /^^ fdA = then / G A(C°^(M)). Since /^//c - 2T:x{M))dA = we have there 
is a smooth (j) satisfying A0 = k — 27rx(M) and hence k^ = e^'^(— A0 + k) = 27rx(M)e^'^ 
as required. 

q.e.d 

With this observation in our x{^) < world we will restrict our attention 
to metrics with strictly negative curvature. It is worth noting that as in the discrete 
teleportation lemma the problem is linear. 

Now given a varying negative curvature metric h if we chose a p = {pi, . . . 
and topological triangulation with p as vertices then we could measure how close to 
uniform h is with the energy of the previous section. Namely we could let 

E/.(p) = ^tHp), 
teKp 

where V^'^(p) is computed using the angle data associate to the "triangulation" viewed in 
the h metric. This of course means connecting the needed vertices of p with h geodesies 
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and measuring the resulting h angles. Notice this makes sense for sufficient h dense 
points since by the Gauss-Bonnet theorem for geodesic triangles the angles in a triangle 
{A,B,C} will indeed satisfy 

k\p) = A + B + C-'iT<0. 

This gives us a measurement of uniformity relative to {p} and Kp and to rid 
this dependency it is natural to average this measurement over all complexes as in section 
|1.1| . To do so weight the point distribution and decide how to form the complexes using 
a fixed metric g and density A, and then compute EA(Eh). It's worth noting that it 
may be necessary to shrink the decision radius down a bit to make sense out of this 
construction, since we'd like points in a triple to be within the injectivity radius of each 



other in both metrics. As with the random variables of section 1.1 this random variable 



is naturally expressed as the sum over all the faces t £ Kp of a function dependent only 



on the data associated to an individual face, so we would expect the computations of 1.1 
to go through. 

It is clear that the higher the vertex density the higher the percentage of tri- 
angulations involved, hence the more sensitive this energy should be to measuring the 
uniformity of the curvature. With this as inspiration we proceed as with as with our 
Euler characteristic computation and take the limit as the density of vertices goes to 
infinity. Let 

E{h) = lim EA(Eh-Eg), 

A— >oo 

and note the second term E^ is independent of h and needed only to normalize the 
computation. We will show that among h = e^'^g with with < Q this formula is 
equivalent to 

E{h) = -[ WV^W"^ + {A(j)-k) ln{A(j)-k) + kin \k\dA. (1.5) 

The energy is clearly scale invariant and k^ = e~^'^(— + k), so we see that 
it is natural to consider the energy as a function on the "energy norm" closure of 

y = {0 G C7°° I [ ^dA = and - A(j) + k < 0}. 
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Using the energy norm allows us to use compactness arguments to gaurentee the existence 
of a maximum, just as in the discrete case (for the details of everything that takes place 
in this section see section |4.2| ) . Also as in the discrete we can easily force the maximum 
into the the interior (where esssup{—A + k) < 0). This requires looking at the energy's 
Frechet differential, which is quite revealing. 

Formula 5 The Frechet derivative of E at (p in the direction ijj is 

DEitp) = - Ai;ln\k^2^\dA. 
Jm 

This formula for Frechet derivative guarantees that kh is constant and (p is 
smooth. The easiest way to understand this is to recall the Fredholm Alternative de- 
scribed in the poof of lemma ^ and observe... 

Observation 2 The above formula insures us that if we knew the critical (j) to he a C°° 
function then ln\kh\ would orthogonal to the image of C°°{M) under the Laplacian, 
so by the Fredholm alternative ln\kfi\ is constant. Hence we indeed see that a critical 
point of E should he a metric of constant negative curvature. 

Just as in the discrete case at this point concavity comes to the rescue to 
gaurentee uniqueness. Since everything is C°° we really can use the fact that the Frechet 
Hessian at </> G ^ applied to (V', V') 

D^E{^, V) = - / I VV'P + j^d^o (1.6) 

is strictly negative at a non-zero ip to gaurentee the critical point is in fact unique. 

So just as in the discrete case we arrive at the uniformization theorem (although 
only when x{M) < 0). 

Theorem 5 (The Metric Uniformization Theorem) Every metric is conformally 
equivalent to a unique metric of constant curvature. 

1.2.3 The Determinant of the Laplacian and Entropy 

The use of energies to solve the uniformization problem is not at all new, and 
it is nice to relate this process to some other more familiar energies. The method of 
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producing other energies is to average over metrics using the varying metric to distribute 
the points. So in this section I will superscript the expected values with the metric used 
to construct the point distribution. 

The first variant I will describe has the advantage of not needing to be restricted 
to a conformal class. It involves a function on the space of metrics I call the entropy 



Hih) = [ khln\hh\dAh. 



The fact is that 



E^{h) = lim (Eh(Eh) - Ef (Eg)) = H(h) - H(g). 

So viewing g as fixed we have an energy on all metrics - which can be checked to optimize 
at a constant curvature metric and be strictly concave down in a conformal class. 
More interestingly is the following energy. 

Theorem 6 There is a constant C such that on the set of all metrics h conformally 
equivalent to g of a fixed area 

E2{h) = Ei{h) - E{h) = lim (Eh(Eh) - Ef (Eh)) = ln(det(Ah)) + C. 

A^oo 



As we shall see in section 4.1.2, the key to being able to make this interpretation 
is the beautiful integral formulation of the ln(det(Ag)) due to Polyakov (see fl^). The 
ln(((iet(Ag)) was found to have the uniformization property by Osgood, Phillips, and 
Sarnak (see Q). In |Q| Yves Colin de Verdiere suggested that energies related to circle 
pattern problems might be related to the determinant of the Laplacian. This procedure 
provides such a relationship. 



1.3 Three-dimensional Dreams 

The entire two dimensional story has a three dimensional wishful analog. 



Discrete Uniformization 

The discrete case is essentially an idea do to Casson for deforming a triangulated 
three manifold into one of constant curvature (as I understand it the approach presented 
here is the Lagrangian dual of his approach). 
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The idea is to suppose you have a topological triangulation of a three manifold, 
let N be the possible dihedral angle data such that the sum of the angles about an edge 
is 27r and such that the data can be used to produce a hyperbolic simplex. Call the 
dihedral angles value at e in t q*. Now we can deform our angle taking the span of the 
transformations We^'^'^ which adds one unit of angle to a^^ as subtracts one unit form 
ag2, for neighboring simplexes ti and t2- Let N be the angles geometrically equivalent 
to X. 

On N we can now place the energy 

teP 

where the sum is over all the simplexes t and Vt{x) is the hyperbolic volume of t relative 
to x's angles data. 

Now, as in the the two-dimensional discrete case, the critical points of this 
energy are a sets of simplexes which fit together. To see this we recall Schafli's formula 
which tells us 

set 

where Z* is edge length of e. 

So at a critical point we have 

dE{w%''^) = ll^-ll\ 

Unlike in the two dimensional case we fail to have good boundary control (or 
Casson would have already uniformized the manifolds of interest) . In particular simplexes 
may collapse and enough collapsing may take place that even the topological type of the 
complex could change. This is where the randomizing may help. 

From the Discrete to the Continuous 

Now just as before we randomly Delaunay triangulate. We use a Poisson point 
process to distribute the points and then assign a simplex to a quadruple on a small 
enough sphere passing through the four points which is empty of other points. With 
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high probabihty this complex will form a triangulation. However not as canonically as 
in the two dimensional case, since there are no natural geometric simplexes and faces 
only natural edges and vertices. 

For any set of vertices forming a triangulation we can construct a point in some 
N, just as in the discrete case. It is essential to note that while we have no faces and 
hence no dihedral angles, we have dihedral angles at both end points (using the geodesic 
directions) so can average them and form the angle data of some y in some N. So we 
may use the energy to once again form a random variable which can be computed. 



Continuous Uniformization 

Just as in the two dimensional case we may now use 



E{h) = lim EA(Eh - Eg), 

A— KXD 



to form an energy on the space of metrics. Currently not much is known about this 
energy, but the hope of course is that it will form an energy with which the techniques 



of section 4.2 can be carried out. 

In the end hopefully one will be able to see that a three manifold admitting 
a metric of variable negative sectional curvature will accept one of constant negative 
sectional curvature. 



Chapter 2 

The Discrete Uniformization 
Theorem 



This chapter is dedicated to proving the generahzation of the Thurston- Andreev 
theorem mentioned in the introduction. It is presented in a potentially strange order from 



the point of view of presentation in section \L2. In section 2.1 the details of the energy 



argument presented in 1.2.1 are given ending in a theorem giving conditions under which 
triangle angles in a topological triangular decomposition can be conformally deformed 
to the angles of a geodesic triangular decomposition of a hyperbolic surface. Also in this 



section a connection between the disk patterns mentioned in section 1.2.1 and certain 
hyperbolic polyhedra is explored and exploited. 



In section 2.2 we explore the linear part of the discussion in section 1.2.1, and 
prove a warm up case of the general ideal convex Thurston-Andreev theorem. This 
warm up case is theorem |l] form the introduction including the possibility of M having 
boundary and has as an immediate consequence the teleportation lemma (lemma |2|) . In 
section the general convex ideal case of the Thurston-Andreev theorem is presented 
and dealt with. 



2.1 The Triangular Decomposition Theorem 
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Figure 2.1: A Covector and a Vector 



2.1.1 Statement and Notation 

Throughout this paper M will denote a compact two-dimensional surface with 
x(M) < 0. By geometry I will mean a hyperbohc structure. Uniqueness of geometries, 
triangulations and disk patterns is of course up to isometry. 

The main theorem in this section really should be stated for the following 
structure which generalizes the notion of triangulation. 

Definition 1 Let a triangular decomposition, T, be a cell decomposition of M that lifts 
to a triangulation in M 's universal cover. 

We will keep track of the combinatorics of such a decomposition by denoting 
the vertices as {vi, ...vv}, the edges as {ei . . . eg} and the triangles as {ti . . . tp}- For 
convenience I will let E, V, and F denote both the set of edges, vertices, and faces and 
the cardinalities of these sets. Similarly for the subsets of E and V contained in M's 
possible boundary, denoted dE and dV. Let {e G S} denote the set of edges on the 
surface in a collection of triangles S, and let {e G v} denote all the edges associated to 
a vertex t; as if counted in the universal cover. The set of triangles containing a vertex 
{t G v} has the special name of the flower at v. 

Note that in a triangular decomposition there are 3F slots {ai} in which one 
can insert possible triangle angles, which we will place an order on and identify with 
a basis of a 3F dimensional real vector space. With this basis choice we will denote 
this vector space as R^^, and denote vectors in it as x = '^A^ai. Further more let 
a' be a dual vector such that a*(aj) = Sj. With this we will view the angle at the 
slot ai as a^{x) = A^. It is rarely necessary to use this notation and instead to use 
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Figure 2.3: The p" Covector 



the actual geometry. We will denote a vector by placing the coefficients in a copy 
of the triangulation with dashed lines and covector will contain its coefficients Ai in a 
copy of the triangular decomposition with solid lines. Thick lines will always denote a 
boundary edge, as in lower left corner of vector and covector in figure ^T^. If in the 
picture we mean the non-specified values to be arbitrary we will surround the picture 
with a loop (see the vector in figure |2.1|) and if we mean the non-specified values to be 



zero the picture will not be surrounded (see the covector in figure The pairing of a 
vector and a covector denoted ^"'cKj) can be viewed geometrically by 

placing the copy of the triangular decomposition corresponding to the vector on top of 
the triangular decomposition corresponding to the covector and multiplying the numbers 
living in the same angle slots to arrive at Yl{ai} (see figure [2.21 ). For a triangle t 
containing the angle slots a^, aj, and let d^{x) = {A^ , A^ ,A^^ and call S{x) the angle 
data associated to t. 

In order to live on an actual nonsingular geometric surface all such angles should 
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Figure 2.4: The /* Covector 

be required to live in the subset of R^^ where the angles at an interior vertex sum to 
27r and the angles at a boundary vertex sum to vr. Let p" be the covector living in the 
flower at v defined as in figure p.3| , this encourages us to choose our possible angles in 
the affine fiat 

V = {xe R^^ I p''(x) = 27r for all v G V - aV and p''(x) = vr for all v G dY}. 

To further limit down the possible angle values we define the covector /* as in 



figure 2.4 and note by the Gauss-Bonnet theorem that 

k\x) = l\x)-Tr 

would be the curvature in a geodesic triangle with angle data d^{x). We will now isolate 
the open convex subset of V where the curvature is negative and all angles are realistic. 



Definition 2 Let an angle system be a point in 

N = {x G V I k*(x) < for all t and a'(x) G (0,7r) for all a'}. 

Note the actual angle data of a geodesic triangulation of a surface with negative curvature 
has its angle data living in this set. 

Observe that from any set of triangle data d*{x) = {A,B,C} with the angles 
in (0, vr) and A + B + C — tt < Owe may associate an actual hyperbolic triangle, call this 
triangle t{x). Suppose the triangles in {t{x) \ t G T} fit together in the sense that all 
the corresponding edges are the same lengths. Then T being a triangular decomposition 
implies every open flower is embedded in M's universal cover and when the edge lengths 
all agree this flower can be given a hyperbolic structure which is consistent on flower 
over laps. So we have formed a hyperbolic structure on M. 
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Figure 2.5: The Vector 

Definition 3 Call an angle system u uniform if all the hyperbolic realizations of the 
triangles in u fit together to form a hyperbolic structure on M . Let T(u) denote the the 
iesic triangular decomposition corresponding to T and u. Further more let 



U = {u G N I u is uniform }. 

In section 2.4 we will attempt to take a point in N and deform it into a point 
of U. Such deformations are located in an affine space and I will call them conformal 



deformations (see the introduction to section 1.2 to motivate this terminology). To 



describe this affine space for each edge e £ E — dE construct a vector We as in figure |2.5| . 

Definition 4 A conformal deformation will he a vector in 

C = span{we I for all e £ E — dE}, 
and call x and y conformally equivalent if x — y € C. 

The first thing worth noting is that x £ V and y is conformally equivalent to x then 



as an immediate consequence of geometrically pairing the covector in figure 2.3 with the 



vector in figure 2.5 we have for v £ V — dV that 



p^{y) =p-lx + Y, B'wA = p^{x) + = 27r + 0, 

\ eeP / eGP 

SO y is also in V. Similarly for v G dV. 

To combinatorially understand the points in N which we may conformally de- 
form into uniform structures it is useful to express a particularly nasty set in the bound- 
ary of N. 
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Figure 2.6: The ipf Covector 



Definition 5 Let t be called a legal with respect to x £ 9N if d!'{x) = {A^,A^,A'^} / 
{0, 0, vr} yet either = or for some i we have = 0. Let 

B = {x £ 9N I X contains no legal t}. 

We will prove the following theorem. 

Theorem 7 If there a uniform angle system conformally equivalent to x then it is 
unique, and for any angle system x with {x + C){^B empty there exists a conformally 
equivalent uniform angle system. 

Much of what takes place here relies on certain basic invariants of conformal 
deformations. To describe them for each triangle t and e £ t we form the covector ipf 
as in figure |2.6| . For each edge e G dV we will denote ifjf as ip*^ while for each edge 
e £ V — dV associated with triangles ti and t2 we will let 

We will call the ip^ covector the formal angle defect at e. Let the formal intersection 
angle be determined by 

when e G E — dE. and 

when e G dE. 



TT 

2 



ijj^x) 
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Looking at the pairing between a covector ijj^ and the vectors spanning C in 
figure we see if y is conformally equivalent to x then 

V /eP / /eP 
and indeed for each relevant edge e we see ij:^ and are conformal invariants. 

It is worth noting the trivial but extremely useful fact that the curvature as- 
sumption gives us some control of the angle discrepancy. 

Fact 1 When x S N we have i/^fii^) ^ f ); when x G 9N we have iljf.{x) S 

r —TT 7r1 
L 2 ' 2j • 

Proof: 

Let d*'(x) = {yl,5,C} and note since B + C<A + B + C = l\x) < vr and 
^ < TT we have 

TT ^^,e 5 + + C vr 

— < — < = < < -• 

2 2 ~ 2 ~ 2 2 

The boundary statement follows from the possibility of these inequalities becoming equal- 
ities, q.e.d 

As we shall see, the fundamental reason why the above triangulation theorem 
is related to circle patterns and polyhedra construction is that at a uniform angle system 
6'^{u) can be geometrically interpreted as the intersection angle of the circumscribing 
circles of two hyperbolic triangles meeting along e £ E — dE. 

Equivalently 9*^ will be realized as the dihedral angle at an edge of a polyhedra 
in the class / corresponding to e. It's high time to describe this class /. 

2.1.2 Polyhedra in the Class I 

In this section we will construct and examine a bit of the geometry of the infinite 
sided ideal polyhedra that will arise in the proof of theorem |^. These polyhedra can be 
constructed out of building blocks each in the form of an ideal prism. 

Definition 6 Place a given t{x) on a copy of C . Let Pt{x) he the convex hull 
of the set consisting of t{x) unioned with the geodesies perpendicular to this C 
going through t{x) 's vertices. See figure \2.% 
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Figure 2.7: The Ideal Prism Pt{x) 

Given a triangular decomposition and uniform angle system u let T{u) denote 
the hft of T(u) to H^. With this notion we may now construct the polyhedra in I. 



Definition 7 Let the class of polygons I be those constructed by placing T{u) on an 



C H"^ and forming UteT ^t{u). 



Notice abstractly a polyhedron P in I is an ideal polyhedra symmetric upon the 
reflection through some plane upon which lives a group of isometries forming a compact 
surface which when extended to are isometries of P. 

Our first observation will be that if we happen to know a T and u forming P, 
then from this data we can easily construct the dihedral angles. We will also refer to 
the dihedral angle associate to an edge of T, despite the fact there is always in fact a 
symmetric pair of such angles. 

Formula 6 If P is constructed from T and u as in definition then the dihedral angle 
at an edge of P associate to the edge eofTis given by 9^{u). 

Proof: Note that the needed dihedral angle is the sum of the angles in Pt-^ (u) and Pt2 
corresponding to e. It is these angles that will be computed. 

To do this simply note that Pj . (n) can be decomposed as into three ideal tetra- 
hedra as in figure |2.8| , where (i*'(n) = {A,B,C} with the angle slot containing the A 



coefficient across form e. The labeled angles in the figure 2.8 are the internal angles 
closest to the label. 

Now use the fact that at the vertex of an ideal tetrahedron the angles sum to vr, 
to form ten linear equations in the labeled unknowns. Solving in terms of A = A^,B = B^, 
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Figure 2.8: A Decomposition of Pt{x) 

and C = one finds 

, TT + A-B-C 
x = l = ^ 

as needed. 

q.e.d 

Note 1 The polygons in the surface with boundary case now are extremely non-convex 
and have faces between the two hemispheres. Note that the {6^{x)} still represent the 
dihedral angles in these faces. 

Let G be a compact surface group extend form the through which P is 
symmetric and under which P is invariant. To be in the class / such a group exists. 
When such a G is chosen we wih be interested in the volume of P/G. Notice given such 
a G we may choose a geodesic triangulation which descends to a triangular decomposition 
T of /G with an associated uniform structure u from which P is constructed as in 
definition ^. Using this triangulation we have. 

Formula 7 The volume of P/G is ^t^T^Vt{u) where 
Vt{u) = K{A) + K{B) + A(G) + A 



vr 



A-B-G 
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^,Ti + A-B-C\ ^['K + B-A-C\ .f-K + C-A-B' 
+A + A + A 



2 

with A is the Lobacevskii function 



a 



A(a) = - / ln(2|sin(t)|)dt. 
Jo 

Proof: First note the volume of P/G can be reduced to the volumes of individual Pt by 
noting the volume is X^^gx ^t{u), where Vt{u) is the volume of Pt{u). 

Its useful to get a formula for Vt{u) Recall if a tetrahedra has angles a, (3, and 
7 meeting at an ideal vertex then its volume is 

A(a)+A(/3)+A(7) 

where A is the Lobacevskii function. 

In the proof of formula ^ we decomposed Pt{u) with d*(n) = {A,B,C^ into 
three ideal tetrahedra, and found linear equations determining all the angles in these ideal 
tetrahedra in terms of A, B, and C. We already wrote down the angle corresponding to 
X and / and can further note 

TT + A-B-C 

y = t; 



TT + C -A-B 
z = 



TT-A-B-C 
y -n = 



m = C. 

Plugging these angles in the tetrahedra's volume formula gives us the needed formula. 

q.e.d 

Notice that when P is constructed form T and an edge e of T is associate 
a dihedral angle tt that P's edge corresponding to e is fake in the sense that the two 
triangular faces meeting at this edge share the same polyhedron face. If a polyhedra's 
face is triangulated then a change in this triangulation will not effect the polyhedra so it 
will turn out convenient to have articulated the topological cell divisions which naturally 
arise form polyhedra in /. 
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Definition 8 Let a polygonation of a surface be a locally finite cell division such that the 
closure of a cell is an embedded polygon and the intersection of two polygons is empty, 
contains a single point, or contains one edge and the two vertices associated to the edge. 
Let a polygonal decomposition P be a cell decomposition of a compact surface which lifts 
to a polygonation in its universal cover. 

To keep track of all the the combinatorics of such a cell division we will use the 
same notation as we did for triangular decompositions. 

At this point it is useful to name what turns out to be the appropriate home 
of the possible dihedral angle assignments. Let {ej} be the set of edges in a polygonal 
decomposition and just as we did with the angle slots let them correspond to the the 
basis vectors of an E dimensional vector space, which we will denote with this basis 
choice. We will be viewing this as the space of possible angle discrepancies. Denote these 
vectors as p = ip'^^ ei . 

Notice the were covectors in the previous section. We can motivate this 
abuse of notation in the case that P is triangular decomposition by letting 

* : R3P ^ R^ 

be the linear mapping given by 

^{x) = Y,r'{x)ei, (2.1) 

and noting by the above lemma that we do indeed hit the dihedral angle discrepancies 
when using a uniform angle system. As a further justifiable abuse notation we let 9^^ (p) = 
TT — e'(p) = TT — -0^' when e E — dE and 6^^{jp) = f — e*(p) = f — V'^' when e* e dE. 

At this point it is convenient associate explicitly the data contained in polyhedra 
in the class / with certain disk patterns. 

2.1.3 Ideal Disk Patterns 

In this section we will discuss the relation of these polyhedra to disk patterns 
specified by combinatorial and topological data. The topological data comes in the form 
of a polygonal decomposition P and the extra data associated to P will be a point 
p G R^. 
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Figure 2.9: An Ideal Disk Pattern 



To begin to articulate the pattern here It's necessary to define a particularly 
nice polygonal decomposition. 

Definition 9 A polygonal decomposition on a geometric surface is called circumscribable 
if it is a geodesic polygonal decomposition where each polygon is circumscribed by circle 
in M 's universal cover. 

Given a circumscribable P we may use a point p G to keep track of all 
angles of intersection between the circumscribing circles of the polygons meeting at e 
(namely at e this angle is 9^{p)). Simply to articulate the sense of angle to be used 
here it is convenient to introduce P* which is a triangular decomposition associated to 
P by triangulating each polygon. It is occasionally useful to have an explicit grip on 
this triangulation so we may assume that we triangulate each polygon with a fan as in 
figure 2.10| . Note this triangulation can be chosen to be a geodesic triangulation if P is 



circumscribable, since all the polygons are then convex. Also if we are given a p E 
relative to P we will let p be the angle discrepancy assignment on P* which is ip^{p) 
on the the edges of P and on the new edges, and hence corresponds to geodesically 
triangulating a circumscribable decomposition. 

Note given a uniform structure u associated to P* that every edge sits between 
two triangles ti and t2, and the sense of the angle used here can be chosen relative to 
these geodesic triangles. Let the intersection angle 9^{p) be a number in (0, vr) if the 
vertex of ti not on e is out side t2's circumscribing disk, in (vr, 2-it) if this vertex is in t2^s 
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Figure 2.10: A Fan Associated to a Polygon 

circumscribing disk, and exactly vr when this vertex is on t2^s circumscribing circle. 

Definition 10 Let an ideal disk pattern be a collection distinct disks on a geometric 
surface whose boundary circles are in one to one correspondence with the circumscribing 
circles of a circumscribable polygonal decomposition. Notice in such a situation using 
a P* we naturally have the triangles needed to associate p € (0, 2tt)^ to keep track 
intersection angles at each of the polygonal decomposition's edges. 

The pattern problem is given a topological polygonal decomposition and a 
suitable p G to assert the existence and uniqueness of an ideal disk pattern. Such 
assertion are equivalent to such assertions about polyhedra in the class /. 

Observation 3 Every ideal disk pattern can be associated a unique polyhedra in I. The 
polyhedra associated to a pattern is convex if and only ifp E (0, vr)^, in which case we will 
also call the pattern convex. To each polyhedra in I and choice ofG as described preceding 
formula there is a uniquely associated ideal disk pattern. Furthermore under these 
correspondences the dihedral angle of the polyhedra associated to an edge e is precisely 
the intersection angle between the circumscribing circles of the polygons sharing e in the 
ideal disk pattern. 

Proof: Note that given a circumscribable polygonal decomposition that we may form 
a geodesic triangular decomposition P*(u) which can be associated an element of / 
via the construction in definition |^ Similarly given an element P in / and a G as in 
the discussion preceding formula |^ we know there is a geodesic triangulation invariant 
under G which descends to a geodesic triangular decomposition of H'^/G. Let P be the 
circumscribable decomposition formed by ignoring the edges where 'ip^{u) = 0. 
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Figure 2.11: A Conformal Mapping 

To see the angle correspondence we will recall a map from a specified C H^^ 
to the upper half of the sphere at infinity S*^ (using the usual conformal structure of 
the sphere at infinity) which sends circles in to circles in S^. To form this map 
send a point in this fixed to where the geodesic perpendicular to this hits 5*2° (as 
in figure |2.11| ). Note any circle can be sent to what we view as the center of via a 
hyperbolic isometry preserving , were by symmetry it is sent under this mapping to a 
circle at infinity. This isometry bringing the circle to what we view as the center induces 
a Mobius transformation on the sphere at infinity and preserves the set of geodesies used 
to form this map's image - so the fact that Mobius transformations send circle to circles 
on now gives us that indeed the image of any circle is a circle. 

So a neighboring pair triangles ti and t2 have circumscribing circles in sent 
under this map sent to circles at infinity intersecting at the same angle and going through 
the ideal points of the neighboring Ptj^^u) and Pt^iu). But these circles at infinity are 
also the intersection of S°° with the spheres representing the hyperbolic planes forming 
the top faces of Pti{u) and Pt^iu). So the intersection angle of these spheres is precisely 
the dihedral angle, which is now seen to be the intersection angle of the circles on the 
sphere at infinity, or finally the original intersection angle of the circumscribing circles. 

The convexity assertion can be seen immediately by looking at the two neigh- 
boring Pt{u) prism's. 

q.e.d 

At this this point there are two natural questions about such patterns and 
polyhedra, namely when they exist are they unique and are there nice way to insure 
existence? The first question will be answered in the next section with the following fact: 
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a pattern/polyhedra is determined uniquely by its/an associated topological polygonal 
decomposition P and p G (0,27r)^. The existence issue will be handled in sections p.2| 
and p.3.2| where necessary and sufficient conditions for p G (0, 27r)'^ relative to P to be 
associated to a convex pattern/polyhedron with this given data will be presented. 

With respect to these existence and uniqueness results I will only deal with the 
terminology of ideal disk patterns from here on out. 

2.1.4 Proof of Theorem ^ 

Now we will prove theorem |^. The proof relies on the energy introduced in 



section 1.2.1 which lives on N, and which at a uniform u agrees with the volume of P/G 
with P constructed from P(u) as in definition ^ and G the surface group associated to 
P(u). We may assume we that either P is a triangulation or form it we have constructed 
the triangulation P*, so we will denote in T. 

Let the energy be the of the volume of the abstract disjoint union IJ^grj, Pt{x) 

given by 

E{x) = Y,Vt{x). 
teT 

Since hyperbolic objects optimize at fat objects we will be maximizing this energy, and I 
suppose to call this an energy (in the physical sense) I really should negate it. However 
I like both the fact that hyperbolic objects like to be fat and the term energy, so I will 
simply warn the reader about this odd terminology. Note we get an explicit description of 
the energy from formula ^. Using this formula we can differentiate to find E^s differential, 
dE^. As usual for a function in a linear space like R^^ we use translation to identify 
the tangent and cotangent spaces at every point with R^^ and (R^^)* and express our 
differentials in the chosen basis. From the formula for the Lobacheski function we have 
dE^ = ^ai£T ^ii^)'^^ ■^i*^ 

sin^(^) sin( "+-^,^-^ ) \ 



Ei{x) 




sin ( T^+B-A-C \ gjj^ ( TT+C-A-B \ ^-^ ( -K-A-B-C y 



2 ; — V 2 -)sin(^ 2 ;/ 
where a*(x) = A, ai G t, and d^(x) = {A, B,G}. After a little bit of trigonometry this 
can be simplified to 

^, , 1, / sm{Af{cos{B + G)+cos{A)) 



2 V (cos(^ + C) + cos(B))(cos(A + B)+ cos(C)) 
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To compute further let o, b and c denote the edge lengths opposite to the 
angles A, B, and C respectively in the hyperbolic triangle determined by {A,B,C}. 
From elementary hyperbolic geometry we know that 

, , cos(C) cos(-B) + cos(^) 
cosh a) = ■ fTj\ ■ (n\ ' 

allowing us to simplify cosh(a) — 1 to 

cos(S + C)+ cos{A) 
sin(S) sin(C7) 

So plugging this in the above formula we arrive at... 

Formula 8 dE^ = X]q,,£x -^«(^)'^* ^^^^ 

^ , , ^ f , f cosh(a) — 1 \ - / cosh(5) — 1 \ , / cosh(c) — 1 
Ei{x) = -- In - In 7, - In ' 



2VV2 J \ 2 J \ 2 

where a*(x) = A, ai £ t, = {A,B^C}, and the a,b, and c determined as above. 

Now let's look at the conformal class of a point x G N, and define 
Definition 11 Let 

Nx = (x + C)f|N 

and call it the conformal class of x. 

Recall that Tp(Nx) = spanjwe | e G E — 9E} from definition |^. Let e be the edge 
between ti and t2 with d*' = {Ai,Bi, Cj}, Ai corresponding to the angle slot across from 
e in tj, and the Oi, bi , and Cj determined as above. Then we may observe... 

Observation 4 From the above formula and the description o/Tp(Nx) a critical point 
y of the energy when restricted to x 's conformal class satisfies 

= dEy{v,) = In ( _ In ^^"^^("^^ " ' 



for all edges e. Hence at such a critical point we have that ai = 02 and y is a uniform 
angle system. 
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With this observation the the existence and uniqueness of such critical points 
is equivalent to the existence and uniqueness of uniform structures conformal to a fixed 



Lemma 4 //Nx contains a uniform angle system this angle system is unique in Nx- 

Proof: This will follow form the fact E is concave down. What we are really asking 
for is that the energy's Hessian is negative in the C directions at all points in Nx- In 
fact we will show something consider ablely stronger, namely that the Hessian is negative 
throughout all of the open subset of R^^ satisfying N's open condition and in all of 
R^^'s directions. (A fact which has several interesting application to the production of 
triangulations with special symmetries.) Note that in this setting E and hence its Hessian 
splits up into a sum of independent functions associated to each triangle - so we only need 
to show the 3x3 matrix corresponding to a fixed triangle is negative definite. Assume 
we are in the triangle with angle slots a, (3, and 7 and = {A, B, C}. Recalling that 
a linear change of coordinate will not effect whether the Hessian is negative definite or 
not, I found it useful to use the coordinates {A + B,A + C,B + C) for this computation. 
To write down the Hessian in these coordinates it is use full to introduce the functions 



one. In particular we may prove the uniqueness assertion in theorem ^ 




and 




With these named a direct computation shows 



Hess = 




-1 



/ G{B,C,A) 
FiC,B,A) 
\ FiB,C,A) 



F{C,B,A) F{B,C,A) \ 
G{A,C,B) F{A,B,C) 
F{A,B,C) G{A,B,C) j 



To see this matrix is negative definite we can use the following easily derived 
condition: a symmetric 3 by 3 matrix a*- is negative definite if < 0, a^al — (03)^ > 0, 
and 



091 191 '-i/l\9 



alialf - alialf > 0. 
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Note through out this computation that the assumption that A;*(x) < or rather A + 
B + C < TT imphes aU the sin functions are evaluated at positive angle less than ^ and 
all the cos functions are evaluated at sums of angle of absolute value less than -1; so all 
such evaluations are positive. In particular the negative sign outside guarantees that all 
terms in the matrix, including the needed ag term, are negative. 

Now we need to satisfy the remaining two conditions; namely we'd like 

GiA, B, C)G{A, C, B) - {F{A, B, C)f = 

cos ( ^"'"f 

sin (4) sin (f) cos (4) cos (f ) 

to be positive, which by the observations made above it is. 
As well as needing 

G{A, B, C)G{A, C, B)G{B, C, A) + 2F{C, B, A)F{B, C, A)F{A, B, C) 
-G{A, B, C){F{C, B, A)f - G{B, C, A){F{A, B, C)f - GiA, C, B){F{B, C, A)f = 

32 (cos (ii+|+g))'cos ( --^+.f+^ ) cos cos 
sin (4) sin (I) sin (§) 

to be positive which once again as observed above is true. 

So indeed the volumes Hessian is negative definite throughout the open subset 
of R^^ satisfying N's open conditions. 

q.e.d 

Prom this we will be able to show the uniqueness statement claimed in the 
previous section. 

Corollary 1 An ideal disk pattern is uniquely determined by its topological ploygonal 
decomposition and associated angle discrepancy system p. An ideal disk pattern can be 
constructed form a polygonal decomposition P a,nd associated angle discrepancy system 
p precisely when there is a critical point of E in f]^- 
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Figure 2.12: The nie Vector 

Proof: Notice the choice of triangulation in forming P* from P at this step is irrelevant 
for both uniqueness and existence, since in an actual pattern any edges with a 9^{p) = vr 
will not be included in the polygonal decomposition's description. 

Using such a triangular decomposition note ^ has rank E since the pairing of 



-0*^* with the vector m-g^ in figure 2.12 satisfies "^{me^) = ej for each j. Further note from 



section 2.1 that the null space contains the E — dE dimension space C and is 

3F - E = 2E - dE - E = E - dE 

dimensional - so C is precisely the null space. In particular all angle systems which could 
conceivably hit a specified set of discrepancy angles {ip'^ip)} is in ^~^{p), which is x + C 
for some x. So the above lemma gaurentees the uniqueness and existence of the geodesic 
triangulation necessary to construct the ideal disk pattern under the above conditions. 

q.e.d 

Now its time to explore the existence of critical points. Given a pre-compact 
open set O and a continuous function on O we automatically achieve a maximum. 
For this maximum to be a critical point it is enough to know that E is differentiable in 
O and that the point of maximal E is in the open set O. 

One way to achieve this is to show that for any boundary point yo that there 
is a direction v, an e > and a c > such that l{s) = yo + sv satisfies 

^(/(0,e)) C N 

and 

lim^E{l{s)) > c, 
s^o ds 
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for all s G (0, e). This works because under these hypothesize E{l{s)) is continuous and 
increasing on [0, e) and i/q certainly could not have been a point where E achieved its 
maximum. 

It is useful to note that the compactness of O gaurentees us that l{s) eventually 
hits the boundary again at yi for some s > 0. So we may change the speed of our line 
and assume 

l{s) = (1 - s)yo + syi 

is the line connecting the two boundary points. So the remainder of theorem |7| follows 
by applying this above criteria to Nx and E and noting... 

Lemma 5 For every pair of points yo and yi in but not in B with l{s) P|N 7^ we 
have 

]imJ-E{l{s)) = 00. 
s^o as 

Proof: By using k''{x) and the previous lemma's notation for the angles in a triangle t, 
we can write the equation for the dual of a's coefficient Ea{x) in formula ^ as 

cos(^ — /c*(x)) — cos(^) 



In 



+ 21n(sin(yl)). 



(cos(S - k^{x)) - cos{B)){cos{C - k^{x)) - cos(C7)) 

Fixing a triangle t let d^iyi) = {Ai, Bi,Ci}; and note with this notation that the contri- 
bution to ^^j^g^^^ coming from the triangle t is given as 

ds 

= {A, - Ao)E^{l{s)) + (Bi - Bo)Ep{l{s)) + (Ci - Co)^^(/(s)). 

Since yo is on the boundary of N and not in B there is some triangle t such 
that {Aq, Bq,Co} 7^ {0,0, vr} however either k^{yo) = or some angle is zero. 

We will show that for any triangle t in this case lim^— >o -D*(s) = 00. Its useful 
to divide the possibilities into the following three cases. 

1. Where d*(yo) contains zeros but k^{yo) 7^ 0. 

2. Where fc*(yo) = and no angle is zero. 
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3. Where k^{yo) = and one angle in d^{yo) is zero. 

In the first case note that the only pieces of Dt which become infinite are of 
the form {Aq — ^i)21n | sin(^)| and further note that if Aq = and l{s) P) N 7^ then 
Ai — Aq > 0. So we indeed have lims_»o-Dt(s) = 00 as required. 

To understand the second case note that if k^{x) tends to zero then then we 
may rewrite La as 

La = -31n(sin(y4)) +ln(sin(S)) + ln(sin(C)) +ln|A;*(/(s))|. 

From this we find the part of Dt{s) that is not bounded is in the form 

-{Ao + So + Co - (^1 + B, + Ci)) In \k\l{s)\. 

Further note when k\yo) = and l{s) f| N / that k\yi) = Ai + Bi + Ci- {Aq + Bo + 
Co) < 0. So once again we have lim^-^o -D*(s) = 00 as required. 

The final case is a combination of the above two were we find the part of Dt 
that is not bounded is in the form 

D\s) = -{Ao + Bo + Co- {Ai + Bi + Ci))(ln \k\l{s))\ + ln(sin(A)) + 4(Ao - Ai) ln(sin 

The arguments above immediately imply the correct derivative behavior. 

Now note that lim^-^o B*{s) is clearly bounded on triangles with angles not not 
satisfying any boundary conditions, hence if it were bounded or 00 on triangles where 
cZ*(a;) = {tt, 0, 0} we would be done. The whole reason the B is bad set is that it is in 
fact finite. Using the same argument as above and Taylor expanding you find 

D\s) = {4k\yi) - 4k\yo)) In \s\ + g{s) = g{s), 

with g(s) bounded. 

So the proof is complete. q.e.d 

2.2 A Warm Up Thurston- Andreev Theorem 

In this section we will prove a warm up Thurston Andreev theorem. This section 
and the next have been made independent of each other so there is a bit of repetition. 
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Corollary ^ tells us that an ideal disk pattern is always unique when it exists, 
and we are now left to deal with the dilemma of finding good existence criteria. Here I 
will describe in detail the strict convex case where p G (0, Tr)^^-'^^) x (0, f relative 
to T a triangular decomposition. Linear conditions on the possible p will be produced 
which are necessary and sufficient for the p relative to T to be the data of an ideal disk 
pattern. 

Let 

D = |xGN| ^(x)g (0,7r)(E-^E) X (0,|)^^|. 

and call this the set of negative curvature Delaunay angle systems. These angle systems 
are remarkably easy to work with such angle systems and in fact... 

Observation 5 Every point of D has unique ideal disk pattern associated to it. 

Proof: This observation relies on the following fact which will be of interest in its own 
right. 

Fact 2 If X £ T> is conformally equivalent to a point in (9N where for some triangle 
d\x) = {0,0, vr}. 

To see this fact assume to the contrary that for some t and c we have d* (x + c) = 
{0, 0, vr}. Let e be the edge of t across form t's n and let ti be t's neighbor next to e if it 
exists. Note by fact § that that the conformally invariant ^/^^(x) G (0,7r) would (even in 
the best possible case when e is not on the boundary) have to satisfy the contradictory 
inequality ^j^x + C) = + tp^^ < 0. 

From this fact we have that if x G D then no element in x + C could possibly 
be in B and the observation follows form theorem |^ 

q.e.d 

Now lets explore certain two necessary conditions on a p = \E'(x) with x G D. 
The first condition is the condition related to the fact that the angles at the internal 
vertex in a geometric triangulation sum to 27r and at a boundary vertex sum to vr. 



42 



This condition is equivalent to the following simple lemma. 
Lemma 6 

= {p G I p satisfies (ni)}. 
Proof: First note that if p = ^(x) then 

So by choosing x £ V we see is included in 

W = {p£R^ \p satisfies (ni)}. 

Recall form the proof of corollary |^ that ^'(R^^) = R^. So we may express any p £ W 
as p = ^{x) and the above computation gaurentees x £ V as needed. 

q.e.d 

The second necessary condition is a global one; namely an insistence that for 
every set 5 of |5| triangles in T that 

M >7r|5|. 

eeS 

Verifying (722) relies on the following formula. 
Formula 9 Given a set of triangles S 

{eeS} t€S ^ ^ e&dS-dE 

with the t in ipfix) term being the triangle on the non-S side of e. 

Proof: 

{ee5} eeS~dE e<^SC\dE 



e(^S-dE eeSr\dE 
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teS ^ ^ eedS~dE 

with the t in ^f{x) term being the triangle on the non-S" side of e. Substituting the 
definition of A;*(x) gives the needed formula. 

q.e.d 

Note for any point x S N that —k^{x) > and form fact ^ that ^ — ^f{x) > 0. 
So removing these terms form the above formula strictly reduces its size and when 
summed up we arrive at (712). 

With these two necessary condition we have our first pattern existence theorem: 

Theorem 8 // 

p£ D= !^q£ (0,7r)(^-^^) X (o, | q satisfies (rii) and (n2)| 

then p is realized by a unique ideal disk pattern. 

By observation ^ above this would follow if we knew the following proposition. 
Proposition 1 

*(D) = D. 

It is this bit of linear algebra we now will tackle. Notice the fact (ni) and (722) 
are necessary gaurentees that ^'(D) C D, and we left to explore ^''s surjectivity. 

2.2.1 The Surjectivity of the Delaunay Case 

To see the surjectivity of ^ let's assume the contrary that that ^'(D) is strictly 
contained in D and produce a contradiction. With this assumption we have a point p 
on the boundary of ^'(D) inside D. Note p = ^{y) for some y G 9D. Furthermore note 
(C + y) n D is empty, since other wise for some w G C we would have {y + w) G D 
which along with the fact that ^ is an open mapping when restricted to V would force 
p = ^{y) = ^{y + w) to be in the interior of D. 

At this point we need to choose a particularly nice conformal version of y, which 
requires the notion of a stable boundary point of D. Before defining stability note since 
D is a convex set with hyperplane boundary if x G 9D such that {x + C) f^T) = (p, then 
(x + C) Pi 5D is its self a convex k dimensional set. 
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Definition 12 A point in x ^ 9D is stable if {x + C) {^T) = (j) and x is in the interior 
of (x + C) Pi 9D as a k dimensional set. Any inequality forming D violated in order to 
make x a boundary point will be called a violation. 

The key property of a stable point is that a conformal change w £ C has 
x + ew £ t)'^ for ah e > or for some sufficiently small e > we have x + ew must still be 
on 9D and experience exactly the same violations as x. The impossibility of any other 
phenomena when conformally changing a stable point is at the heart of the arguments 
in lemma and lemma ^ below. At this point subjectivity would follow if for a stable 
X G 5D we knew that ^'(x) could not be in D, contradicting the choice oi p = ^'(x) as 
needed. 

We will prove this by splitting up the possibilities into the two cases in lemma 
1^ and lemma ^. 

Lemma 7 If x £ 9D is stable and = for a* in a triangle where k*{x) < 0, then 

^{x) is not in D. 

Proof: Look at an angle slot which is zero in triangle to satisfying k^°{x) < 0. View 
this angle as living between the edges cq and ei. Note that in order for x to be stable 
that either ei is a boundary edge or the ewg^ transformation (with its positive side in 
to) must be protected by a zero on the — e side forcing the condition that ewei ^ > 
or else for small enough e we would have x + ewe^ would be a conformally equivalent 
point on (?D with fewer violations. When ei is not a boundary edge call this neighboring 
triangle ti and when it is a boundary edge stop this process. If we have not stopped let 
62 be another edge bounding a zero angle slot in ti and stop if it is a boundary edge. If 
it is not a boundary edge then there are two possibilities. If /c*^ < repeat the above 
procedure letting ei play the role of eo and 62 the role of 6i and constructing an 63 in a 
triangle t2- If k^^{x) = conformally change x to 

X + eWei + 6^62 • 

Notice no triangle with k^{x) = can have two zeros by fact |2|, so for the initial zero 
violation to exist there most be a zero on the — e side of ewe2- Once again we have 
determined an 63 and t2- 
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Using this procedure to make our decisions we may continue this process form- 
ing a set of edges {ej} with the angle between and ei-\-i, A^''^^^(x), always equal to zero. 
Since there arc a finite number of edges either we stop at a boundary edge or eventually 
in this sequence will have some k < I such that = e/ and e^^i = ei^i. (This by the 
pigeon hole principle since some edge e will appear an infinite number of times in this 
list and among its infinite neighbors there must be a repeat). 

In the case the sequence never stops we can produce a contradiction. To do it 
first note if Cj and e^+i are in ti then ^*'*+^(a;) = i/j^I + V'tj^^- So for the set of edges 
{eiltfe we have 

l-l l-l 

i=k i=k 

our needed contradiction. 

In the case the sequence did hit the boundary perform the construction in the 
opposite direction. If we don't stop in this direction we arrive at the same contradiction. 
If we did then this computation still produces a contradiction on the path with the two 
boundary edges, since for a boundary edge in the triangle t we have ipfix) = '4'^{x) G 

(o,f)- 

q.e.d 

Lemma 8 If a stable x satisfies the condition that if a^(x) = then ai is in a triangle 
t with k^{x) = 0, then ^{x) is not in D. 

Proof: In this case, in order for x to be a boundary point of D for some t we have that 
/c* = 0. We will be looking at the nonempty set of all triangles with fc* = 0, Z. The 
first observation needed about Z is that it is not all of M and has a non-empty internal 
boundary (meaning dZ — dM). To see this note 

^ k\x) = ^A' -ttF = irdV + 27r(F - dV) - ttF 

= 2TrV - (ndV + SttF) + 27rF = 2ttV - 2'kE + 2'kF = 27rx(M) < 0, 
so there is negative curvature somewhere. 



46 



By the stability of x once again there can be no conformal transformation capa- 
ble of moving negative curvature into this set. Suppose we are at an internal boundary 
eo edge of Z, call the triangle on the Z side of the boundary edge to and the triangle 
on the non-boundary edge t_i. Since t_i has negative curvature the obstruction to the 
eweg transformation being able to move curvature out of Z must be due to to- In order 
for to to protect against this there must be zero along eo on the to side. 

Now we will continue the attempt to suck curvature out with a curvature vac- 
uum. Such a vacuum is an element of C indexed by a set of Z edges. The key observation 
in forming this vacuum is once again fact Q telling us if an angle in t is zero and k^{x) = 
then there is only one zero angle in t. Let ei be the other edge sharing the unique zero 
angle along eo in to and if ei is another boundary edge we stop. If ei is not a boundary 
edge use We^ + Weo to continue the effort to remove curvature. Continuing this process 
forms a completely determined set of edges and triangles, {ej} and {tj}, and a sequence 
of conformal transformations e Y27=o G C. 

We will now get some control over this vacuum. Note a vacuum never hits itself 
since if there is a first pair k < I such that tk = ti then t^ would have to have to have 
two zero and zero curvature, which fact ^ assures us is impossible. So any vacuum hits 
a boundary edge or pokes through Z into Z'^. 

In fact with this argument we can arrive at the considerablely sponger fact that 
two vacuums can never even share an edge. To see this call a vacuum's side boundary 
any edge of a triangle in the vacuum facing a zero. Now simply note if the intersection 
of two vacuums contains an edge then it contains a first edge e^ with respect to one of 
the vacuums. There are two possibilities for this edge. One is that tj+i has two zeros 
and k^{x) = 0, which we showed was impossible in the previous paragraph. The other 
is that Ci is a side boundary of both vacuums. In this case we have an edge facing zero 
angles in both directions in triangle with zero curvature, so this would force tp'^{x) = vr, 
a contradiction. So either case is impossible, and indeed no distinct vacuums share an 
edge. 

Let S be the removal from Z of all these vacuums. First I'd like to note that 
S is non-empty. Note every vacuum has side boundary. Since vacuums cannot intersect 
themselves or share edges with distinct vacuums, S would be nonempty if side boundary 
had to be in Z's interior. Look at any side boundary edge e of a fixed vacuum. Note 
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e cannot be on dZ — dM since then the vacuum triangle it belonged to would have at 
least two zeros and A;*(x) = 0. Furthermore e cannot be on dM since then ip^ix) = ^. 
So indeed S is nonempty. 

Now lets observe the following formula. 

Formula 10 Given a set of triangles S 

{eeS} e&dS-dM t&S ^ ^ e&dS-dM 

with the t in ipf{x) term being the triangle on the non-S side of e. 
Proof: 

{eeS} eeS 
eeS-dM eedS 
teS ^ ^ eedS-dM 

= E(-^)- E 

te5 ^ eedS-dM 

q.e.d 

Now every edge in dS — dM faces a zero on its S"^ side in a triangle with 
= 0, so 

eedS-dM 

Similarly each triangle has zero curvature so from the above formula we have 

E ^'(^) = i^i'^ 

{eeS} 

violating condition (n2). So we have constructed a violation to (n2) and ^(x) cannot be 
in D as need. 



48 



q.e.d 

It is worth noting that nothing prevents us from extending the main theorem of 
this section (and the next) to the case were p^{x) 7^ 27r and in particular to the version of 
this theorem where (x) = and the resuhing uniform surface is a finite area hyperbohc 
surface with cusps. 

2.3 The Ideal Thurston- Andreev Theorem 

In this section we prove and sate the general ideal convex Thurston- Andreev 
theorem, which is simply the polygonal decomposition case of the theorem of the previous 
section. The proof here is done in detail when the surface has no boundary, and dealing 
with the boundary can be accomplished exactly as in the previous section. 

2.3.1 The Statement and Reduction to Linear Algebra 

Corollary Q tells us that an ideal disk pattern is always unique when it exists, 
and we are now left to deal with the dilemma of finding good existence criteria. In fact 
in the convex case, (when p G (0, vr)^ relative to P) we will produce linear conditions on 
the possible p which are necessary and sufficient for the p relative to P to be the data 
of an ideal disk pattern. 

The first example of a necessary condition on p is the condition related to the 
fact that the angles at the vertex in a geometric triangulation sum to 27r, namely 

(ni) Y.^rip)) = 2vr. 

This and all the mentioned necessary condition will be demonstrated as such in the next 
section. Another necessary condition is a global one (though often localizable) , namely 
an insistence that for every set 5 of |5| polygons in P that 

(n2) j;r(p)>7r|5|. 

eeS 

With these two necessary condition we have our first pattern existence theorem: 
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Figure 2.13: A snake and a balloon 

Theorem 9 That (ni), (^2), and p G (0,7r)^ hold is necessary and sufficient for p 
relative to P to be associated to a unique convex ideal disk pattern. 



Note 2 This theorem works equally well for surfaces with boundary. The only modifi- 
cations is the obvious one that at a boundary vertex Y2e<^vi'^'^ip)) 



TT. 



Such pattern data is extremely common; for example the data determined by the circum- 
scribing circles in any random Delaunay triangulation of a varying negative curvature 



surface (see section 1.2). This fact is the main reason I've chosen to use the terminology 
of the disk pattern construction rather than polyhedra construction in this section. 

One way to prove theorem |9| allows some understanding of the non-convex case 
as well. Namely we will find some necessary criteria on angles when p G (—11,71)^. 
In order to articulate these conditions we need certain snake and a loop concepts in a 
triangular decomposition. 

Definition 13 A snake is a finite directed sequence of edges {cj}'^^ directed in the 
following sense: if k < I we start with the edge et between t^-i and tk, then we require 
Ck+i to he one of the remaining edges on t^ . Then letting t^+i be the other face associated 
to Cfc+i we require ek+2 to he one of the other edges of tk+i and so on until some tail 
edge ei and tail face ti are reached, and if I < k we reverse the procedure and add rather 



than subtract from the index. See figure 2. IS for examples. A loop is a snake {ej}^^; 



where Ck = ei and tk = ti, see figure \2. 14\ for a pair of examples. 
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Figure 2.14: A loop and a barbell 

It is a condition on snakes and loops which allows one to articulate the remaining 
necessary conditions. However as defined there are then an infinite number of such 
objects and it nice to first isolate a finite sub-set that does the job. 

Definition 14 A set of edges {ei}'^j. is called embedded if ei ^ ej. A snake {e^}^ is 
said to double back on itself if we have a pair of non-empty sub-snakes with {ci}^ and 
{^i}k-m containing the same edges. A barbell is a loop which doubles back on itself and 
such that {ei}'^^/{ej}"^^ is embedded. A balloon is a snake which doubles back on itself 
with {ei}\^i^/ {ei}2^^ embedded and such that ei = e^. 

With this terminology the remaining necessary conditions are 

(ns) Yli=k ^iip) < 1^ ~ when {ei}[^i^ is an embedded loop or barbell, 

and 

('^4) Yli=k^iiP) < (1^ ~ ^1 + ^)'^ when {ei}[^i^ is an embedded sake or balloon. 
With these conditions let 

N = {p E (— TT, tt)^ I p satisfies (nj) for each i}. 



To each ideal disk pattern we may produce an uniform element of a N by choos- 
ing a geodesic P* associated to the patterns circumscribing P. So the above necessary 
conditions would follow if ^'(N) C N. In fact in the next section we shall prove.... 
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Theorem 10 

^'(N) = N 

With this resuh we are in a position to prove theorem ^. 

Poof of theorem ^ from theorem 10: First we will show that p £ N relative 



to a chosen P* (where we assume polygons have been triangulated as in figure plop . To 
do this we need that conditions (n^) and (n^) are satisfied. Since 9 G (0,vr] (71,4) is 
automatic and (77-3) can only be false if there is a loop or barbell {ei} on which 0*^* = vr. 
Note that by our choice of P* (though any other choice would in fact still work with 
a slight modification) we see that the only snakes {ej} containing all 9^'- = it edges are 
snakes with edges contained in some polygon's fan and in particular can never loop up 
or form a barbell. So by theorem |l^ we see the point p described in theorem |9| has a 
preimage which intersects N non-trivially. 

At this point all we need is that a point y in this preimage y + C satisfies the 
conditions of theorem 0. Namely we will suppose that {y + C) B is not empty, i.e. 
{y + w) £ B with w £ C, and produce a contradiction. In particular this assumption 
gaurentees there is some triangle to with (i*°(y + w) = {7r,0, 0}. Let the edge ei of to 
be the edge with the {0,0} of on it. The fact d*"(y + w) = {7r,0,0} allows us to 
control the d*^ (y + w) data in the other triangle containing the edge ei. This follows by 
observing the natural and what will prove be useful decomposition of ip^'^ = ip^^ + ip^^ 



as in figure 2.15, and the following trivial but useful fact... 



Fact 3 When x E N we have ipt^i^) ^ (~^) f )? '3^'^^ when x G 9N we have ipfii^) ^ 

r —77 TtI 

L 2 ' 2J • 

reason for the fact: 

Let d*'(x) = {A,B,C} and note since B + C<A + B + C = l\x) < vr and 
^4 < vr we have 

vr ^ B + C-A B + C vr 

< < thf = < < -. 

2 2 ~ 2 ~ 2 2 

The second statement follows from the possibility of these inequalities becoming equali- 
ties. 
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Figure 2.15: The decomposition ip^ = ipf^ + 

Back to the proof. From this fact we have that 

e^^{y + w) = TT + ^ - iP'^iy + w), 

forces V^'Hy + ^) = f ) and hence d*^{y + tt;) = {Ajir — A,0} with the zero opposite 
to ei. Now the assumption we are in B imphes d^^{x) is not legal and A = vr or 0. 
Let 62 be the {0, 0} edge ti and continue this argument hence forming a snake of edges 
with O^"- = IT. Note by finiteness of the triangulation there must be a first / and k < I 
where ti = tf^. Note when this happens that forms an embedded loop with all 

its =71", contradicting condition (^3). 

So {y + C) B is indeed empty and we have our need triangulation, hence our 
needed ideal disk pattern. 

q.e.d 

2.3.2 Proof of Theorem |Tg 
Injectivity 

Our first goal in proving theorem |l^ is to show 

^(N) C N. 

For starters note the fact p G (— 7r,7r)'^ follows immediately from fact |l| in the 
previous section. Now we need to verify the conditions (ni) through (714) hold in ^'(N). 
Condition (ni) is equivalent to the following simple lemma. 
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Lemma 9 

^{V) = {peU^\ ei(p) = 27r}. 

{eiev} 

Proof: First recall that we know ^ is surjective, so we may express any p N as 
p = ^{x). L'So letting {e G v} denote the set of edges at a vertex v we have 



Y^e^{p) = Y,i^''{x)=p^x). 



So in particular the affine flat 

W = {peK^\ ej(p) = 27r} 

{eiSv} 

is precisely ^{V)- q.e.d 
Verifying condition (71-2) relies on the following formula. 

Formula 11 Given a set of triangles S 



with the t in ipfi^) tsrm being the triangle on the non-S side of e. 
Proof: 



{ee5} ee5 eeS 



t&S ^ ^ e&dS 



E 

te5 ^ eeds 



q.e.d 

Note for any point x G N that < and so with this and observation || we 

have Ylees(^ ~ ''P^i^)) > ^I'S'I and in particular condition (n2) is necessary. 
(71.3) and (714) rely on certain a pair of related formulae. 
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and that 



Formula 12 Let A*'*"*"^ be the angle slot between and Cj+i in a snake {ej}-^^. We 
have 

I l-l 
Y,0'^{x) = \l-k\n- Y.A^^'^\x) - ^11^ - i^ll 

i=k i=k 

and if {ej} is a loop 

l-l l-l 

^r»(x) = \i- -^A^''+^{x). 

i=k i=k 

Proof: Simply note both that 

and sum up. q.e.d 

The second formula immediately implies (ns) immediately, and from fact |l| we 
have x/jf > ^ allowing the first formula to demonstrate (n^). 

Surjectivity 

In this section we will finish the proof of theorem |lO| by showing ^' maps N onto 
A^. To do it let's assume the contrary allowing that that ^'(N) is strictly contained in 
and produce a contradiction. With this assumption we have a point p on the boundary 
of ^(N) inside N. Note p = ^{y) for some y £ (9N and that (C + y) P| N empty, since 
other wise for some w £ C we would have {y + w) € N hence forcing by the openness of 
^ p = ^'(y) = ^'(y + to be in the interior of ip{N). 

At this point we need to choose a particularly nice conformal version of y, which 
requires the notion of a stable boundary point of N. Before defining stability note since 
N is a convex set with hyper plane boundary if x € such that (x + C) P| N 7^ 0, then 
(x + C) Pi (9N is its self a convex k dimensional set. 

Definition 15 A point in x £ (9N is stable i/ (x + C) P| N = and x is in the interior 
of (x + C) Pi 9N as a k dimensional set. Any inequality forming N violated in order to 
make a stable x a boundary point will be called a violation. 
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The key property of a stable point is that a conformal change w £ C has 
x + ew ^ for all e > or for some sufficiently small e > we have x + ew must still be 
on 9N and experience exactly the same violations as x. The impossibility of any other 
phenomena when conformally changing a stable point is at the heart of the arguments 



in lemma 10 and lemma 11 below. At this point subjectivity would follow if for a stable 
X G 9N we knew that ^{x) could not be in N , contradicting the choice of ^ = ^{x) as 
needed. 

We will prove this by splitting up the possibilities into the two cases in lemma 



IC and lemma 11 



Before starting lets define ... 

Definition 16 The end of a snake {cj}-^^ is said to have a head with respect to x if 
d*''{x) = {0,0, tt} with the pair of zeros located at the angle slots of ti along ei. 



Lemma 10 If x £ dJSS is stable and a*(a;) = for ot in a triangle where k^{x) < 0, then 
^(x) is not in N . 

Proof: I will suppose that G N and produce a violation to the (n^) or (714) 

conditions. 

Now look at an angle slot which is zero in triangle to satisfying k^°{x) < 0. 
View this angle as living between the edges cq and ei. Note that in order for x to be 
stable that the ewe^ transformation (with its positive side in to) niust be protected by 
a zero on the — e side forcing the condition that ewei G , or else for small enough e 
we have x + eWei would be a conformally equivalent point on 5N with fewer violations. 
Call this neighboring triangle ti. If we see a pair of zeros and a tt facing the to from ti 
we stop. Otherwise let 62 be another edge bounding the zero angle slot in ti and repeat 
the above procedure if /c*^ < 0. If k^^{x) = form 

X + ewe^ + ewe2 ■ 

Note that the only way this construction could have difficulty is precisely the case in 
which ti was a head - in which case we already stopped. 

So we may continue this process forming a snake {cj} with = until we 

hit a head. Notice we can also make the same construction in the other direction. 
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Note as such if his snake formed a embedded loop, a barbell, an embedded 
snake with two heads or a balloon with a head then the ^*'*+^ = condition would 



contradict one of the formulas in formula 12, hence violating (n^) or (ns), and we would 



be done. It will be shown that one of these cases must occur. 

To produce the needed snakes note that by finiteness in the positive direction 
there is a first time when some tk = ti and A: < Z or we terminated at a head before 
such an over lap. If this sequence terminated in a head look at the snake in the negative 
direction and if it terminates in a head then we are done. If not we have the same 
situation as the positive snake not terminating in a head, i.e. there is a first when 
the one headed snake hits itself. If we hit the head we have our needed embedded loop. 
If not we have two possibilities either that e^+i can be chosen to be e/+i in which case 
we have our need embedded loop or Ck+i must be ei. In this case we can reverse the 
construction going form eo to e^to form the needed barbell. 

Suppose in the positive direction we experienced our first moment when for 
k < I tk = ti. Then as above if it hinges so ek+i can be chosen to be e^+i we have 
our embedded loop. If not we will form a chain in the opposite direction starting with 
6^ = 60. Now we are searching for the first m < n when e„ hits {ej} at or hits some or 
{ej}'^^ at tp or terminates in a head. If it terminates in a head we can form our needed 
balloon as 



{e,}t^lj{e,}U,lj{e,}^^,. 



If our snake hits itself at tm and e„+i can be chosen as im+i we once again get our 
needed embedded loop. If it hinges such that e„+i must then we can double back to 
form our needed barbell 



m=o[J{^^}lm-l[J{e^} 



k 

i=l- 



The other possibility is that {ei} hits first {ei} at tp with e„. As always there 
are two possibilities for how they hinge and in either case we can form one of the following 
embedded loops 



I 

i=p 
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or 

So in any case violation to (714) or (ns) can be produced. 

q.e.d 

Lemma 11 If a stable x satisfies the condition that if a^{x) = then ai is in a triangle 
t with k*{x) = 0, then ^{x) is not in N . 

Proof: In order to be a boundary point of N for some t we have that fc* = 0. We will 
be looking at the set of all triangles with = 0, Z, which is not all of M and has a 
non-empty boundary To see this note 

^ k\x) = ^A^ -ttF = 2'kV - SttF + 2ttF 

= - 2'kE + 27rF = 27rx(M) < 0, 

so there is negative curvature somewhere. 

By the stability of x once again there can be no conformal transformation 
capable of moving negative curvature into this set. Suppose we are at a boundary edge 
of call the triangle on the Z side of the boundary edge to and the triangle on the 
non-boundary edge Since t-i has negative curvature (and hence no tt angles) the 
obstruction to the ewe^ transformation being able to move curvature out of Z must be 
due to to- In order for to to protect against this there must be zero along cq in the to 
side. 

Now we will continue the attempt to suck curvature out with a curvature vac- 
uum. Such a vacuum is an element of C indexed by a snake. The first edge in the snake 
is the boundary edge e^. If in to eo faces a tt we stop. We say that we stopped at a head. 
If not let ei be the other edge sharing the unique zero angle along eo in to and if e\ is 
another boundary edge we stop. If e\ is not a boundary edge use to continue the 
effort to remove curvature. Continuing this process n steps forms a snake {e^}™^" and 

X + e ESl" € C. 

Suppose a vacuum hits its self and tn = tm- Then we must have an extra zero in 
tra in which case we have a d*"" (x) = {0, 0, tt}. Note it is not a head with respect to either 
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direction and it fact now form a vacuum loop. In the conformal change associated to this 
loop consistently changes the angles with value it to having value vr — 2e, a contradiction 
to stability. 

So any vacuum in fact pokes through Z into Z"^. In fact this argument shows us 
something slightly stronger, namely if a vacuum hits a triangle with = {0,0, vr} 

then it is a head. If not after we poke through we could still reduce the tt to vr — e 
conformally. This because under the lemma's hypothesis, there can be no zeros in Z^ 
protecting the vacuum form consistently sucking. 

From this note one vacuum can never pass through another since this would 
force a = {0,0, vr} triangle which is not a head for at least one of the vacuums. 

Similarly the outside edges of a vacuum always face zero angles in the vacuums and to 
be in A'^ two zeros can never face each other since then 

TT — TT TT — TT 

= — + — = 0. 

So all the edges associated to distinct vacuums are distinct. 

Now simply let S be the removal from Z of all these vacuums. If S is non- 
empty then every boundary edge of this set faces a zero in a triangle of zero curvature 



so formula receives all zeros from the boundary terms. Similarly each triangle having 

so 



exactly zero curvature gives us exactly a vr for each internal triangle in formula 11 
See5^'^(^) ~ '''"I'^l ^ needed to violate (^2). 

So we are reduced to seeing that 5 is nonempty. Since two vacuums can never 
border each other, this is reduced to seeing that every vacuum has a Z internal edge. 
Well suppose not then our vacuum would be an embedded snake with all boundary edges 
having a zero along them and all internal angles being zero and zero curvature. This 
forces our vacuum to have only triangles t with = {vr, 0, 0}, so could only be a pair 
of heads. When two heads face each other at edge e we have 6^ = 2vr contradicting that 
fact we are in A^. So S must be non-empty and we are done. 

q.e.d 



Chapter 3 

From the Discrete to the 
Continuous 



This chapter is dedicated to setting up the geometry and probabihty needed 
to compute the random variables discussed in section 1.2.2 of the introduction. In 
particular we fill in the details to all the steps in the probabilistic proof of the Gauss- 
Bonnet theorem sketch in section 1.1. In section 3.1 we develop all the geometric tools 
necessary to prove theorem [l| from the introduction. This includes section 3.3 where we 
examine some properties of Delaunay triangulations on surface independent of the rest 
of the thesis but of interest to anyone wanting get a feel for these triangulations. 

In section 3.2 we explore random Delaunay triangulation. In 3.2.1 we develop 
all the ideas need to prove the theorem |2| and Euler-Delaunay-Poisson formula. In 
sections 3.2.2 and 3.2.3| we develop a formula for computing random variables on the 



space of random Delaunay triangulations (or complexes), and in particular prove the 
Euler-Gauss-Bonnet-Delaunay formula form the introduction. Section 3.2.4 contains 
some particularly boring facts concerning the measurability of certain function and sets 
which arise in the first two sections of this chapter. 

3.1 Delaunay Triangulations 

This section is dedicated to the exploration of Delaunay triangulations. The 
technical backbone for all that occurs in this proof is theorem Q from the introduction, 
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which is dealt with in section 3.1.2. The geometry continues in section 3.1.3 where 
theorem |l] is proved (with the help of certain "inflating families" also dealt with in 
section 3.1.2). 



Then in section 3.1.5 we prove some interesting properties about Delaunay 
triangulations. These facts will not needed in the rest of thesis but are of interest in 
showing how certain facts about Euclidean Delaunay triangulations carry over to sur- 
faces. The facts explored include that Delaunay triangulations are local, several local 
facts, justification of algorithm constructions, as well as results showing that the Delau- 
nay triangulation of a dense p is in several ways optimal amongst dense triangulations; 
where a triangulation T is called dense is if each t £ T has its vertices on a ball of radius 
less than 5. For the Euclidean versions of essentially all these facts with quite different 
proofs see Throughout this section I'll assume 6 is simply min{|,r}. To entice 
the reader perhaps I'll mention now what these optimality properties are. For our first 
optimality property we have. 

Property 1 Among dense triangulations associated to a dense p the energy 

E{T) = ^ length{e) 

{eeT} 

is minimized precisely at the Delaunay triangulation. 

In section 4 we will see that the "gradient flow" of this energy tells us how to 
deform a dense triangulation associated to a dense p into its Delaunay triangulation. 

Notice the smaller this energy the squatter the triangles. That Delaunay trian- 
gles minimize this energy is one reflection of the fact that they prefer fat triangles. There 
are many realizations of this fact, another is the fact that the Delaunay triangulation 
attempts to minimize the sizes of associated spheres. To articulate this we first must 
acknowledge that by lemma |^ of the next section the vertices of a triangle t in a dense 
triangulation lie on a uniquely associated ball which will be denoted Bf. 

Property 2 Among all dense triangulations T associated to a dense p the Delaunay 
triangulation minimizes 

maxrad{T) = max{radius{Bt) \ t £ T}. 
In section 3.1.5| we will also see the sense in which this is locally true. 
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3.1.1 Some Geometric Reminders and Notation 

To understand this paper one must be aware of geodesies, the exponential map 
expp : TpM — M, and this map's imphcit interaction with balls and spheres. Let the ball 
of radius r at p be expp{Br{0)) - where -6^(0) is the open ball of radius r in = TpM, 
and denote it Bj.{p). Let the sphere Sj.{p) be its boundary. The first half of the needed 
results can be summed up in the following lemma. 

Lemma 12 (Geometric Reminders) Assume M is a compact Riemannian surface 
then: 

1. (Normal and Convex Neighborhoods) i and r are greater than zero. 

2. For any p we have that expp is diffeomorphism of Bi{0) onto Bi{p); and if p and 
q satisfy d{p, q) < i, then there is a unique geodesic of length less than i between 
them. 

3. (Gauss's Lemma) The unique unit speed geodesies from p to points in Bi{p) are 
given by 7(r) = expp ^'^jj^^ for some v; and any such geodesic is orthogonal to 
Sr{p). 

It is worth explicitly reminding the reader that given an orthonormal basis 
{ei, 62} at p G M we have the lovely normal coordinates: 

N{p, zi, Z2) = expp(ziei + 2362) : 5^(0) C ^ M. 

Sometimes it is useful to think in terms of angular and radial coordinates. Let v[9) = 
cos(^)ei + sin(^)e2, then as alternate coordinates we have the geodesic polar coordinates: 

G{p, r, 9) = expp{rv{e)) : (0, i) x ^ M. 

In order to vary p in the above it is necessary to have smoothly varying or- 
thonormal frames and they will be denoted / = {61,62}. We can always construct one 
on, say, a convex set; and can even globally have one on M — {points} f. This is ac- 
complished by Graham- Schmidting a pair of generic vector fields, where {points} f is 
the finite set of points where the vector fields are not linearly independent. In the the 
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presence of a frame we have a canonieal choice for a a ^ rotation field Q; by using the 
fact v{6) parameterizes the tangent spaces we can let 

e{v{e)) = v{9)^ = - sm{9)ei + cos(^)e2. 

The last frame idea used is that of a geodesic frame at p. Fixing an orthonor- 
mal basis {61,62} at p, let the geodesic frame be the the frame given by the parallel 
transport of this orthonormal basis of TpM along the geodesies spitting out from p. 

The other bit of geometry used are some basic Jacobi field results. Recall that 
a Jacobi field is a vector field along a geodesic 7(r) satisfying 

with the initial condition J(0) = V and -^(0) = W . (The choice of r here, as opposed 
to the usual t, stems from the fact that our Jacobi fields will be thought of as living 
along geodesies parameterized by r in some geodesic polar coordinates.) 

Jacobi fields have the wonderful property of being in 1-1 correspondence with 

smooth one parameter families of geodesies (in the standard notation J = ^ where 
r(r, s) is a geodesic for each s G [sq — 6, -Sq + e]) and r(r, sq) = 7(r) . On a surface 
they come in four flavors. To taste these flavors first one notes Jacobi fields with initial 
conditions perpendicular or parallel to ^(0) remain as such for all time. Also the 
equation is a second order O.D.E., hence linear in its initial conditions - so a Jacobi field 
can be decomposed into its component along 7 and its perpendicular component simply 
by decomposing its initial conditions as such. 

So we have all Jacobi field are linear combinations of the following types (and 
some examples of corresponding r(s,r)): 

• J(0) = ^(^o) and "^(0) = (The two parameter family corresponding to this case 
is the friendly r(s,r) = expp{{s + r)v{6o)).) 

• J{0)=v^{9o) and ^(0) = 
. J(0) = and ^(0) = v{9o) 

• J(0) = and -^(0) = v^{6o) (The two parameter family corresponding to this 
case is the friendly r(s, r) = expp{rv{s)). In particular this Jacobi field is precisely 
N4rv^{9o))=G. (^). 
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Sometimes one starts with a T(s,r) and wants to understand the associated 
field - a well known example, that will prove relevant to us, can be constructed by fixing 
a geodesic a{s) and a frame giving the vector field v{9){s) along a{s). The Jacobi field 
along expa[so)'''v{6) corresponding to expa{s){^v{9)) is the the one with initial conditions 

^(0) = S(o), ^(0) = 

Here are some facts we will be needing about Jacobi fields, and in particular 
what they look like in normal coordinates. 

Lemma 13 1. The Jacobi field along expp{rv{9Q)) with initial conditions J(0) = 
v{9o) and -§^(0) = in normal coordinates is J{r) = v{9o). 

2. The Jacobi field along expp{rv{9Q)) corresponding to J(0) = and ^ = V'^{9o) in 
normal coordinates is rv^{9o). Calling \ \rv-^{9o)\\M = j6»o(''); '"'fi have je^^ 's Taylor 
expansion is je^ir) = r(l — + 0{r^)). 



3. The Jacobi field along expp{rv{9Q)) with initial conditions J(0) = v-^{9q) and 
■^(0) = in normal coordinates can be written J{r) = hQg{r) . k V-^{9o). 

4- Here are a few immediate consequences of part two - the area of a ball at p function, 
a{r), satisfies a{r) = r^(7r — j|r^ + 0(r^)); and the product jQ^jo^je.^ =r^(7r — |r^ + 
o(r^)). 



It is necessary to have certain global estimates of the above o(r^) functions, 
resulting from the fact we are on a compact surface. 

Lemma 14 // M is compact, then there is a Cm > such that when r < i we have the 
O(r^) functions in 3 and 4 above all satisfying \0{r'^)\ < Cm^^ globally. 



Proof of lemma 14: First I'll find the constant related to 4(c) above. The 
Jacobi fields are the solutions to an o.d.e. - so vary continuously with initial data - which 
is indexed by {p,v) S UTM (the unit tangent bundle). In particular je^ir) = j{p,v,r) 
varies continuously with initial data. Now recalling from above the Taylor expansion in 
the radial variable at r = we have j{p,v,r) = r + lMeI _|_ rO{r^){p,v,r); and form 
Taylor's formula the third term in this sum can in fact be represented as (assuming the 
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metric is smooth) 

uip,v,r) = — {l-tf-^{v,p,tr)dt. 

In particular this term is continuous even after dividing by r^. Now UTM x [—i, i] is 
compact so we can feel free to take |^^^^^3^|'s maximum over this set for our Cm- 

Now observe that the other O(r^) are directly related to this one and the 
curvature function, and since the manifold is compact st/,p|A:| exists. Choose the Cm in 
the lemma to be the biggest of the constructed constants among these O(r^) functions. 

q.e.d 

Now we will prove some potentially less familiar geometric facts. The first of 
which will be lemma || from the introduction. My original proof of this fact was quite 
inelegant and can be found the the appendix. I'd like to thank Albert Nijenhuis for 
sharing his beautiful proof with me. 

The next facts we will need concerns the notion of an inflating family of circles 
through a pair of points {p,q}. As a set, this family will be all circles of radius r < 5, 
passing through both the points p and q. The following lemmas justify the fact that this 
set can be thought of as the continuous inflating family of circular balloons to the left 
or right "sides" of the geodesic through p and q, as in figure |3.l| . 

Too articulate this given a continuous curve Cpq{t) let Dpq{t) be the continuous 
family of closed disks centered at Cpq{t) of radius d{p,q) + \t\, and let dDpg{t) be the 
corresponding family of circles. For starters we have a lemma gaurenteeing the existence 
of inflating families. 

Lemma 15 For each pair of points p and q such that d{p, q) < 6 there is a curve 
Cpq '■ id{p, q) — 5,5 — d{p, q)) ^ M such that every circles of radius less 5 going through 
p and q is dDpq{t) for some t. 

Notice that the radius increases monotonically as |t| does. We need a lemma 
giving us another sense of monotonicity. To articulate it we first develop a little notation. 
If d{p, q) < i let pq be the unique shortest length geodesic connecting p and g, let inid{pq) 
be its midpoint, let Br{p) be the ball of radius r at p, and let dBr{p) be its boundary. 
To articulate the next lemma note that the geodesic connecting p and q removed form 
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Figure 3.1: The Inflating Family 




Figure 3.2: The Monotonicity of Area 

Bi{mid{pq)) decomposes Bi{mid{pq)) into two open sets, which wih be referred as the 
decomposition determined by pq. Furthermore since 5 < t any ball of radius less than 5 
with p and q on its boundary is also divided into two such pieces. 

Lemma 16 Cpq{—{5 — d{p.q)),0) is in one half of the decomposition determine by pq 
(call it ) and Cpq(0, 6 — d{pq) in the other (H^ ). If d{p, q) — 5<c<d<5 — d{q, q) 
then Dpq{c) fl C Dpq{d) fl and Dpg{d) f] C Dpg{c) f] with all the subsets 
proper. 



The essence of this lemma is that figure 



3.2 



is accurate. 
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3.1.2 Small Circle Intersection 

My original proof of lemma |l| was a bit long winded and I would like to thank 
Albert Nijenhuis for showing me the elegant proof presented here. Actually the proof 
presented here is a much less elegant modification of the one Nijenhuis showed me where 
the needed injectivity radius is chased through the argument. Any errors or any realiza- 
tion that the injectivity bound is not as sharp as possible is solely my fault. Note that 
the lemma will follow if one could show that two circles of radius less than | intersected 
in at most two points, which is proposition ^ below. Another proof of this can be found 
in [12 1 where the sharper bound of | is also demonstrated. 

If a and b have distance between them less than the injectivity radius, d{a, b) < 
i, denote as ab the unique minimal length geodesic connecting them. Let x and y be the 
centers of two intersecting circles of radius less than | . Let Cy be the circle centered at 
y, and xyi be the open geodesic segment containing xy with midpoint x of length i. 

Lemma 17 xyi intersects Cy in exactly two points. 

Proof: Note d{x, y) < j and any point on Cy has a distance less than ^ < § form x, 
so the diameter of Cy intersecting xyi is in fact included in it. So xyi intersects Cy at 
least twice. However a geodesic ray form y of length less than ^ < i can hit each circle 
of radius r < i centered at y only once, so xyi can only hit Cy twice. q.e.d 
Consider f : Cy ^ R, f{z) = d{x,z), so / measures the distance between x 
and the points of Cy. 

Lemma 18 If z,w G Cy, z ^ w, and f{z) = f{w), then f has (at least) one critical 
point on each of the two circular arcs zw. 

Proof: This follows from a standard min-max argument. q.e.d 

Lemma 19 If p is a critical point of f , then p £ xyi. 

proof: Since p is a critical point and has a distance at most ^ from x, the 
tangent to the segment xz at z is perpendicular to Cy. So by Gauss's lemma the distance 
to y along this geodesic or its potential | or less length extension is less than | form x, 
and we have p C xyi. q.e.d 
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Proposition 2 Any two distinct circles of radius less than | have at most 2 points in 
common. 



proof: Suppose not. By lemma 18 there would be at least 3 distinct critical 



points, which would all lie on xyi by lemma 1£. But lemma 17 assures there are only 
two such critical points, the need contradiction. q.e.d 



3.1.3 Inflating Families 

I will construct the family by first describing the point set along which the 
centers of the circles in the family live with a different parameterization than that of 
lemma Note a point set this curve consists of points satisfying d{p, z) — d{q, z) = 
- with d{p, z) with less than 5. We will always denote as pq the most sensible connected 
extension of pq. For the following lemma let it be the connected extension of pq in 
^i{p) U Bi_{q). With this notation we have: 

4 4 

Lemma 20 Ifd{p, (?) < I the point set described by d{p, z)—d{q, z) = in Bi_ (p) |J Bi_ (q) 

4 4 

can be described by a curve Cp^q{t) with t £ {c,d) satisfying 

1. c < < d, Cp^g(O) = mid{pq), and mid{pq) is the unique point of pq on Cp^q{t) 

2. d{cp^q{t),p) and d{cp^q{t),q) strictly increase as the parameter \t\ increases. 

Proof: To see that the point set is a nicely parameterized curve it is useful 
to note that it can be described as the integral curve of a vector field. Let Dp denote 
d{p, x) and let VDp denote its gradient. Note the solution to the equation Dp — Dq = 
are integral curves of the vector field QiVDp — VDq) where © is a | rotation field. 

To understand these integral curves we will first deal with the the uniqueness 
of mid{pq): suppose a point I ^ mid{pq) is on pq. Then I is within | of p; hence d{p, I) is 
determined by the length of the segment of from p to I, similarly for q (using |). Now 
note that as we move from mid{pq) toward, say, p that Dp decreases while Dq increase - 
so Dp — Dq ^ Q aX another point oi pq. When I is on pq/pq, say above p, the segment of 
pq from q to I in fact covers the shorter segment from p to I - forcing Dp — Dq ^ once 
again. So / cannot satisfy Dp — Dq = 0, forcing mid{pq) to indeed be the unique point 
of pq on Dp — Dq = 0. 
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Now we we will see that we indeed get a union of curves by noting that that 
the vector field has no zeros in this set. In fact the triangle inequality tells Bi{p) C 

4 

Bi_ (p) n Bi_ (q) and in this region we will prove the stronger fact that V-D„ 7^ cVDq for 

22 

any c. First note that VDp is unit length, with integral curves the geodesies emanating 
from p. So at mid{pq) we have Q{VDp — VDq) is length 2. To finish the assertion assume 
at some point p not on 7^^^ that we have VDp = c\JDq. First note from the that fact 
that VDp is unit length we have c = ±1. There are two cases, first we'll deal with c = 1. 
Since the geodesies satisfy a second order O.D.E they are uniquely determined by their 
position and tangent vector, so when c = 1 we have both the geodesic from p and the 
geodesic form q are the same curves. Without loss of generality p is further away than q 
and this point lies along the same minimal length geodesic (of length less than |) which 
connects p and g, i.e. pq. In the case c = —1 we can follow the geodesic form p to the 
point and then from the point back to q forming a geodesic of length less than i - which 
then must by the definition of the injectivity radius be the unique such one, i.e. pq. 

To finish off the first part we need that our curve has only one component. This 
is intimately related to the second part. To see why we first look at the component of 

Cp,q{t) in Bi_{p) and note any component of Dp — Dq = would have to have a point 

' 4 

closest to p. This closest point is tangent to a sphere emanating form p. The same sort 
of phenomena must take place for the distance function to have a critical point; namely 
if a point z along any integral curve of Q{SJDp — VDq) is a critical point of the distance 
function D(p,-) then either Wp = VDg or a circle is tangent to the solution curve. In 
the tangent case Q{VDp — VDq) = cQVDp, or rather VDp — VDq = cVDp-, so both 
these situation have forced the case VDp = cVDq. so we may use the above observation 
to note that the point where this occurs is on pq; but from above to be on Dp — Dq = 
and pq means you must be exactly mid{pq). So we have both that every component of 

Dp — Dq = {) m. Bi_ (p) contains niid{pq) , and that the distance to p parameterized by t 

4 

can have no critical points except at mid{pq) (similarly for q). 

q.e.d 



We can now prove lemma 15 by noting by part (b) of the above lemma we 
may reparameterize as claimed and that any circle of radius less than 5 has its center 

contained in Bi. (p) \J Bj^ (q) so this new parameterization can be chosen on (and beyond) 

4 4 

{d{p,q) -6,6 - d{q,q)). 
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Now we shall prove lemma 16 



Proof of lemma |l^ : Notice the first part follows from part (a) of the above 
lemma. To prove the rest of it it is useful to isolate a sub-lemma. 

Sub- lemma 1 The containments between the halves must switch. Precisely ifDpq(c) P| C 
Dpq{d)f]H^ then Dpg{d)f]H^ C Dpg{c)f]H^ and visa versa. 

Proof: Suppose that Dpg{c) f]H+ C Dpq{d) f] H+ . Note that either Dpq{d) ^H' d 
Dpq{c)f]H~ or Dpq{c)f]H^ C Dpq{d)f]H~ since a violation of this inclusion would 
result in a third intersection of two circles of radius less than 5 - contradicting lemma 
|l[ From this observation, to violate the above choices would mean that Dpq{c) f] C 
Dpq{d) Pi H~ . Such an inclusion would force the circles to be tangent at there intersec- 
tion points - and hence the centers of both disks to be on pq (since the curve orthogonal 
to the tangent is a geodesic heading to the circle's center by Gauss's lemma). But then 
d{p,ci) = d{p,C2) forcing both the centers and the radii to be the same. So the disks 
would be identical contradicting distinctness. 

q.e.d 



We may finish the proof of the monotonicity of infiation half of lemma [1^. By 
the above sub-lemma we are left to explore three cases. 

The first case is where c or d is zero. If c = then the geodesic from mid{pq) 
to Cp^q{d) to the boundary of Spq{d) is strictly larger than d(mid{pq),p) (by part (b) 
of lemma So a switch of containment is impossible in this case. For the d = 
case note by the above sub-lemma one is contained in the other, and as just noted 
Dpq{{)){^H~ C Dpq{c){^H-. So it must be the case Dpq{c){^H+ C Dpq{{)){^H+ as 
needed. 

If < d then from the above we know Dpq{{}){^H^ C Dpq{d){^H^ . From 
this the fact Dpq{d){^H^ C Dpq[c)^H'^ would by the continuity of Dpq{t) and the 
uniqueness of the circles in the pencil (lemma pO| ) force there to be a < / < c (or 
c < / < 0) such that dDpq{f) f^dDpq{d) contains a pont in H~^. This produces extra 
intersections of distinct circles contradicting lemma |l[ 

To handle d < recall form the first case that we have both Dpq{0) P| C 
Dpq{c) f] and Dpq{0) f] C Dpq{d) . So by the sub- lemma we must have 
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Figure 3.3: That Which Cannot Occur 

Dpg{c)f]H+ C Dpg{0)f]H+ and Dpg{d)f]H+ C Dpg{0)f]H+. Now we can contra- 
dict Dpg{d)f]H^ C Dpg{c)f]H^ using same argument as in the second case, with 
the one variation being that this time we construct an / such that d < f < where 
dDpg{f) Pi Pi dDpg{c) P contains other points. 

q.e.d 

3.1.4 The Existence of Delaunay Triangulations 

Here we prove theorem |l] from the introduction. To do so first we needed to 
actually construct R. Part of i?'s construction is canonical; namely the mapping of the 
1-skeleton. This because \l^{pi,...,p„}\ is in an affine space and we can let the edges of 
l-^{pi,...,p„}l onto there corresponding unique geodesic segments by factoring with an 
affine map through the geodesic's unit speed parameterization. Now we need to extend 
this continuous mapping of the 1-skeleton to a continuous map of |-?^{pi,...,p„}|. 

Let q = {qi}f^i lie on a disk D of radius less than 6 and note the existence of 
R requires only an identification of the correct triangle qiq2q3- Throughout this proof a 
bold face letter will always denote a triple forming a face. 

Since D's radius is less than r for each pair of distinct points a and b on dD we 
have that D is split into its two halves. If q lies in one half let Dg^ denote the closure 
of that half, the region enclosed by the bold dashed line in figure |3.3| . Notice since the 
radius is less than i that by the uniqueness of small geodesic this set is convex. 
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Figure 3.4: A Missing Region and Balloon Popping 



Lemma 21 //q lies on a disk D of radius less than min{T, |}, then the triangle qiq2Q3 = 
Pi £)^ii3 Pi /j^ags ^ convex topological disk bounded by qiq2, qiq3 and 9293- Further 
more for any a and b on dD i/q C D'^^ then qiq2q3 C D^^ . 

Proof: Since (?i(?2'?3 is the intersection of convex sets it is convex. Each 
contains q so by convexity each contains all the qiqj. Now the only possible boundary of 
this set is the q^qj or the boundary of D, but other than q all the points on L>'s boundary 
have been eliminated. So 91^293 is a convex set with boundary {qiqj} contained in the 
convex disk D, hence itself a topological disk as needed. 

Now for the second assertion. qiq2q3 is an embedded disk so if the assertion 
fails then one of the qiqj crosses ab. By continuity it will hit ab in at least a pair of points, 



as in figure |3.3| . So this pair of points is connected by two distinct geodesies of length 
less than i. Hence we have a contradiction regarding the the uniqueness of geodesies of 
length less than the injectivity radius. 

q.e.d 

Now that we have the mapping R let us prove R onto. Since R{\K^p-^^ p^y\) is 
closed (a union of closed sets) if it is not onto it must miss an open set, and the boundary 
of this open set must be composed of the edges of triangles in (see bold region 

in figure 3.4). So we would be done if the notion of edge and the notion of an edge 
belonging to a pair of faces becomes interchangeable in the presence of 5-density. 

To see this suppose we have an edge connecting a pair p and q, then by its 
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very definition there is a k £ {d{p,q) — 6,5 — d{p,q)) with the property intDpqD{k) is 
empty of points. We may right off the bat use our assumption of (5-density to see that 
this empty disk must correspond to a radius d{p,Cp^q{k)) < 6. Now start inflating to the 
left of k. When t = d{p, q) — 5 the radius is 6 so intDpq{d{p, q) — 5) contains a point. By 



the monotonicity lemma (lemma 16) once you hit a point moving leftward you cover for 
it all future time, so there is a unique c G {d{p, q) — 5,5 — d{p, q)) such that dDpq{c) first 
contains a third point to the left of k. You may view this process as the blowing up of 



a circular balloon as dynamically represented in figure 3.4, and note we proved that it 
must pop. Similarly there is a c? > /c to the right were the balloon pops. Since we are 
assuming that there are never four points on a circle we have c < d. So dDpq{c) and 
dDpq{d) correspond to the unique left and right faces as needed. 

Now let us deal with R being 1-1. Since the individual faces are embedded this 
would follow immediately from the following lemma. 

Lemma 22 Two distinct faces P1P2P3 and qiq2qz can intersect only in a vertex or an 
edge. 

Proof: Suppose P1P2P3 and 91^293 intersect and have associated now intersecting circles 
X and Y respectively. If X and Y intersect in one point then p and q being the only 
points oi P1P2P3 and qiq2q3 on dX and dY respectively forces P1P2P3 Cl Q1Q2Q3 to be a 
vertex. So we may assume by lemma |l] that X and Y intersect in precisely two distinct 
a and b. Furthermore as in lemma ^ X and Y are both split by ab and since no disk 
contains four points the interiors of Xp^ and Yq^ are on opposite sides of ab. So by 
the second part of lemma ^ P1P2P3 and qiq2q3 can only intersect along ab. Since the 
boundary of the say qiq2q3 is {qiqj} for this to occur either qiq2q3 H o& is precisely a or 
b, some q^qj is tangent to ab, or some qiqj contains two points of ab. By the uniqueness 
of small geodesies in the last case q^qj = ab. In the tangent case the fact that geodesic 
are the solutions to a second order O.D.E. once again gives us qiqj = ab. So in any case 
distinct faces can only intersect by sharing a vertex or an edge as needed. 

q.e.d 

So we have i? is a bijective continuous map from a compact space hence a 
homeomorphism, and theorem |l| has been proved. 

q.e.d (theorem ||) 



73 



3.1.5 Basic Properties of Delaunay Triangulations 

The first property worth exploring is the fact that a Delaunay triangulation is 
a local phenomena. 

Definition 17 Call a dense triangulation Delaunay at an edge e if the vertex forming 
the face to left side of e is out side the circle associated the right side 's face. 



Notice by the monotonicity lemma (lemma 16) that the property of being De- 



launay at e is equivalent to the same property with the sides reversed. 

Lemma 23 A triangulation is Delaunay at each edge if and only if it is Delaunay. 

Proof: Clearly Delaunay implies locally Delaunay. To see the converse suppose it were 
not true and there is a triangle t with an extra point p in its associated disk. Then p is 
in a region of the disk to one side of an edge of t. However since each edge is Delaunay 
this vertex cannot be the third vertex of the face on this side or be in its associated 
face. Further more by the monotonicity lemma it is also in this new triangle's associated 
disk. So the same situation persists for this new face. Using this observation one can 
now construct a sequence of such triangles and edges with each new edge clearly closer 
to the point. So we produce an infinite sequence of distinct edges containing points in a 
neighborhood of p contradicting even local finiteness of the triangulation (let alone the 
fact it is globally finite.) 

q.e.d 

There is another very basic local property of the dense triangulations. I will 
call it the switching property. 

Lemma 24 // two triangles abc and abd of a dense triangulation are not Delaunay at 
ah then abed is convex. The triangulation formed by switching the diagonal inside abed 
is Delaunay at cd and length{cd) < length{ab) . 

Proof: First we prove the convexity assertion. Note the triangles both are inscribed 
in Babe n Babd which is convex and have interiors on opposite sides of ab. The key fact 
is that a shortest length geodesic connecting two points in abc and abd must cross ab. 
From this observation abed is convex, since if it were not then then we could construct 
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a shortest length geodesic form an interior point of a triangle to ab which must cross 
through another triangle side, contradicting the fact that the triangle is convex. 



From the monotonicity lemma (lemma 16) the ball B with diagonal ab is con- 



tained in BabclJ Babd and contains BabcClBabd, see figure 3^. Furthermore lemma 16 



gives us that the closures of Babe Pi ^abd and B intersect only at a and b. In particular 
the ball with diagonal cd is strictly contained in B with a and b now out side and on 
different sides of cd. So we immediately see the the length of cd is less than the length 
of ab, and we may use the inflating family associated to c and d to hunt down a and b 
demonstrating that the triangulation is indeed Delaunay at cd. 

q.e.d 

It worth noting that the converse holds among convex Delaunay parallelograms, 
although the convexity assumption is now necessary. 

The first corollary of the locality and switching lemmas is the proof of property 

B 

Proof of Property |l|: Note the set of possible triangulations is finite (it is 
certainly less than or equal to the cardinality of the vertices choose 3), so there is at least 



one minimizer of E. If the minimizer is not D then by lemma 23 there is an edge which is 
not Delaunay, so by lemma ^ we may rearrange the triangulation at the parallelogram 
living at this edge. Note the only edge which is changed is the diagonal and its length 
is strictly decreased hence the energy is decreased contradicting T's minimality. So the 
minimizer must he D. 

q.e.d 

To see property ^ it is useful to first prove the following local lemma. 



75 



Lemma 25 Let Dbet^a triangle in the Delaunay triangulation and let be the center 
of t^ 's associated circle. If T is a dense triangulation relative to the same points and t 
is a triangle in T such that both Cd G t, then rad{Bt) > rad{Bd)- 

Proof: Assume t^ ^ t, since the result is trivial otherwise. First note all of 
the points including the vertices of t are in B^. By the triangle inequality if B^ is a 
proper subset of Bt we are done. So we may assume the bounding circles, dBf and 
dBt^, intersect in two points (lemma |l]) with lemma |l| assuring us the vertices of t all on 
one side of the decomposition determined by these circles (or this side's closure). From 
lemma ^ t and hence q will also lies on this same side side of the decomposition (or its 
closure). Now look at the inflating family determined by the intersection points of the 
bounding circles. By the monotonicity lemma (lemma [T^) in order to contain greater or 
equal area on the side of the pencil where the vertices of t and Cd live rad{Bt) must be 
greater than rad{Bd) as needed. 

q.e.d 

Proof of Property [2]: First take a triangle td ^ D with Bd having maximal 
radius in D. For any dense triangulation T there is a triangle t of T containing the center 
of Bd- Now lemma ^ guarantees that the radius of Bt is at least that of Bd proving the 
corollary. q.e.d 

Its also worth noting that the geometrically minimal spanning trees (GMST) 
make sense on a surface as well as relative neighborhood graphs (RN) and Gabriel graphs 
(GG) and for the same exact reasons as in Euclidean space we have 

GMST c RNG CGGCD. 

So to construct all these things it would be nice to see that we can construct D com- 
putationally, and that in fact the basic algorithm for changing a triangulation T into D 
works just fine. 

The algorithm is essentially gradient flow of the energy E; namely switch edges 



which are not Delaunay as in lemma 3.5. This procedure terminates in certainly fewer 
than the number possible triangulation steps since each triangulation has an associated 
energy and by lemma 13.51 the energy is always decreasing under this "flow" . By lemma 



23| this flow must terminate in the Delaunay triangulation. 
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In fact this procedure must end before the number of vertices choose two steps. 
This fohows immediately from the following lemma telling us that an edge once switched 
can never be formed again. 

Lemma 26 // one dense triangulation contains ah and fails to he Delaunay at ah then 
there is no dense triangulation containing ah which is Delaunay at ah. 

Proof: This follows from the monotonicity lemma (lemma [l^ ) which guarantees that 
no disks in ab's inflating family can ever be empty. q.e.d 



3.2 Random Delaunay Triangulations 

As described in the chapters introduction here we tackle a the details needed 
about random Delaunay triangulations and the needed random variable computations. 

3.2.1 Basic Facts About Random Delaunay Triangulations 

Poisson Point Process Reminders: In order to set up some notation and con- 
vince the unfamiliar that nothing deep is occurring, we will now construct from scratch 
what little we need of the Poisson point process, i.e. the probability space from the 
introduction. 

To get started it is useful to look at the space of sets of ordered sets of n points, 
with a small set conveniently removed. This small set is the union of three closed measure 
zero sets: the set where some Xi = xj for i ^ j, the set where four points land on a circle 
of radius r < min{|,r}, and the set were three points land on a circle of radius exactly 
5. (These removals are merely a technical convenience - so I will not index P with a. 6 - 
though a bit of 6 has been programmed into it.) When I say measure zero I mean using 
the Riemannian product measure dA""; that these sets are closed and measure zero is 



straight forward though bit boring; detailed proofs can be found in section 3.2.4. Let 
x"M_ denote this full measure open subset of x"M, and let dA^ be the Riemannian 
volume element restricted to this open set. The sets of points of interest to us can now 
be expressed as points in P = Unez+ 

The measure, P'^, on this space is given by weighting dA^ on each component 
by ^e~^^. This is a probability measure since the measure of x"M_ under dA^ is A^, 
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and so the size of P is Yl'^=o ^nf = e~^^e^^ = 1. Lastly, the measurable sets B 

will be the Borel fi-algebra. 

On the index we have what is usually referred to as the Poisson distribution, 
namely P^(x"M_) = e~^-^ ; which justifies (in this model) the computation of 
Ea(F) = AA from the introduction. As further warm up from this view point it is useful 
to explore the probability that some chunk U of area A(U) in M is empty of points. 
Strictly speaking of course x G P is not a set of points in M. Being careful about is 
perhaps less burden then any confusion resulting form it, so I will introduce the mapping 
set : x"M_ — > 2*^ defined as set{pi, . . . ,pn) = {pi, ■ ■ ■ ,Pn} C M. So we are trying to 
find the size of the set where set{x) D B = 0,or rather IJnez+ {^^U^ f] P). For each for 
each n using the Riemannian volume element this set has precisely the measure of x'^U^, 
which has size {A — A(U))^ ^-^^ — ; and now we can sum them up to find the needed 
measure is 

n=0 

So we arrive at the last property of this model need for this proof. 

Now we would like to beginning proving theorem || from the introduction. We 
start with: 

Proposition 3 E;^(1t5) > 1 — ce'*'*', with c and d greater than zero. 

Proof: The heart of the proof is that points when distributed as above land 
very densely which is what being in D5 means; and we know by theorem ^ that D5 C T5. 
Note by the triangle inequality any ball of radius | covering the center of a ball of radius 
5 must be contained in it. So to force every ball of radius 8 to contain a point it is 
sufficient to cover the surface with balls of radius | , and then force this finite set of balls 
to contain points. So cover the surface with the | balls and pick out a finite sub-cover 
with say c elements. Let d be the minimum area amongst these c balls. Now let E5 C P 
be the set where none of these c balls is empty. By the above observation C D5, so 
we have C C E^. Now E^ is precisely the union of the sets were some individual 
of the c balls is empty - so clearly measurable. By the sub-additivity of measures, the 
warm up computation, and the choices of c and d we now have E^ has size less than 
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^^-dA\^ So assuming is measurable (in fact it is open, to be seen in the final section), 
we have < ce~*^^^, as needed. 

q.e.d 

From this we can get the needed estimates regarding the F, V, and E finishing 
of theorem |2| from the introduction (up the measurability which is in fact trivial since 
the functions are continuous on P - a fact whose details can be found in section p.2.4 ). 



Proof of Theorem ^: First let's do it for the faces. If there are n points then 
the number of faces is usually much less than (g) . So by the above proposition even the 
worse case for the area of is better than when ce~'^'^^ worth of area is crammed into 
the part of the space where (3) is the largest (or (2) for the edges or n for the vertices). 
In other words letting S = {B C B : E;s^(1b) = ce^'^'^} we have 

E.(QlT,)<supE.((^)lB 

So really it is this quantity we estimate. The monotonicity of (3) in n indicates that we 
can realize this sup with any set size ce~'^^ which fills up all the x"M for n > Nx along 
with some subset of size ce~'^^ — m e~^^ A"X" -^^ x^^M. To make sure we can 
perform this construction it is useful to observes as a scholium to the above proposition 
that we have d < A. So we naturally may stay away from the e~'^^ point mass at x'^M, 
and can indeed realize a set of the needed size. So for this Nx we have 



A—/ T7.I 



and 



ni 

n=Nx+l 



E 



n! 

n=Nx 



So we now have 



E(F|CS)<«< ± (3)e--^ 



- ^ ^ (n-3)! ~~6~ ^ ^ 

n=Nx ^ ' n=Nx-3 
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- 3! ^ + 6 

The second term decays faster than any polynomial, so we are reduced to 
seeing the terms in the form (^A) e~^ ^^j^'^_^y decay quickly, with / and k non-negative 
integers. 

To approach these terms first we use that our upper bound gives us (jv^_^_^)! < 
ceiA-d)X^ or rather AX < {{Nx + l)!ce(^-'^)^)^. 
In particular 



< 



Applying this to our term gives the quite manageable estimate 

(AA)^e-^^^^^ < (iVA + iy^\AX)'e-''\ 

Now we see these terms would decay faster than any polynomial if the Nx grew 
(with ^A) no faster than a polynomial. To see this is so, note that our estimate in the 
other direction gives us Y1'^=Nx ^^nT" ^ ce^"^"*^^^. Now d < A so we might hope must 
remain quite small for this to be so large. In fact this is easy to show; we may even 
observe a rather extreme fact that Nx cannot have a subsequence grow even as fast as 
(^A)^ with 



lim > ; — = lim > — - = 

A=^-»oo n! m-»oo ^-^ n! 

. _ 1 

To see this note that Stirling's formula tells us that n! > ^n"e "'(27rn)2, so 
assuming m > e we have 



lim 

m— >oo 



— < C lim > (27rn) 2 
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<c lim y -] <c lim y (-) . 

n=m^ n=m^ 

This goes to zero since the geometric series at ^ converges. 
The E and V cases are virtually identical - simply do the same thing with (2) 
and n rather than (g) . So once again up to the measurability we are done. 

q.e.d 

Let be the set of ordered triples in x^M which are on circles of radius less 
than 5 , and let r be a measurable function on x^M with its restriction to Vs in L^iVs) 
and which is symmetric under any permutation of the coordinates. Let i? be a random 
variable on which is given as a -R(p) = ^^gXp ^(^) '^h^re by t we mean the ordered 
triple in p corresponding to the face t G K^. (that this function and the functions 



introduced below are measurable can be easily seen; see section 3.2.4 if there is any 



confusion). For such a random variable we have its expect value given by: 
Theorem 11 With the above notation we have 

l3 



E,(R) = ^ / r(y)e-'^(y)MA=^(y) 



Proof: Given a set s = {ii.,i2-,i-i : i\ < i2 < "is} C {!,... ,n} let vr^ be the 
projection mapping of Mi x • • • x M„ onto M,^ x Mi^ x Mjg. Now note that every triple 
will be uniquely represented as set ■ vrs(p) by one of these s. Define a function fs{x) 
to be one if p G x^Af„ and the triple set ■ vrs(p) is on a disk of radius less than 5 and 
has its uniquely associated open disk is empty of points in the configuration set(p); and 
zero otherwise. Let Rn{v) = Ss=i ^('''"s(p))/s(p) on x^M, and note the above random 
variable R is precisely the function defined on P = Unez+ which is Rn on each 

x"M_. 

This formula buys us a better look at Ea(R). First noting that x^M- and 
x"M differ by a measure zero set and breaking up the integral into the pieces over the 
disjoint pieces of the space we have 



E,(R) = V / R„(p)dA"^e-^^ 



n=3 ' 

Using our formula for i?„ and the linearity of the integral we have 
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oo (3) , ji 

E(I^) = Y.Y. r(vr,(p))f,(p)dA"-e-^\ 

Now by symmetry each r{TTs{v))fs{v) (foi' ^ fixed n) has the same integral so we may 
write E(F) as 



^"^V E 7 ^ / r(vr,(p))/.(p)dA«. 



To do this integral note ?'(7rs(p))/s(p) is a function which is zero off the set with 
two properties. First of all, the triple set{'Ks{v)) is on a ball of radius less than 5 - call 
it B[set{'Ks{p)))] and, secondly, all of the other points in set(p) are in B{set{'Ks{'p))Y ■ 
So this set can be given by Uj/gV* ^ '^-^(y)'^)) with y denoting the coordinates 
onto which tTs projects. Note that on this set our random variable is r{'Ks{x)) (at least 
off a measure zero subset where it is zero) and so is constant on each of the disjoint 
y X (x"^'^i3(y)'^) pieces of the above set. So Fubini's theorem (with the r G L^iVs) 
assumption) tells us 

/ r(^,(p))/,(p)d^" = / r{y) ( [ l,.-^BiyrdA^-') dA^{y) 

= / r(y)(A-a(y))"-W(y), 
JVs 

where a{y) is the area of the ball function for y with set{y) containing three distinct 
points and zero other wise. 

Plugging this into the above computation, and using Fubini's theorem once 
again, we have 



A3 



E(F) = e-^^^ / - a(y))"-MA=^(y) 

6 ^3(n-3)! Jv, 



e 



^ r{y)e^^-^^y^^^dA\y) = — I r{y)e-<y^^dA^{y). 



6 Jvs 6 Jy^ 



Q.E.D. 

Note in particular when r = 1 we have equation (^]5) from the introduction. 
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3.2.2 The Geometric Coordinates 

Now its time to carefully set up the coordinates form the introduction, and to 
derive equation (|1.2| ) from the introduction. Here we will parameterize a full measure 
subset {Vs)- C Vs via the set 

Ls = x^^l X (0,6) X {M - {points}}. 

This mapping will be based on the exponential mapping expp : TpM M (see [^]), 
and its definition requires a orthonormal frame / = {61,62} - which can be defined at 
all but a finite number of points {{points} f) of M. Identifying x^S"^ with the triple of 
variables (^i mod 27r,02 mod 2ti,6^ mod 27r), denoted as 6, we can explicitly define 
our mapping <I>j : — > (14)- as 

^f{e,r,z) = {expz{rv{ei)),expz{rv{92)),expz{rv{e-i))), 

where the points excluded from M are precisely {points} f . 

Clearly we land inside the triples on circles of radius less than 5 and so are well 
defined. The differentiability of the exponential map guarantees <^>/ is a differentiable 
on Li. We would of course like to say more than that we have a differentiable map. 
As we know by lemma Q a triple uniquely determines its disk when the radius is less 
than 5, so the z and r coordinates are uniquely determined. Also set^^{set ■ ^ f{9,r, z)) 
corresponding to the 6 permutation of a distinct triple is hit precisely by the 6 distinct 
images under $j of {set~^{9i,9i,6i},r, z) with r and z fixed. So is injective. Also 
$j's image is almost everything since any triple in (Vs)- not centered at z G {points}/ 
is hit; a set which his of course measure zero (see section |3.2.4| ). 

Being bijective (to a full measure set) and differentiable, $ / forms a re-parameterization 
of a full measure subset of Vs; and with it we may continue the computation from the 
previous section this section finding 
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The remainder of this sub-section will be dedicated to finding an explicit formula 
for <5j((iA^). Since we are pulling back a Riemannian volume form, understanding how 
the metric pulls back will determine how the the volume pulls back; and this will be our 
strategy. Pulling back the metric in an arbitrary frame is fortunately not needed; in fact 
the expression for the volume form at {9,p, r) when pulled back by <5j depends only on 
the orthonornial basis choice at p. Letting dA'^ be the the volume form on Vs (since Vg 
is an open set in x^M) we can express this as: 

Lemma 27 Suppose the frames f and g agree at p E M , then 

^) [dA") {e,r,p) = % [dA^) 0,r,p). 

Proof: First notice that given two smooth orthonormal frames / = {61,62} and g = 
{di, ^2} defined on a simply connected open set of M (call it U), they can be orthogonal 
compared via a smoothly varying element of 0(1). U being simply connected allows us 
to lift this mapping from 0(1) to its cover R-*-. Otherwise said: there is a differentiable 
function Of,g{z) such that 



ei 




COS{e fg) - Sm{e fg) 




di 


62 




Sm{e fg) COS{e fg) 




d2 



Let F be the mapping 

F : y.^S\ X (0,(5) x {[/ - {points]] x^Sl x (0,5) x {U - {points]] 
such that 

F : {e, r, z) = (^1 + Ofg, 02 + 9fg, 63 + 9fg, r, z). 

Using this mapping note that our mapping satisfies = $g • F on x^Sl_ x (0, 5)x {U — 
{points]]. 

Now simply note that using d{6fg{z)) AdA = we have 
F*{v{9, r, z)dA AdrA d9i A d92 A d^g) 



= v{9i + 9fg, 92 + 9fg, 93 + 9fg, v, z)dA A dr 
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Ad{0i + Ofgiz)) A d{e2 + efg{z)) A d(03 + f{z){z)) = 

viPi + Ofg, 02 + 6fg, 6*3 + 9fg, T, z)dA hdr h d9i A de2 A dO^. 

Assuming the frames agree at p, we have dfg{p) = 0. So by the above formula 
F* acts as the identity at this point - giving 

{dV) ie,r,p) = F*^; {dV) i0,r,p) = ^*g {dV) ie,r,p), 

as needed. 

q.e.d 

To actually compute the inner products it is useful to canonically identify all 
the tangent spaces near p as as and r vary with a frame determined by the triple 
product of the normal coordinates at p, N{p,6,r), via 

Lp{e,r) = x^iV, : x^E^ ^ T^^(,-r,p)V5. 

While not an orthonormal frame a certain amount of the Euclidean geometry in 
x^E^ is preserved by this frame. For starters each of the copies in x^E^ is orthogonal 
to the others, since is an open subset of x^M in its product metric, and this mapping 
is respecting the product structure. Denote the vectors in this x^E^ as {wi,W2,W3). 
Further more by Gauss's lemma each of these Wi components decomposes orthogonally 
(when < r ^ i) into v{6i){p) and v-^{6i){p) - with respect to the both the Euclidean 
metric and the surface metric. So for each triple of angles we may in fact represent the 
vectors with the following orthogonal decomposition 

{aiv{ei){p) + birv^{ei){p),a2v{e2)ip) +b2rv^i92)ip), 

a3v{93){p)+hrv^{93){p)). 

To understand the lengths of these vectors in this decomposition it is necessary to remind 
our selves about Jacobi Fields. 

We can now get a grip on our needed vector lengths. 

Lemma 28 Suppose f is a geodesic frame at p then using the above notation we have: 
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Similarly for 62 and 9^ . 



■ (^/), (ei) = (cos(0i)t;(0i), 005(^2)^(^2), cos(e3)t'(e3))+ 

(-^ sm{e,)vie,)^, -!:!^smi92)vie2)^, sin(03)i^(e3)^). 
Jeo Jeo Jeo 

Similarly for 62 . 

Proof : The idea for all these computations is identical. For each component one finds 
a curve T{s) such that j^^{s) = vrj-v and then one notes that in fact T{s) = T(r,s) for 
some two parameter family of geodesies - forcing our vector, v , to be Jacobi fields. At 
this point one uses the Jacobi lemma (lemma 13 in section |3.1.1| ) to observe the above 
formulas. 

For example for (<I>j)^^ is by our Jacobi field observations is 
-^{expp{rv{s)),expp{rv{92)),expp{rv{e3)))\0^. 



By the discussion preceding lemma |13| this is the Jacobi field with initial conditions 
J(0) =0 and ^(0) = v-^{9i) - located in the in the first component of M^. So we have 
a description of it from part 3 of lemma ^ as {rv{9i)-^{p),0, 0). 

Now let's play the same game for the radial direction and note that the image 

of (^/)*f is 

-^{expp{{r + s)v{9i), expp{{r + s)v{92), expp{{r + s)v{93))\s=o, 



which by part 1 of lemma 13 is J(r) = {v{9i){p), {v{92){p), {v{9^){p)). 



86 



Finally to compute the vector in the ei direction use the geodesic pointing in 
that direction and vary along it. This is the situation of the example preceding lemma 



13, and we have (^/)^ei described as 
d 



ds 



{expexppise,) (rviei)), expexppise,) (?'^^ (^'2) ) , exp^xpp (se,) (rv (6*3) )) . 



Since the frame we are using is geodesic at p we have —^-^ = so -^(0) = (by 
the example preceding lemma ^ - once again) . By the discussion preceding the Jacobi 
lemma p^, we may decompose this field by decomposing J(0) = ei{p) into v{6i) and 
v{6i)^ in each component; so we have this Jacobi field has one summand corresponding 
to 

j(o) = {<v{ei),ei>v{ei) 



, < v{e2),ei > v{e2), < viOs), ei > u(^3))(p) 

and -§^(0) = 0. This field is dealt with by the first part of lemma |l^, and gives the first 
set of vectors. The other summand in the decomposition corresponds to 



j(o) = (< viei)^,ei>viei)- 



,<v{e2)^,ei>v{e2)^,<v{e^)^,ei>v{e^)^){p) 

with -§^-(0) = - which is dealt with in by the second part of lemma These formulas 
are precisely the last of the needed inner-product relationships. 

q.e.d 

With the use of these lemmas we accomplish our goal of computing the pullback 



of the volume form in V5 coordinates and derive equation (1.2). 
Proposition 4 

^){dA^) = y{e)jeMesde AdrAdA 
with 1^(6*1) the area of a triangle with vertices at {v{9i)} in the Euclidean plane. 
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Proof : Let's compute the form at a point {6,r,p). Using the normal coordinates for 
the M variables note we have ^ = e^. We need to find {det{gij)'^ dzi A dz2 A d9 A dr. 



To do so first note by lemma 27 that we may use the geodesic frame at p with the initial 
vectors {ei{p), e2{p)}- Recalling that < v{6i),v{6i) >= 1, and that from lemma [I^ in 
section 1.1 < rv-^{6i),rv-^{6i) >= Using the above observation we now have can 
easily compute the determinant. 
Letting 

iy{ei) = I sm{e2 - Oi) + sin(03 - O2) + sin(0i - ^3)! 



■ sm 



sm 



we find 



Take its square root to find 



sm 



det{gij) = v\e)jljl^jl^. 



^*{dA^) = v{e)jejeje-ide Adr A dzi A dz2. 

Now dzi A dz2 is precisely dA at in these coordinates (at p) - so we have our 
needed equality. 

q.e.d 

Actually it is useful to witness the matrix forming the determinant in the above 
proposition and hence find that it is relatively easy to compute. Note 



< 



> 



^ JL 3 - 
^ ar • asi ■ 

^ _§_ a 



^ 3 3 ^ 



< 

> < 



a o ^ 
ar • 383 ^ 
a a 

3^1 ' 363 

3 3 

3^2 ' 3^3 



> < 

> < 



<' 3 

















.2 




























3 


-iOi ftfli sin(ei) 




-iflj^ej sin(e3) 




Jeihei cos(Si) 


J»2 ™=(^2) 


302^^2 ™''(*2) 





-je^he^ ain(e3) 

Efe = l(''9(,i^in2(efc) + cos2(efc)) 

El = l(sin{efc) cos(Sfc) - h2 sin(Sfc) cos{et.)) 



J9i 'i^i coa{ei) 
3^2^62 cos{e2) 

oe^hg^ cos{e3) 

ELiCja. cos2(Sfc) +sin2(Sfc)) 
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so one can pull the jg^ out of the matrix and the symmetries of the remaining matrix 
make the computation of the determinant relatively transparent. (Trying to see this first 
hand proves useful when one begins exploring the three dimensional computation.) 

3.2.3 The Euler-Gauss-Bonnet-Delaunay Formula 

It is now time to compute the needed random variables. All the random vari- 
ables encountered have expected values expressible as 

Ea(R) = — / r=^f (r, p, ^)e-^('''P'^)Mrd^dA 

where 

fir,P, 0) = {Hp, 0) + h{p, 9)r+ Mp, ey ln(r) + Mp, 9^ + 0(ln(r)r3)) 

with the fi bounded and measurable functions on L^. In this expression and throughout 
out this paper remainder a function denoted 0{f{r)) will mean a function 0(p,6,r) on 
M X [0, 27r]'^ X (0,(5) satisfying \0(p,6,r)\ < Cf{r) for some constant C. 

To compute the above expected value first note that one can Taylor expand in 
r the area of a ball of radius r at p as a(r, p) = vrr^ — ^^^^ + O(r^) (in fact this follows 
immediately form the o.d.e. describing the v^{6) Jacobi fields and the fact the surface is 
compact). This formula allows us to choose constants < e < 1 and pM > such that 

g-Aa(p,e,r) ^ g-A7rr2(l-e) 

for all pm > r > and at all p G M. 

Theorem 12 Using the above notation, letting c = \n{r)r^e~'^^ dr , and using a de- 
cision radius smaller than pM we have E;^(R) equals 

6 7m 751x51x51 27r2 6 JmJs^xS^xS^ IStt^ 

f f / 15/i(p,g) ^ cf2{p,e) ^ 105k{p)fi{p,9) 

6 Jm 751x51x51 V 167r2 ^1 12 -32713 

1 / / ^ tmM] S4A + 0(ln(A)A-4). 

6 7m 751x51x51 V TT'^ 47r'^ J 
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Proof: As a preliminary observation, the above fact about a{r,p) and the mean value 
theorem gives us 

Using this fact E;^(R) equals 

^ / {Mp, e) + hip, 9)r + f2{p, ey ln(r) + h{p, 0)r^) e-^'^^drSdA, (3.1) 

+ Y / (/o(p, 0) + e)r) ^^^M^le-^r^^drSdA, (3.2) 

+ Y / {f2{p, ey ln(r) + h{p, ^^^^^^e-^^'^drdOdA, (3.3) 

+— j r^^^^^^^^^f{r,p,e)0{Xr^e-'''-^^)drd9dA, (3.4) 

+^ / r3^^^^^/(r,p,^)O(AV«e-'-'^(^-^))drd0"'dA (3.5) 

6 / r 12 



PM 

+4 I r3 ^''y''' Q(r^ln(r))Q(AV^e--'^(^--))drdg"'d^ (3.6) 

6 ./r_ 12 



6 A 12 

^PM 

The first thing to observe at this point is that for any 5 

f e-^^'\''dr = ""t+i + 0(A-°°) 
•''0 (7rA)~ 

and similarly with the ln(r)r'^ integral. This is because 

1 



[ e >^^'^^r'^dr= f e \^s''ds 

Jo Jo (7rA)~ 



poo 2 roo 

= / e-* ^s'^ds - / 1 - ,+1 



e r-rS^ds 



'0 (7rA)~ ^5(7rA)2 (;rA)" 



"4Tr + ^Mrrf'(<^(-A)^) 



(ttA)— (ttA)" 
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Where d^{x) = e~* s^ds - which by Phoptial's rule decays faster than any 
polynomial. 

Armed with this observation wc may explicitly integrate the first two terms in 
the above function with the only expense being a 0(A^°°) term. After doing so we find 
the non O ^ln(A)A~2 j terms in E;^(R) to equal equal 

U f (/o(f.»)^ + /*.«)^+ 

o JmJs^xS^xS^ \ 7r2 

hip, 0)1 ^ + H^X)X^mA ^ ^^^^ ^^rnA ^.^^ 
\ 7r2 7r2 y TT y 

i> Jm Js^xs^xs^ 1^ \^ TT'^ 7r2 y 

where nik = J^r'^e~^^dr. Plugging in mz = \,m2 = = 1, my = 3 and 

mg = ^°3^ we arrive at the above formula. 

Likewise upon integration we see the remaining four terms are indeed 0(ln(A)A~2 ). 

Q.E.D 

We now can compute this for the random variable F counting the number of 

faces. 

Corollary 2 With 5 < pM o-nd calling J^s^i i^{9)d6 = u, we have the expected number 
of faces in the configuration can he expressed as 

Proof : From the theorem in the previous section the f{r,p,6) corresponding to the F 
random variable is 

Using the o.d.c. describing the Jacobi field and the fact M is compact we now im- 
mediately have je^s Taylor expansion is jegir) = r(l — + 0(r^)), and as such 
jeijeijda = ^^(^ - l?^^ + o{r^)). So we have 

f{r,p,e) = i.ie)(l-^ + 0{r')\, 
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as needed. Q.E.D. 

Let us now use the Euler-Delaunay-Poisson formula (formula [l|) to figure out 
the constant ( if you prefer, the integral it represents is easy to compute). We know 
the Euler characteristic of the flat tours of area one is zero; so as A ^ oo we have 

= E(x(M)) = A - ^E(F) = (1 - ^)A, 

forcing v = 247r^. 

Plugging this into the formula in the corollary we have the needed equation 

(0). 

3.2.4 Silly Continuity Results 

We have several functions which we need continuous, and several sets we need 
either open or closed and measure zero; here we deal with these issues. Let us start by 
collecting and dealing with the needed closed measure zero results. 

Fact 4 The following sets are closed and measure zero where indicated: 

1. {(xi, • • • , Xn) e x"M : Xi = Xjwithi / j} in x"M. 

2. ^f{LsY in Vs. 

3. The set of points x in x^M where four points in set{x) lie upon a circle of radius 
r < mm{|, r}. 

4- The set of points x in x'"'M where three points in set{x) lie upon a circle of radius 
exactly 6. 

Proof : We will simply realize these sets as finite unions of the differentiable images of 
lower dimension sets, hence measure zero. They will all be clearly relatively closed in 
the specified sets. 

Let D : M ^ M X M be D{x) = The set we are removing is precisely 

the union of the nl index permutations of the set D{M) x (x"~^M), as needed. 

For the next part I'll prove something a bit stronger. Let V^^ri. be the set 
of all triples on circles of radius less than or equal to mm{|,r} with repeats allowed. I 
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will show the stronger result that $j acting on x^S^ x [O, min{|, r}] x M — {points} j 
hits all but a measure zero set of V^j„|i ^y- To see this let 5 be a frame defined at each 
point of {points} f, and then we have missed exactly 



0, mm{-,T} 



X M — {points} g){z = p) 



(J i^g{x^S^ X 

p&{points) f 

This is a finite union of the differentiable images of 4 dimensional sets so measure zero 
and clearly relatively closed in the 6 dimensional space. 

For the third part first look at the image of the the 7 dimension space 

T" 



X (M - {ponts}f) X 



x^S^ X 

under the mapping given by 

{^f{e,r,z),expz{rv{ei)). 

It is measure zero as a subset of x^M. Union it with the same set using the above g. 
With this unioning it is clearly relatively closed in x^M. Now product this with the 
other n — 4 components and take the union over the n\ permuted copies of this set in 
x"M and remove them; hence removing all the possible quadruples on the same circles. 
This amounts to a finite removal of closed measure zero sets, as needed. 

Similarly for the last part, where we use instead the the image under (and 
^g) of x^S^ X (5 X M in x^M. 

q.e.d 

Now we need some functions to be continuous and some sets to be open. 
Fact 5 is open in P; and F , fg, E and V are continuous in P. 

Proof : All these cases could we handled simultaneously if one could perturb 
a set of points and not change any of the assignments given to a triple, a pair, or a 
singleton. We may restrict or attention to a set coming form x"M_, since P is a disjoint 
union of such sets. Clearly the assignment is discontinuous on x"M; and what is being 
claimed is that all the discontinuities occur on the removed measure zero sets. 

That any singleton is still a distinct vertex is true due to the removal of the set 
corresponding to part one of the above fact. In fact this removal guarantees all points 
in set{x) can simultaneously be separated by open sets (including any pair or triple). 
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To deal with the stability of face assignment to a triple, set-iTs(x), first observe 
the removal of the fourth set guarantees Vs and are open. So set ■ 7rs(x)'s relationship 
to 6 is open, finishing off the stability argument if the point fails to be on a circle of 
radius less than 6. In the case where the triple does lie on such a circle, note that 
TTsix) G ^f{Ls) for a suitable /. So the disk associated to a triple varies continuously 
with the triples position (r and z in the parameterization). Now suppose set ■ tTs{x) 
is not assigned a face due to a point in its associated disk's interior. This violating 
point must have a neighborhood in the disks interior. So by the above continuity of disk 
position, there is an open set Ui x U2 x U3 about TTgix) and an open set Up about p 
such that Up is contained in the disks associated to set(y) for each y £ Ui x U2 x U^- 
So the triple having no associated face is indeed an open condition. Virtually this same 
argument guarantees that in the case when the disk is empty that there is an open set 
about set{x) — set ■ iTs{x) such that no point in this open set is in the disk associated 
to set{y) for y in some Ui x U2 x U3. The only difference is that we must note that we 
removed the possibility of four points on a circle (part 3 above) - so indeed each point 
in set{x) — set ■ iVsix) is in an open set separating it form the closed disk, as needed. So 
off the above measure zero sets the notion of face is indeed stable in an open set. 

The edges require a small amount more thought. An edge pq existing implies 
there is a A; such that Dpq{k) is empty. We may attempt to deform the inflating family 
to the left and to the right of A;. If we cannot it is because a third point lies to the left 
or right side of the circle. Now since we may assume there is not a fourth point on the 
same circle, in one of these directions so we indeed can deform our inflating family; and 
in fact in both directions if the circle's radius is 6 - since no triple lives on such a circle. 
In particular, choosing a different k if necessary, we now have a circle of radius less than 
S which is empty, and is contained in an neighborhood empty of other vertices. Now one 
can proceed exactly as above to note that E is constant on a neighborhood. The notion 
of no edge requires a more delicate use of our inflating family ideas. The key observation 
is that from the monotonicity lemma a point is in intDpqQi) P| intDpq{k) if and only if it 
is in each intDpq{c) for c € [h, k]. Prom this the idea is to find a = cq < ci < ■ ■ • < c„ = 6 
such that Dpq{ci — e) Dpq{ci+i + e) contains a point in its interior. Then each of these 
sets will, as above, satisfy this property when the points are perturbed - and the notion 
of no edge will be stable. (Note at a and h we need not use the e - since there is a point 
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in the interior of each of these disks which we may use.) To construct the q start at 
o and take a point inside it. Move the inflating family rightward until at some di this 
point fails to be in the disk (if such a point does not exist we may use cq = a and ci = b.) 
Now since there is no empty disk intDpq{di) contains some other point. Being in its 
interior this is in fact true for the parameters in [di — ij, di). Let q = — ^ and e = f • 
Now continue this process making e smaller if necessary. There are a finite number of 
distributed points - so eventually one must hit b or an empty disk. An empty disk is 
impossible, since no edge was put in; and we are done. 

q.e.d 



Chapter 4 

The Continuous Uniformization 
Theorem 



This section accomplishes two things. For starters in section 4.1 we use theorem 



12 of the previous chapter to actually calculate the energies on the space of metrics of 
interest form section 1.2.2 and 1.2.3. This involves first a careful look at how the angles 
in a geodesic net deform under a conformal breeze followed by the actual computations. 
Section 4.2 contains an actual proof of the uniformization theorem for surfaces with 
x(Af) < 0, mimicking in a precise sense the discrete uniformization proof. 

4.1 The Random Energies 
4.1.1 Triangulation Deformations 



The key to computing the need expected values in section 1.2 is to compute how 
an angle in a triangle deforms under a conformal change of metric. The method used 
here is to solve the the boundary value problem for the geodesies in the e'^'^g metric, and 
compare the initial directions. Of real interest is the case when we have the geodesies 
forming the edge of a triangle, so I will phrase the results in this language. Using normal 
coordinates at the point p in the g metric let the points on our triangle are labeled 
rv{9i) = rv^ = (r cos(0j), r sin(0j)). Further let Vij = v{6j) — v{9i) for any i and j and let 
vjj be the left handed ^ rotation of Vij in the Euclidean metric. Given that the surface 
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is orientable pick an orientation and let a(i,j,k) be one if {vi,Vj,Vk) is ordered in the 
"clockwise" direction and —1 if not. Further more let k be the Gaussian curvature at p 
and H be Hessian of (p in normal coordinates (which are denoted via {x,y)). 

Lemma 29 In the above notation the initial direction of a geodesic from to in the 
e^^g metric is given by 



(^(%-V^) + 0(r2)) 



k 



Proof: For the purposes of readable notation let Vij = v = {v,w), and note we 
may rewrite our initial starting direction as v + 1. 

The new metric in normal coordinate is up to order O(r^) given by 



1 — ^ky"^ ^kxy 



hkxy 



1 - ifcx^ 



Recall the geodesic equation is 



Using this a metric up to order 0{r) the Christoffel symbols are found to be 



and 



-1- rr 



(t>x + <l>xxX + (t>xyy + <i>yyy + <l>xyX - \ky 

(t>y + (l^yyU + <t>xyX -Iky -(^x + (l)xxX + (t>xyy) + ^kx 



-{4>y + (l>yyy + (t>xyx) + ^ky (j)^ + 4>xxX + (l)xyy - \kx 

(t^X + (t>XxX + (t>xyy - IkX (f)y + (f>yyy + (f)xyX 
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All the above <p derivatives are evaluated at the center of the circle, p; and I am 
letting k denote the curvature a p. We are solving the boundary value problem where 
wc start at rv^ and ending a rv-' , i.e. introducing the notation that {p,q) = {x,y) — v 
we have 

p + v 

q + w 

{T\^{p + vf + 2T\^(jp + v){q + w)+Tl^{q + wf) 
(r}i (jp + vf + 2T\^ (jp + v){q + w)+ Tl^ (q + wf) 

Since we know the ChristofFel symbols up to 0{r), to analyze the solution up 
to order we may simply linearize the problem form the point of view of the origin 
(this can be seen immediately by looking at the power series expansion). So we need to 
solve the linear o.d.e. 



X 

y 




p 




_ q _ 





X 



X 

y 
p 

q 



h + Ax = h + 



/ 

-C -2D 



X, 



where b 



V 

f 



V 

w 



and 



C 



-{(Pxiv^ - W^) + 2(f>yVW) 



.{v^ — iiP') + 2(f)xyVW + ^kw'^ 



v^(v^ • V(/)) - v(v V(/)) 



f^xy{v^ — w'^) + 2(f)yyvw — ^kvw 



— (f)xy{v^ — 'Up') + 2(j)xxVW — 'ikvW —^yy{v^ — w'^) + 2(t)xyVW + Ifcv 



D 



bxV + ^yW —(f)xW + (pyV 
bxV + (l)yW 



v(V0)*'- + v^(V</.^) 



tr 
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For future use it is convenient to note 

Df = V ((V<^ • v^)2 - (V.^ • v)2) + 2v^ ((V.^ • v^)(V<^ • v)) 
Now this o.d.e can be solved via variation of parameter via 



X = e"" I I e ^^hdt ] + e*"xo 



with initial condition 



xo 



rv 
1 



r cos(0i) 
r sin(6'i) 
I 

k 



So up to order o(r^) we have 



^ = {t+^-A + jA'^)h + {I + tA + ^^2)xo. 

Now we need to find the I and k so that when this is evaluated a r it is at the 
point rv^ . 

So plugging in r we find we need 

(7 - rD)l = ^ ((Cr«* + ^Cv) - - ^D)i 

Inverting to isolate 1 we find (up to order 0(r^)) that 

1 = ^ {rC(v' + - (, + 5z))f) = + ^ (C(v' + iv) - i«f) . 

Simply plug into this expression the above formulas to arrive at the claimed 

formula. 

Q.E.D. 

Now we are capable of computing the needed angle. It is convenient to denote 
the angle in the triangle at the point rv* in the Euclidean coordinates as Ej. With this 
we have... 
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Lemma 30 The angle in the triangle at rv'^ is 



where 



= E, + a{ij, k) {-\{vkj^ ■ V<^) + + O(r^), 

~ (((V(/. • Vji^){V(t) ■ Vji)) - ((V(^ • Vfe/)(V(^ • Vki))) . 

Proof: For this problem denote the solution to the previous problem as 1 + v = 
(1 + mij)v + Uij-v-^. First note that the angle, r]ij, between this initial vector and v has 
its cosine given by 



cos{rjij)^ <(l + v),v><; 



|v|U)(||l + v| 



Hence using the notation of the previous problem and letting ci = TpTrf , cs 
± 2 II^IIe 

j!^ .ji^ , and C3 = ^rprp^, we have up to cubic order that 



cos{rjij) ^ ^ '—^ 



+ ci)(l + 2mij + mfj- + n?ci + 2mijCi + 2nijCz) 

Multiplying out and using + x = 1 + — + j^x^ + . ■ ■ and = 
1 — X + x'^ — x'^ . . . up to the third order in r we have the nice fact that this expression 
(up to r"^) is independent of the 50 rnetric and equal to... 



n2 



as... 



cos{r]ij) = 1 + ( — ^ + ruijUi^). 
Using cosine's power series this allows us to isolate the angle up to second order 

WA = WijW'^ - 2mij = \nij\(l - rriij), 
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Further note that if one would hke this angle to positive if it contributes to the triangle's 
internal angle and negative if not, then 

Vji = (^{hj, k)nij{l - ruij). 

With this computation out of the way we are left needing to measure the angle 
between Vij and Vik in the g metric. The trick will be to measure the angle from 
—v^ to Vij at rv^ and the angle a^i from — to v^i at rv^, both with the proper signs as 
contributers to a^, and then sum them up. Let afj be the angle in Euclidean coordinates 
and (•, •) be the Euclidean innerproduct. 

To compute these angle note at rv^ we have Wv^lg = 1, = 1 — + 

o(r3), \\v^ - v^Ie = \/2Vl - {v\v^). and that 



cos(aij) 



< V , V 



{\\V%){\\V'-V^\\g) 



1 - {v' 



^(1 - ivi,vj)y + i{vi)^,vjyil - V2^yi-{v\v^)^ 



V2 V2 V 12 

Now note that 



cos(a,^) 



E Vl-iv\V^) 



and from this by keeping track of the necessary sign we have 
sin(|aj^|) = a{i,j,k) 



So the above expression is precisely 



cos(a,^.) - sin(|a,^|) (=^^^{{vYy)) = cos (a^ - '-^^{{vY,vn) ■ 
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So we have that aij = afj — {{v^)-^,v^), and from this the needed 

a; = Ei + a(i,j,k)^^((vy,vJ)+a(i,k,j)^((V)^,v'^) 



Ei + a(i,j,k)^^(vi,v') 



Now summing up to get ■0^ = Vji + Vki + ci^^, and so 
= Ei + cj(i,j,k) — (vjk,v') +cr(i,j,k)nij(l - rriij) + cr(i, k, j)nik(l - rriik) 

= Ei + o-(i,j,k) f — (Vj4,v') + nij(l-mij)-nik(l-mik)j . 

Plugging into the formula from the previous lemma now finishes the computation. 
Q.E.D. 

4.1.2 The Energy Computation 

Using the angle formula from the previous section I will now compute the the 
expected value of the energy, deriving formula L5 and theorem |6[ In the end we will 
arrive at the formula 

Formula 13 For a negative curvature metric g and h = e^'^g we have 

Ef (Eh) = DoAA + Dix(M) ln(A)Ai + D2x(M)A^ + D3x(M) 

+ / ||V(/)|p + (A(/.-A:)ln(A(/)-A;)(i^ + o(ln(A)A-5), 
with the Di constants. 

Proof: To do the computation it is necessary to find the f{r,p,0) function 



used in theorem 12 with respect to the volume induced random variable. 
By changing to the (j) coordinates we reduce E^ (Eh) to 



6 



/ / r yip. r, e)e-^''\{ey{l - r^^ + o{r^))drd9dA, 
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where V{p, v, 9) is the prism volume associated to the triangle withe angle data given by 
the triangle on this surface formed with this data. So we need to expand 



^ ^ k 
V{p, r, e)v{9)r'^{l - r^- + o{r^)) 

in the r variable. To do this it is nice to give our smah objects some names; let A 
TT — V'rf,! Ist 5i = ipi — Ei,and let = + 4 ~ -^i- '^^^^ which is small since 



TT + Ei 

Note that at this point 



2 

TT + Ei - Ej - Ek 



V{r,p, e)=Y, A(Ei + Si) + A(Ei + ^i) + A 
The power series expansion of minus the Lobacevskii function A about a positive 



E is 



A{E + 5) = A{E) - ln(2| sin(^)|),5 - cot{E)5^ + O ( r-^.^o , ■ 

\ (sm(£/jj^ / 

Note from the formula for tpi in lemma ^ that the 6 and 6i functions are divisible 
by sin(£'j) = sin( ^^ ^^'^ ) so this series when applied to the first six terms in the above 
expansion for V{r,p,9) gives us a power series in r with bounded continuous coefficients 
and a remainder of order O(r^). 

Note that the remaining term is in the form 

= ^-^ln(^) + 0(^^). 
Given this expression it is extremely useful to get a grip on the expression 

Lemma 31 



A = ^^Ptf,-7T = (-A0(p) + kyu{9) + 0{ 
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Proof: 

Let i < j represent a pair where cr{i,j, k) = 1. 

We can begin with the observations that X]i<j ^ij ~ that since both 

i < j and j > i include all ordered pairs 

i<j i>j 

^((V</. • Vji^)(V</. • Vji)) - ^((V</. • Vki^)(V0 • Vki)) = 

i<j i>j 

and 

E - E = E i^^THvf - vfHvt-) 

i<j i>j i<j 

i i<j 

To compute explicitly the remaining terms the following fact useful: 
Fact 6 If A is a symmetric 2x2 matrix then 

J^iv^y^Av' = - = tr{Ayiy{9). 

Using these observations, the above fact, the fact that the angle in the Euclidean 
triangle sum to zero, and the fact that tr(H) = A^, we find the sum is 

(-2^+4)rM^-*) + 0(r3), 

as needed. 

Q.E.D. (lemma) 

It is worth noting as a confirmation to the previous sections computation, this 
is lemma is exactly what one expects from the the Gauss Bonnet formula when = 0. 

With this in mind we see that in fact the expansion for — A (^) has terms of 
the form /r' for i > 2 and g\n.[r)r'^ for z > 2 with continuous coefficient functions. So 
up to order ln(r)r^ V{r,p, 6) is of the form 
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y{r,P, 0) = Y, (A(Ei) - ln(2| sin(Ei)|)(5i + k) - cot(Ei)(5i + k)'' 



i=l 



+ ^ - ^ + 0(ln(r)r3) 



With this we are in position to apply theorem 12. Note that 6i + Si = 2il}\ + 



^ — 2Ei. So using the notation of the previous section we have the fi{p, 9) from theorem 
and we get 



UpJ) = Y^ki^-MO) (4.1) 

i=l 

fiiP,0) = - j;in(2|sin(Ei)|)(v(0k)-v(0j))-VMe) (4.2) 



f2{p,9) = (AcP -k){u{9)y (4.3) 

^ k ^ /\(h — k 

h{p,9) = -A(Ei)^.(0) + ^— ln(2|sin(Ei)|)(z.(0))2 (4.4) 

-2&iln(2|sin(Ei)|)i/(^) (4.5) 

-\{^fk-y<t^fcom)m (4.6) 

+^(A0- A;)ln(A(^-A;)(z^(^))2 (4.7) 

Now we attempt to compute Ii = /;y^/[0 27r]3 fiiP^ ^)d9dA. The first one clearly 
can be integrated to a constant times the surfaces area, Iq = CqA, where the constant 
is independent of any of the geometry or topology. 

Happily enough Ii = 0. This follows immediately form the first part of the 
following integral vanishing lemma ( after noting sin(Ei) = sin( )): 

Lemma 32 Let f and g be either cos or sin functions then 

J^^ ^ In (21 sin (^^^) 1) f{9kH9)d9 = 

^^In (^2|sin (^^Y^) 1) f{9k)fi9M9)d9 = 
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J In ^2| sin (^^^y^) sin(6'fe) cos{ek)u{9)de = 



[0,27r; 



[0,27r 



3 



In (^2| sin (^^y^) l) + fi9k)9{e,)Me)d0 = 



In ( 2| sin ( ^ ) I ) i{f{ek)9ie,) + = 



Proof: The idea is simply to note in each case that there are symmetric regions 
of the 9 cube where the function has opposite signs (the finiteness once again fohows form 
the integrabihty of the logarithmatic singularity). I'll simply indicate the symmetries. 

One key observation is that u and In ^2| sin ^^^-2"^^ |^ ai^e invariant the trans- 
formation, Tc, were for each i 6i goes to + c modulo 2tt; and the transformation N 
sending all 9i to —9i modulo 27r negates v while of course leaving In ^2| sin ^^^y^^ |^ 
invariant. 

Now note for the first integral either N of N ■ T^^ will produce the needed 
symmetry. For the second integral N alone will work. The Remain integrals have the 
order four transformation producing the four points - two of each possible sign. 

Q.E.D. (lemma) 



In I2 note that after integrating the Laplacian term in equation 4.S integrates 
away and the k term integrates to an Euler characteristic - so in the end we get a constant 
times the Euler characteristic I2 = C2x{^)i with C2 depending on none of the geometry 
of g or topology of M . 

can be broken up into the four pieces as in its formula (equation [4.4| - [4.7] ). 
The first piece (equation |4.4| ) can as above be integrated out to give C^x{^)^ with C3 
depending only on the Euler characteristic. 



The last two equations (l.C and 4.7) can be explicitly integrated to give 



-evr^ j I |2 + (A(/. - k) ln(A(/) - k)dA. 
Jm 



This leaves equation 4.5 involving the 6j. Looking at the expression for bi we 
see that the integral vanishing lemma immediately gives us that these terms integrate 
away to zero, with exception of the terms in the form 



106 



i^k + (P,,) j^^ ^ In (^2| sin (^^y^) l) ((" ^HOk) + sin(%)) zos{di)u{e)d9 

{\k + <^yy) j^^ ^ In (^2| sin (^^y^) l) ((cos(efc) - cos(0j)) s\n{eiy0)d9 = 0. 

However the last integral in the vanishing lemma tells us that these integral 
are equal. So as in the previous piece we are left with a term which integrates out to 
C3X(M). 

So using theorem |l2| our energy is now indeed in the claimed form. 
Q.E.D. (formula |l|) 

Now we can compute the energy in formula simply by noting that the terms 
in the formula with Di constants in them cancel (in the prelimit even) and that we are 
left precisely with the needed terms. Also noter the formula for Ei is immediate. With 
these we can verify the formula for E2 and prove theorem ^. 

Proof (theorem ^) It is immediate that amongst metrics of the same area, 

A, that 

E2 = Ei-E = - I kh ln{-kh)dAh + [ \ \V4>\ 1^ + (A0 - k) ln{A(j) - k)dA 
Jm Jm 

= - [ {A4>-k) In (e~2'^(A0 -k))dA+ [ 1 1 V0| ^ + (A<^ - k) \n{A(f) - k)dA 

= [ 2(f){A(P-k)dA+ [ \\V^\\'^dA = - [ \\V(p\\'^ + 2k(l)dA. 
Jm Jm Jm 

We now state the beautiful formula do to Polyakov (see [13| or [^) for the the 
determinant of the Laplacian 

ln(det(A;,)) = -^(l [ \V^\^dA + [ k(f>dA] + In(^) + C. 
6vr V2 Jm Jm ) 

Plugging in above we arrive at the needed 

E<yK) = 127rln(det(A,,)) - 127rln(A) + C. 

which is the claimed formula. 
Q.E.D 
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4.2 A Proof of the Uniformization Theorem 

The goal here is to give the proof of the metric uniformization theorem indicated 
in section |1.2|. It is a direct analog in infinite dimensions of the proof of the angle system 



uniformization presented in section 2.1.4. To see this it is useful to organize the finite 



dimensional proof into 4 basic steps, which will be mimicked amongst metrics. 

4.2.1 The Discrete Uniformization Proof Reviewed 

The negative curvature Delaunay angle system uniformization proof: 

Step 0. We defined what what we called the negative curvature Delaunay 

angle systems conformal to a fixed angle system y and called it Ny, which turns out to 

be nice convex set. Then we placed upon it an energy E 

teP 

which was continuous on Nj,. 

Step 1. Now we observed that at least one point of maximum energy must 
exist, here by the rather trivial observation that space Nj^ upon which our continuous 
energy lives is in fact compact. 

Step 2. At this step we showed that any x where E assumes its maximal value 
is in fact in Ny. To accomplish this recall we took a point x on the boundary of Ny 
and constructed a line l{s) satisfying /(O) = x and that there is an e > such that for 
all s G (0, e) we have that l{s) G Ny and E{l{s)) is increasing in s. Clearly now the 
continuity of E{l{s)) on [0, s) makes it impossible for E to assume its maximum value 
at any boundary point. 

Step 3. At this step we verified that any point x G Ny where E is maximal 
is in fact a uniform angle system. This followed immediately by examining what the 
differential of E vanishing at x implied about x (see observation Q). 

Step 4. We proved uniqueness of the point x where E achieves its maximum. 
To do this we took a second point y and connected it with a line l{s) to our x and note 
that 

E , , . 
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for all s such that l{s) remains in N. So y could not also be maximal (see lemma 

4.2.2 The Indiscreet Proof 

The proof in the previous section can be carried out in the indiscreet world step 

by step. 

Step 0. We need the correct analogs of an N and N on which to interpret the 

energy 

E{(l)) = - [ \ I V(/)| p + {A(f> - k) log(A0 - k)dA. 
Jm 

Namely which (p do we use to conformally change our initial metric. Clearly a rescaling 
will change nothing both since it effects no angles in the discrete model and since E is 
clearly scale invariant on C°° . In fact with this in mind we should feel free to re-scale the 
initial metric and for convenience let's assume it area is — 27rx(M). One way to eliminate 
the possibility of rescaling is to demand 

/ (t)dA = d. 
Jm 

We'd also like to restrict to (p with negative curvature, i.e. k^j, = e~^'^(— A0 + k) < 0. So 
the first guess at a reasonable function space might be 

V = {(l)e C°°{M) I [ (fxlA = and A(j) - k > 0}. 

However the energy on here would make using a compactness argument in step 1 rather 
difficult. The trick to producing a place where a compactness argument will work is 
to take the closure of V with respect to a norm which E interacts with in a sensible 
way. The energy being convex in fact means we can essentially close V under the energy 
viewed as a norm. Such Banach space are well studied and called Sobolev-Orlicz spaces. 

A BRIEF introduction to Sobolev-Orlicz spaces 

In this section we recall several well known theorems concerning Banach spaces 
and in particular the Sobolev-Orlicz spaces. 

One of the key uses of Sobolev-Orlicz spaces (introduced below) is to produce 
Banach spaces B where we have no control over L*' growth for p > 1 yet are still able to 
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represent B as the dual of a second Banach space, i.e. B = D* . Recall that Gel'fand's 
theorem tell us that even L^([0, 1]) cannot be realized a the dual of any Banach space, 
so would not do. 

The reason we would like B = D* is that such a relation gives us sensible 
notions of compactness in the weak topology on B. 

Theorem 13 (Alaglu's Theorem) A closed and hounded set in the norm topology is 
weak compact. 

To identify such set's we will recall on of the most well known convexity theo- 
rems, namely... 

Theorem 14 (Mazur's Theorem) A closed and convex set in B is also closed in the 
weak topology. 



An Orlicz space is an type space using a different convex function than 



= tP. In fact its good to specify the class of convex functions on [0, oo) (the Young 
functions) of use here. A young function is 



where (/>(s) satisfies 

1. (t){s) > for s > 

2. (j){s) is left continuous 

3. (/>(s) is non-decreasing on (0, oo) 

4. (j){oo) = oo. 

We would like to form a norm which behaves something like 



For proof's of these results see |11]. Now let's produce the spaces with the 



B = E* property to which these theorems will apply. 
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In fact the space of function with satisfy p^{f) < oo wiU be essential to us, and we wiU 
denote it L$. 

To actually implement this we take what you can of the inverse of (p, namely 

let 



'ip{t) = sup s. 

</.(s)<t 



Then let 



{t)= j\{s)ds 
Jo 



and call it the the Orlicz conjugate of ^{t). 

Perhaps the most important example in this context are the following conjugate 
relationships.... 







Q<t<l 



log(t) + 1 t>l 







< t < 1 



tlog{t) t>l 



m = { 



1 Q<t<l 



e*-^ t > 1 



t 0<t<l 



nt-l 



t > 1 



>. 

With these notions let 

WfW^ = sup^eU J \fix)v{x)\dA. 

Then this forms a norm on the space of measurable function with 
(call it L*). 



# < oo 



Ill 



The ^ come in two flavors the happy ones which satisfy the existence a T > 
and A; > such that 

and the sad ones which don't. Notice tlog'^(t) is happy and it conjugate e^^^ is not. 

Let be the closure of C°° in under this norm. For happy $ we have that 
L$ = L<j> = is a separable Banach space and that the norm interacts with nicely. 
For example if /„ / then p^{fn) P^if) (this is a special consequence of what is 
known as mean convergence) For sad (p we have L<j, C L$ C Eq, and that these inclusions 
are always proper, also fails to be separable and fails to even be a vector space. 

and Lq, are reflexive if and only if $ and ^ are happy. However it is alway 
the case that -L$ = (Exj/)*, the key property discussed in the first paragraph of this 
section. 

With these space the Sobolev-Orlicz spaces are easy understood. I will only 
present and need a very special case, but everything here works in complete generality 
(see the very nice §). We wih embed C°° into L'^{M) x L^{T{TM)) x L^i„g+(^-,(M) via 
I{u) = (n, Vn, An) and take its closure in the Banach norm. We arrive a Banach space 
B. Just as above even though is bad we can realize B as E* were E is the closure 
of C°°(M) in L2(M) X L'^{r{TM)) x L^t-i{M) under the same embedding. 

Its worth noting at this point that we in fact have certain obvious continuous 
inclusions of the classical Sobolev spaces. Let H^'^ be the usual Sobolev space where 
we control the norm of the first k derivatives. With these we clearly have have the 
following continuous inclusions 

c B C iJ^'^. 

This fact gives us some nice functions in B. Namely we have the Fredholm Alternative 
assuring us that 

A(i?2.2,) = {/ G l2 I / fdA = 0}; 

so for any mean zero L? function / we can construct a function \n g B such that 
A/ = 5- 
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The last fact is a certain set inclusion, namely 

In particular the mapping r/ such that r]{u) = e" takes H^'^ into for all p and H^''^ 
into H'^~^''^. I will refer to this fact as Trudinger's inequality. 

Step zero continued 

With our introduction to the needed spaces out of the way we may proceed 
with step zero by letting N be the closure of F in 5, and N the subset of this where 
esssup{/S.(f) — k) > 0. One key property of the space B in this context is that B = E*. so 
by Mazur's theorem the convex set N is closed in both the norm and weak topologies. 
The remainder of step one can be summed up in the following lemma assuring the 
continuity of E. 

Lemma 33 E is continuous on N in both the norm topology and weak topologies. 

proof: To see the norm topology case let (pn 4> implies A0„ converges to A(/> in the 
II ■ llLiog(L) norm. So A(/)„ — k will converge to A(/> — as well, since k € C°°(M). In 
particular since the L log^(L) norm is happy we have Jj^^^(A„0 — k) log^ | A„(/) — k\dA — >■ 
Jj^j{A(t)-k) log+ \A4>-k\dA. Now observe that on N PLiog+{L) = !m^^^~^) log"^(-A(/)- 
k)dA differs from Ehy a continuous and bounded function. So by the fact that the norm 
convergence implies convergence wc have from the dominated convergence theorem 
that Jj^{An4> - k) log \An(t> - k\dA ^ J {A(j) - k) log |A0 - k\dA as needed. 

The weak topology assertion follows from the fact that E is convex on N hence 
a convex function on its closure, so in the norm topology E~^{[a, b]) is closed and convex 
hence by Mazur's theorem closed in the weak topology. 

q.e.d 

Step 1. Now we would like to proceed as in the discrete case and use a 
compactness arguments to assert the existence of a function achieving the maximum. 
Since N is not quite compact in this case, we must do a little work to see that it is 
compact enough. The first thing to note is 

Lemma 34 sup^ E <Q. 
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Proof: The boundedness of ii^ on N follows form the fact that for any (j) £ B 
that Jj^j A(/) = since A</> is the limit of C°° functions with this property. So we have 
that by Jensen's inequality that 

dA 



< / {-A(j) + k)log\- A(l) + k\- 



2Trx{M) 
dA 



-2vrx(M) ■ 

This along wit the obvious fact that — Jj^j WVcpW^dA < gives us the needed bound. 

q.e.d 



Denote the finite number sup-^E as m. From the lemma 33 on the continuity 
of E in the weak topology K = E~^{[m,m + a]) is weak closed and convex. In fact K 
is weak compact. To see this it is enough by Alaglu's theorem to see that the norm is 
bounded on K. By Poincare inequality ||0||2 < C'||V(/)||2 so 



B 



<Ci(||V0||i + ||A</>||^i,g+(^)) 



<C2m(t>\\l+ PLio^+^L)i^(t> -k))<C2{E{cP)+C3)<C2{m + a + C3). 

Now just as in the discrete world we have a continuous function, E, on a 
compact set K and hence we have at least one point achieving the maximum value. 

Step 2. Now just as in the discrete case we need to control the boundary. 
Suppose a maximum occurs on the boundary at (p. Just as in the discrete case we will 
construct a direction ip and a line l{s) = + s^p is contained in N such that '^^^j^^^^ > 
for all s G (0, e) for some e > 0; hence contradicting the maximality of (p. 

(j) being on the boundary implies that 

M5 = {xeM\A^-k<5], 

has measure > for all 5 > Now since Jj^^ AcpdA = — 27rx(M) > there is 
certainly an interval [a, b] such that a > e and the set S = {x £ M \ A — k & ['^j ^]} has 
measure s > 0. For each 6 let 

, XMs XS 

J5 = • 

ms s 
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Note by the Fredholm alternative that there is a tps ^ H^'^' C B such that A{ips) = 
fs, and further more by altering this function with a constant that one can assume 
Jj^ ifjidA = 0. Using this V's £ H^'^ direction we see that for small enough s that indeed 
l{s) is in N and that 

dE{l{s)) 



[ fs-l{s)-fslog{-A{l{s)) + k)dA. 

J M 



ds 

Since log \x\ tends to — oo as x tends to zero we see that for small enough that 5 that 

/ -f5\og{-A{l{s))+k)dA 

J M 

can be made a large as we'd like. Note that in Ms we have A(Z(s)) is bounded and 
hence by the Green's function representation of l{s) we have that • l{s) is bounded 
simultaneously for all small enough 5. So indeed 5 can be chosen so "''^^^g^^'*'* > for small 
s. So boundary maxima are impossible, and we have an internal maxima. 

Step 3. So now we have that our point of maximal energy is internal. Note at 
such a maxima > A; in essential supremum so ln(A;^) = —2(j) + log(A0 — A;) € L^(M). 
Using l[s) as above we see for each / G C°°(M) such that fdA = we have 

dE{l{s)) 



[ flogik^)dA = 0. 

JM 



ds 

So log(A;^) is a constant as an function. In particular exponentiating we see that A;<^ 
is a constant. In fact note 

and by Trudinger's inequality e^'^ G so <^ G i?^'^ by elliptic regularity. So e^^ G -H"^'^ 
and by elliptic regularity again ^ G Continuing this ^ and hence are in fact in 

all H^'^ and hence by the Sobolev embedding theorem in C°°. 

Step 4. Now we can easily get uniqueness in the x(M) < case, exactly as in 
the discrete case. Take two now C°°(M) solution ^\ and ^2 and note form the second 
Prech'et derivative that the line Z(s) in C°° connecting them satisfies 

for all s with l{s) in N. So as in the discrete case the point of maximum energy is 
unique. 



Chapter 5 

Spheres and Tori 



This chapter is dedicated to exploring the ideas of the previous chapters in the 
cases of x(-^) ^ 0- The cases of primary interest are the torus and sphere cases. As in 
the xi^) < case the fundamental object needing exploration is the class of polyhedra 
related to a triangulations, and the volume energy associated to this class. In both cases 
there is an intrinsic class of such polyhedra but the story presented here in the xi^) < 
case experiences difficulties. The nature of the difficulties is very different for the torus 
and the sphere. 

In the toroidal case the class of polyhedra needed has been studied in the 
literature in |^] and [14|, and the issues in chapter two have already been essentially 
dealt with. I will remind the reader of the issues involved in section ^.1| , and for now 
only highlight the differences with the < case. In the torus case the energy 

is significantly easier to deal with and analogs to theorem |^ and theorem ^ exist and 
can be proved with the same methods. However the story is different in the important 
respect that the linear part of the problem sends one immediately into the land of zero 
curvature world, a phenomena which occurs in both the discrete and continuous cases. 
So in the continuous metric case there is no non-linear metric story at all. This case 
in fact demonstrates the important fact that perhaps a better continuous analog of the 
discrete world would be a connection or an affine like structure. This because, as we 
shall see, the discrete uniformization does not produce flat structures, but rather affine 
structures. In any case the ideas of chapter three breaks down at a rather fundamental 
level. 
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Figure 5.1: The Vector 



The spherical case is quite the opposite, and on some level the ideas don't 
break down at all. As in the xi^) ^ cases there is an appropriate class of "intrinsic" 
polyhedra, and the one gets a natural volume energy which has critical points precisely at 
uniform structures. Unlike in the toroidal case the randomization goes through perfectly 
to produce an energy on the positive curvature metrics. Every thing looks good except 
now the energy is no longer nice at all. Analogs the theorem and ^ exist (see |^5|) but 
appear much harder to prove with the intrinsic methods used in the xi^) < case. In 
particular the energy fails to be convex and the boundary behavior becomes very difficult 
to control. In fact both the discrete and continuous proofs of the previous chapter fail 
in fundamental ways. Of course intrinsic proofs may still exist, but as is often the case 
with spherical uniformization (see and 1^]) will involved significantly more drama. 

5.1 The Toroidal Case 

In the discrete toroidal world one must first replace N with angle systems which 
have zero rather than negative curvature. This of course gives us a significantly smaller 
space of angles and in particular the conformal deformations preserve not only the 9^ 
but also this curvature condition, and are the span of the Wy vectors in figure |5.1| over 
all the vertices. 

To construct the polyhedra first view the Euclidean plane as the boundary at 
infinity of hyperbolic space in the upper-half space model. The polyhedra are now con- 
structed by taking the union of the ideal simplexes over each of the Euclidean triangles. 
Note the energy formula in the xi^) < agrees with this construction and becomes 
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twice the sum of the vohimes of all the ideal tetrahedra over a fundamental domain. As 
before the volume is simply the sum of the volumes of the ideal tetrahedra corresponding 
to the individual Euclidean triangle angles in the angle system. In other words 

teP 

where if = {A, B, C} we have 

Vtix)=AiA)+A{B)+AiC). 

This energy remains convex and boundary controllable and all of chapter two 
carries over with the most interesting point being why at the critical points of E fit 
together. 

Observation 6 At a critical point of E the triangles fit together to form an affine struc- 
ture on the torus. 

Proof: As far as I'm aware the idea in this proof has its origin in Bragger's Q. 
The above formula tells us that from the formula for the Lobacevskii function we have 
dE"" = Y.a,^wEi{x)a' with 

E,{x) = - In (sin(A^)) 
So at a critical point with n faces ti in its flower at v and angles labeled A\. we 

have 

= dExiWv) = m — ToT TTTT T-TT 

^ ^ Vsin(^^) sin(^i) sin(^L 
Now lets attempt to fix our edge lengths. We will denote the edge length opposite to 
^!|_ as a!|_. Each triangle can be scaled with its angles preserved since we are in the 
Euclidean plane. So fix the size of ti. Now scale t2 so that = a\. Continue this until 
the size of a" has been fixed. 

Now from the law of sins = dE{wv) gives us 

2 1 1 ' 

al a_ a_ 

So the entire flower fits together. 

Now we have open sets which are affine related in overlaps, so an affine structure. 

q.e.d 

With this observation the other ideas essentially work out in the same way. 
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Figure 5.2: A Hyperbolic Tetrahedron 

5.2 The Spherical Case 

The initial discrete set up for the sphere is identical to the xi^) < case 
except the use of positive rather than negative curvature. The class polyhedra is very 
simple to construct. Given a triangulation of the sphere take the convex hull of the 
vertices. The volume formula turns out to once again agree with the negative curvature 
case. This is quite a nice fact. To see it view the sphere at infinity in the ball model 
from the origin. Note form this view point the angles you see in the ideal polyhedra are 
precisely the angles in the triangulation as the sphere understands them. Now for each 
triangle on the sphere at infinity with = {A,B,C} form the three ideal vertexes 

tetrahedra by taking the convex hull of the these vertices at infinity with the origin, see 
figure 



Fact 7 The volume of the above tetrahedra Vt{x) is given by 

2Vtix) = AiA) + AiB) + AiC) ^^/vr-^l-S-C 



^ 'tt + A-B-C\ ^{'k + B-A-C\ ^f-K + C-A-B 
+A + A + A ' 



2 J \ 2 J \ 2 

Proof: To see this extend the geodesies in the tetrahedra and take the convex hull of 
this arrangement. We get an ideal octahedron. Using the three new points at infinity 
and the origin note we have in this octahedron a symmetric copy of our original three 
ideal vertexed tetrahedra. 

Each edge Cj of the octahedron corresponds to an ideal tetrahedra Tj, see figure 
|5.3| where three particularly relevant tetrahedra are labeled. Now simply note that by 
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using both copies of the three ideal vertexed tetrahedra that the needed volume can be 
expressed as half of Vol{Ti) — Vol{T2) + Vol{T3). Now the fact that an ideal tetrahedra 
has its volume given by summing the Lobacevskii function over the angles meeting at a 
vertex gives the needed formula. q.e.d 
Note by spherical geometry that 

, , cos(C) cos(-B) + cosf^d) 
sm(i?) sm(Gj 

so the same exact computation as in the x{^) < case tells us a critical point is uniform. 
However convexity and boundary control are both lost, and since we may use the same 



formula as in section 2.1.4 this is easy to see. The random computation goes through 



as in the x{^) < case, and we can see immediately from formula L6 that we lose 
convexity. 

It is worth noting that we should have expected problems, at least with unique- 
ness. Namely there is a sort of Gauge group sitting around, and in the metric world it 
corresponds to the fact there are conformal transformation of the standard sphere which 
fail to be isometrics. In fact given a uniform metric there is a three dimension space of 
distinct (j) which remain constant curvature. This is fun to witness in terms of triangula- 
tions where it indicates that we may expect distinct sets of angles to be be conformally 
equivalent and still fit together. To see it fix the vertexes of a triangulation and move 
what we view as the origin of hyperbolic space in the above construction away form the 
balls origin in the model. Then moving it back to the model's origin produces a topolog- 
ically equivalent triangulation with distinct angles which clearly fits together. Although 
this observation is unfairly mixing our two notion of conformal change it is still indicates 
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that uniqueness of a uniform structures in the discrete world should not be common. 



Chapter 6 

Appendix: A Less Pleasant Proof 
of Lemma ^ 



Here I will present an alternate proof of lemma |l|, and arrive at the slightly 
stronger condition of needing circles to only be on disks of radius less than min{^,T} 
rather than on circles of radius less than | as in section 3.1.2. 

The trick to this proof of lemma |^ is to understand the curves satisfying d{p, z) — 
d{q,z) = 0; with d{p,z) with less than mm{|,r}. This because any point on such a 
curve corresponds to the center of a circle going through both p and q, and if a triple 
{p, q, o} lives on a circle then the corresponding curves for each pair in the triple must 
intersect at the the point corresponding to the center of this circle. 

Before getting started there are a few basics pieces of notation convenient to 
introduce here: if d{p, q) < i call jp^q the shortest length geodesic segment between p 



and q (it is well defined by lemma 12), and let jp^q C jp^q be the connected component of 
the geodesic contained in any set we happen to be exploring with 'jp^q in it. For example 
relative to Bi{p), 7p_q is the geodesic splitting Bi{p) into its two distinct "sides" (simply 
look in normal coordinates). Denote as {%^q}'^ the two components of {7p,q} —'lp,q- We 
will also find it useful to name the midpoint of ^p^q - called it m. (see figure |6.1| for 
periodic notation reminders). For the remainder of this section denote d{p,z) as Dp{z), 
since the differential is usually represented with a d. 

To get started it is in fact useful to consider the more general curve of the type 
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Figure 6.1: The Notation 



Dp — Dq = c, were c G R. By the implicit function theorem when d{Dp — Dp) 7^ and 
is defined, the solution to this equation is locally a curve z{t) with z{t) 7^ 0; where z{t) 
satisfies d{Dp — Dp){z{t)) = 0. In fact any such curve will parameterize a solution. It is 
useful to rephrase this one form business in terms of its dual object the gradient. Recall 
the gradient of a function / is the unique vector field V/ satisfying < V/,?; >= df{v) 
for all V at every point. So we may rewrite our differential relation d(Dp — Dp){z{t)) = 
as < VDp - VDg,z >= 0. 

To get a grip on this differential relation it is first useful to spend a moment 
contemplating V-Dp. 

Sub-lemma 2 (The Distance Gradient Sub-lemma) Let d{p, q) < ^ then: 

1. In Bi_(p) — {p} we have VDp is unit length with its integral curves the unit speed 
geodesies. 

2. In Bi{p)P\Bi{q) we have VDp = —VDq on %,q, VDp = VDg on {%,q}^ ; and 
outside 7p^g we have VDp 7^ cVDg for any c G R. 

Proof : For the first part note the distance function's level sets are the spheres, so 
by Gauss's lemma its integral curves are some re-parametcrizations of geodesies from 
p. Now observe in geodesic polar coordinates that the unit speed geodesies 7 satisfies 
7 = G4|:),and 

< VDp,^ >= dDp (g. (|;)) = G*{dDp) = dr (|:) = 1. 
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So indeed the integral curve of VDp are precisely the unit speed geodesies. 

The equalities in the second part follows immediately from the first part and 
the fact 7p^g is a geodesic. 

To prove the last piece of the second part assume at some z E Si (p) Pi _» (q) 

2 2 

we have VDp = cVDq. First note from the above we have c = ±1. There are two 
cases, first we'll deal with c = 1. Since the geodesies satisfy a second order O.D.E they 
are uniquely determined by their position and tangent vector, so when c = 1 we have 
both the geodesic from p and the geodesic form q are the same curves. Without loss of 
generality p is further away than q and this point lies along the same minimal length 
geodesic (of length less than |) which connects p and q, i.e. 7p^g. In the case c = — 1 we 
can follow the geodesic form p to the point and then from the point back to q forming 



a geodesic of length less than i - which then must by lemma 12 be the unique such one, 
i.e. 7p,q. 

q.e.d 

Back to the relation < VDp — VDq,z >= 0. The first observation is that we 
can express a solution of this relation via a vector field. Using ^ rotation field Q (from 
section 1.1) we see that the solution to the differential relation < VDp — VDq, z >= are 
re-parameterizations of integral curves of the vector field @{VDp — VDq). Fortunately, 
as with VDp, we can say quite a bit about this vector field. We are most interested in 
its integral curve corresponding to Dp — Dq = 0. Let Cp,q{t) be the connected component 
of the integral curve of Q^VUp — 'VDq) passing through m in any set of interest to us; 
and assume its parameterization satisfies Cp^g(O) = m. 

Lemma 35 (The Distance Difference Gradient Lemma) Assume d{p, g) < | then: 

1. In Bj^{p) we have that Cp^q{t) is the unique component of Dp — Dq = 0, and m is 
the unique point of ^p^q on Cp^q{t). 

2. In Bi.{p)[jBi{q) we have d{cp^q{t),p) and d{cp^q{t),q) strictly increase as the pa- 
rameter \t\ increases. 

3. No geodesic from p in Bi{p) is tangent to Cp^q{t) or cuts trough Cp„{t) twice on the 

4 ' ' 

same side of^p^q{t). 
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4- In Bi{o) — o it can never he the case that VDq — VDp = c(VZ)o — VDq) ^ at a 
point on Cp^q{t). 

Proof : First for the uniqueness of m: suppose a point I ^ m\s 7pg. Then / is within | 
of p- hence d{p, I) is determined by the length of the segment of "fp^q form p to /, similarly 
for q (using |). Now note that as we move from ra toward, say, p that Dp decreases 
while Dq increase - so Dp — Dq ^ d a,t another point of ^p^q. When on {7p,g}^ say above 
p the segment of ^p^q from q to I va. fact covers the shorter segment from p to I - forcing 
Dp — Dq ^ once again. So / cannot satisfy Dp — Dq = 0, forcing m to indeed be the 
unique point of 7p^g on Dp — Dq = 0. 

The remainder of the first part and the second part are intimately related. To 
see why we first look at the component of Cpq{t) in Bi_[p) and note any component of 
Dp — Dq = would have to have a point closest to p. This closest point is tangent 
to a sphere emanating form p. The same sort of phenomena must take place for the 
distance function to have a critical point; namely if a point z along any integral curve of 
QCVDp — VDq) is a critical point of the distance function D{p, •) then either VDp = VDq 
or a circle is tangent to the solution curve. In the tangent case @{VDp — VDq) = cQS/Dp, 
or rather VDp — VDq = cVDp] so both these situation we have forced the case VDp = 
cVDq. The triangle inequality tells Bi_{p) C Bj_{p) Bj_{q) so we may use sub-lemma |2| 

4 2 2 

to note that the point where this occurs is on jp^q] but from above to be on Dp — Dq = 
and jp^q means you must be exactly m. So we have both that every component of 
Dp — Dq = in Bi_ (p) contains m , and that the distance to p parameterized by t can 

4 

have no critical points except at m (similarly for q). 

For the third part note that if a geodesic cuts twice on the same side of 7p,g(i) 
then some other geodesic must be tangent (via the mean value theorem in geodesic polar 
coordinates - see the see figure |OD . So we are reduced to the tangent case. Let the point 
where this tangency occurs be called z. This tangency implies @{VDp — VDq) = cVDp, 
or rather QVDq = QVDp — cVDp. But VDp and QVDp are orthogonal, with QV Dp 
also of unit length (by sub-lemma |2|) - so taking norms we have 1 = V 1 + ; forcing 
c = 0. So at such a point @{VDp — VDq) = or VDp = VDq. Once again using 
Bi{p) C Bi{p)f\B2_{q) the above sub-lemma applies and forces z onto {^.g}*^; and 
now part one kicks in to eliminate this possibility of Cp^q hitting {7^,5}*^. 
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Figure 6.2: The Hunt for the Tangent 



To prove the last part we will see first if VDq — VDp = c(VDo — VDg) 7^ 
at z, then c = 1 and VDp = VL'g at z. This because the VDj are all unit vectors (by 
sub-lemma ^ part one); and among unit vectors v and w ^ v the direction of w — w is 
uniquely determined by w (since we are re-parameterizing the radius one circle at located 
a v). Recall from the above sub-lemma that V-Dp = VDg at z implies z G {7p,g}'^. As 
above, no point of Cp^q{t) can be on {7p,g}*^, as needed. 

q.e.d 

The last set of issues concerns not the Cp^q curves themselves but the home of 
the curves. Suppose we have three points p, q and o on the boundary of circles of radius 
r less than |; then using jp^o and 7g^o we can split B^{o) into four cones (as in figure 
|6.3| ). Let the forward length cone (FLC) be the cone containing 7p^g , the backward 
length cone (BLC) be the forward length cone's mirror image (see figure |6.3| and the 
following lemma to see this is well defined). One nice thing about this cone notion is 
that it eliminates certain regions where the center of a small circle might have wanted 
to live - and it is at this point where we find the only place in the argument where using 
the strong convexity radius is necessary. We shall explore these facts (and clear up an 



issue in the introduction) with the following corollary to lemma 12 in section 1.2: 



Lemma 36 (The Length Cone Lemma) Given three points p, q and o on the bound- 
ary of a circle of radius r less than | centered at c: 

1. The forward length cone is well defined. 

2. The intersection of the three forward length cones is a triangle contained in Bi (c) ; 
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Figure 6.3: The Length Cones 

which, incidentally, can be parameterized in a way compatible with the abstract 
gluing to form |i^'{pi,...,p„}| (see figure |6^ . 

3. The Bi_{c) is dissected by the curves {7p,q, 7p,o) 7g,o} 'i-nto seven regions - the unique 



FLC, three distinct BLCs, and the regions separating the BLCs (see in figure 6.4 j. 



4- If in addition the balls radius is less than r, then circle's center cannot be in any 
of the backward length cones, and the simplex is contained in Br{c). 

Proof: In the course of the this proof I use the triangle inequahty in the sense of 
p G Ba{q) then Bi,{p) C Ba+b{q)- I will simply call it the triangle inequality. 

First we attack well definedness of the FLC. Since the radius of the circle is less 
than I we have by the triangle inequality p and q are in i?2i (o), and that d{p,q) < ^. 
Now note 7^.^ is contained in Bm (o) - since if it left Bm (o) then by the triangle inequality 

6 6 

it would be a curve from p to q leaving Bi{p) and then afterwards entering Bi{q), so 

6 6 

would have length greater than ^ ; contradicting < ^ . We would now like this 

geodesic segment to stays in one cone section. If it did not then it would cut a bounding 
geodesic, ^o,p or 7^^^, in at least two points on the cone's boundary. Now the diameter 
geodesies are of length less than i, so between ci and C2 we now have contradictory 
distinct geodesies of length less than i connecting them. So the cone is well defined. 

To get started on the second part first simply parameterize the needed simplex. 
Note we may choose 61 and 62 such that 7o,p(t) = expo{tv{6i)) and ^o,q{t) = expo(tv{92)) 
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bound the FLC, and by the first argument ^p^q C Bm (o) so any point on ^p^q has its short- 

6 

est length geodesies from o described by expo{tv{6)) with ^ € 62]. The fact there are 
no double intersections of short geodesies now demonstrates that each expo{tv[9)) with 
6 € [0i^92\ hits exactly one point on 7^^^, so we can parameterize ^p^q as expo[r{9)v{9)) 
- with r{9) continuous. Now we have a homeomorphism (once again since there are not 
double intersection between small geodesies) of the region trapped between the geodesies 
with the wedge 9 G [^1,^2] and r < r{9). This wedge is clearly homeomorphic to the 
disk. However to parameterize the triangle we would like a parameterization which is 
compatible with the abstract gluing, allowing R : — M to be continuous. To 

accomplish this first parameterize the edges with the unit interval using a d(p, q) speed 
parameterization of the geodesic. Now for each triangle in the surface we will map on 
an equilateral triangle by first mapping the boundary of our triangle to the boundary of 
our simplex M such that the linear unit length sides are mapped onto the geodesies at 
a speed of d{p, q). We have just found that the image is in fact itself the boundary of a 
disk, so this mapping can be extended from a homeomorphism of the equilateral trian- 
gle's boundary to a homeomorphism of the equilateral triangle onto the triangle in M. 
So by its very construction the mapping on the boundary of the equilateral triangle is 
afHne related to the mapping of the interval, and we indeed have the needed compatible 
par ameterizations . 

We also need for this second part that this is the same simplex for all three of 
the points, but by the very definition of the FLC with respect to a point being the cone 
section containing the opposite 7 we have that the FLC must contain this simplex. So 
the simplex is in the intersection - and precisely the intersection since the respective 7^,^ 
curves all cut their respective Bi{p) balls in to two distinct sides - and by the triangle 
inequality this entire discussion is taking place in any one of the radius i balls. 

Lastly for the second part we need the simplex is contained in Bi (c) ball. To 
see this note from above that every point, k, in the simplex is on a geodesic, from say 
o, of length less than ^ which hits jp^q. If k is within ^ of o then k is with in | of c. If 
not it is within g of 7^,0 • The geodesic segment 7^,0 has length less than ^, so the point 
where the geodesic from k hits 7^,0 is within | of either p or q; so k is once again with 
in ^ of c - as needed. 

For the third part note that by the triangle inequality a geodesic between two 
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Figure 6.4: The Seven Regions 
points in Bi{c) is contained in the ball must be of length less than i (use the sense 

2 

of the triangle inequality introducing this proof at a point were the geodesic first hits 
Si^{c)). Now extend all the geodesic segments bounding the simplex out to the boundary 

2 

in B±{c); and note, having length less than i these segments can intersect each other 

2 

only once - at the vertices of the simplex (by part 2 above the simplex is in there). So 
noting the definitions of all of our objects we indeed end up decomposing U into the 
seven components pictured in the figure. 

For the final part, the fact that the simplex is contained in Br{c) is a trivial 
consequences of the 7 curves having interior in this ball and never intersecting each 
other with in the ball. To handle the center's placement, assume the contrary that the 
center, c, is in the backward length cone with respect to say o. Then take the shortest 
length geodesic from c to o and continue it through the forward length cone until it 
hits 7p_g at h. Note this is a geodesic of length less than ^, so in particular it is the 
unique shortest length geodesic between these points. By assumption p and q are on the 
boundary of a disk at c with radius less than r, so by strong convexity jp^g C Br{c). In 
particular the shortest length geodesic to h never leaves Br{c), contradicting the fact it 
hit o before h. So the center must not be in the backward length cone. 

q.e.d 

We are finally in a position to prove the small disk uniqueness theorem form 
the introduction (theorem |l]) . 

Proof of Theorem |^: As above call the points in the triple p, q and o, and 
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recall 6 < min{|,r}. The whole discussion takes place with Bi_(c), since d{cp^q{t),c) < 
d{cp^g{t),p) + d{c,p) < I- In particular the whole discussion takes place at points where 
we can use lemmas ^ and Note a circle of radius less than | passing through 
all three of these points corresponds to an intersection of Cp^g and Cp^o (by part one of 
lemma ^5|). Further note if Cp^q and Cp^o intersect at z then so does Cq^o, since then 
d{z,p) = d{z,q) = d{z,o). 

At these intersections we would like that the curves must cut through each 
other; i.e. not be tangent. If a pair were tangent then there are two possibilities: one 
where neither of the tangent curves is stopped when thought of as a parameterized 
integral curve, and the other when one has stopped {cp^q = 0). If neither has stopped 
then (up to index permutation) Q{S/Dp — VDg) = c@{VDp — VDq) ^ 0, or rather 
{VDp — VDq) = c{WDp — VDo) 7^ 0. This case implies by lemma |2| that we are on 
{Tg.o}'^- If one is stopped then (once again up to index permutation) VDp — VDq = 0. 
So by sub-lemma ^ the point of intersection is on {7^,^}'^. Be warned - I will refer to this 
same case analysis later in the argument. To finish off this impossibility argument simply 
note that all the curves intersect at this point so (up to index permutation) {7p,g}'^ and 
Cp,g would intersect. But by part one of lemma ^this point must then be TUp^q which is 
not on {%,q}'^ ■ So the intersection must include no tangencies. 

Now suppose we had two distinct circles of radius less then | passing through 
these triples. This non-uniqueness of circles would correspond to at least a pair of 
intersections of Cp^q, Cp^o-, and Cq^o - at zi and Z2. Note without loss of generality there 
are no further intersections of any pair of c curves between zi and Z2; since intersections 
occur only as triples at isolated possible distances by the above no tangency result. Hence 



we are forced to have a picture something very much like figure S.5, at least in the sense 
that we have one of the curves trapped between the other two; without loss of generality 
assume this curve is Cq^o- 

Note @{VDq - VDo) = e{VDp - VDo) + G(VD, - VDp) , so G(VL>p - VDo) 
and QiVDq — VDp) are always on opposite sides of QiVDp — VDo); in the sense that 
the VDq — VDo component of these vectors must have opposite signs. Now since the 
curves must cut through each other at zi and Z2 and no where in between, the side of 
Cq^o which Q{VDp — VDo) and Q{VDq — VDp) live on must switch roles while one moves 
along Cq^o from zi to Z2- In order for these roles to switch and the vectors to always be 
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z 

2 



Figure 6.5: The First Possibility 

on opposite sides, we are forced (as with our no tangency result above) to have at some 
point along Cq^o either @{VDp - VDo) = cQ{VDp - VDg) / 0, VDp - VDo = 0, or 
V-Dp — VDq = 0. Immediately, the last part of lemma |3^ kicks in eliminating the first 
case. 

The remaining case (up to index permutation) is were VDp = VDq - which 
means this point lies on {jp^o}^ (by the sub lemma). We will be justifying each of the 
two pictured cases - depending on which component of {%,oY the point lies (figure |6.5| 
and figure |6.6| ) . The idea in each case is this same - we will use this condition to force 
one of the centers Zi into a backward length cone with respect to either o or p, and then 
use the previous lemma to observe this is impossible. 

Observe by the first and third parts of lemma ^ that Cq^o{t) cannot be tangent 
to or cut twice any geodesic in the picture; and Cq^oif) cannot sneak around a geodesic 
in the picture (since the decomposition in the third part of lemma |6.3| has each of the 
geodesies heading all the way out to Sj^{c)). So in the first case the Cq^o is stuck to the 
right of 7g^o and must cross 7p^o to the right and can never cross it again. So one of the 
Zi is stuck in the backwards length cone of o as needed. The other case is similar to this 
except we must justify why we must come from above and cross into the backward length 
cone with respect to p. This is because to come from below would mean to cut across 
7p^o first to get there producing a double cut of Cq^o by 7p,o to left of 7g^oi contradicting 
part 3 of lemma |3^. So indeed once again we are forced into the contradictory setting 
of having a center in a backward length, this time with respect to p. 

So the points zi and Z2 cannot simultaneously exist - and there is no non- 
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Figure 6.6: The Second Possibility 



uniqueness among small circles. 



q.e.d 
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